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SUMMARY 

The  most  popular  large  commercial  building  chilled  water  system  today  is  the  primary-secondary 
design.  This  system  separates  the  generation  zone  (chiller/evaporator)  of  a  chilled  water  system  from  the 
distribution  zone  (air  handlers).  The  generation  (primary)  loop  maintains  a  constant  chilled  water  flow  rate. 
The  distribution  (secondary)  loop  varies  flow  rate  in  response  to  cooling  load  fluctuations.  This  variable 
flow  rate  in  the  secondary  loop  permits  energy  savings  over  older,  constant  flow,  single  loop  systems 
(Rishel,  1996).  Recently,  new  chiller  controls  have  been  implemented  that  allow  variable  chilled  water 
flow  rate  in  chillers.  As  a  result,  a  new  chilled  water  system  design  has  become  available.  This  design  is  a 
single  loop  system  with  a  single  set  of  pumps.  Like  the  primary-secondary  system  the  single  loop  design 
has  two-way  valves  in  air  handlers.  This  variable  flow  rate,  primary  system  also  eliminates  the  primary¬ 
secondary  bypass  loop.  Existing  primary-secondary  systems  can  be  converted  to  this  design  by  simply 
valving  off  the  bypass  leg. 

The  purpose  of  this  study  was  to  compare  the  performance  characteristics  of  the  primary- 
secondary  system  with  the  single  loop,  variable  flow  rate  design.  To  accomplish  this,  individual  system 
components  were  modeled.  These  component  models  were  compared  with  manufacturer’s  data  in  order  to 
validate  their  accuracy.  The  modeled  components  include;  air  handler,  chiller,  pump  and  cooling  tower. 
These  component  models  were  combined  to  form  primary-secondary  and  single  loop  system  models. 

Comparisons  of  these  two  systems  were  made  to  determine  relative  system  efficiency  and  cooling 
capacity.  Different  designs  and  cooling  load  conditions  were  applied  to  the  primary-secondary  and  single 
loop  systems.  These  situations  include: 

1)  Single  and  multiple  chiller  systems 

2)  Constant  and  variable  speed  pumping  schemes 

3)  Various  load  profiles 

4)  Under  and  oversized  air  handlers 


5)  Chilled  water  coil  fouling 
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Each  comparison  showed  the  single  loop  system  to  have  greater  efficiency  and  cooling  capacity  than  the 
primary-secondary  design. 

Different  control  strategies  and  designs  were  employed  in  an  effort  to  improve  primary-secondary 
system  performance.  None  of  these  changes  allowed  the  primary-secondary  system  performance  to  rise  to 
that  of  the  single  loop,  variable  flow  rate  design.  As  a  result,  this  investigation  concludes  that  the  variable 
flow  rate,  primary  (single  loop),  chilled  water  design  is  advantageous  compared  to  the  popular  primary¬ 


secondary  system. 


13 


CHAPTER  I 
INTRODUCTION 
Background 

In  the  recent  past,  due  to  rising  energy  costs,  manufacturers  have  invested  great  time  and  money  to 
increase  HVAC  component  efficiency.  However,  it  is  impossible  to  predict  entire  system  effects  by 
studying  only  a  single  component.  Certainly,  increasing  a  component’s  efficiency  will  affect  overall 
system  performance.  For  example,  a  chiller’s  efficiency  is  increased  by  adding  heat  transfer  surface  area  or 
augmentation  to  the  condenser  and  evaporator.  This  raises  the  chiller  efficiency  by  increasing  heat 
exchanger  UA.  However,  the  addition  of  area  or  augmentation  can  increase  the  heat  exchanger  pressure 
drop,  and  with  it  the  pump  power  requirement.  Without  analyzing  the  entire  system,  it  is  impossible  to 
determine  the  effect  that  this  change  has  on  system  performance.  This  study  compares  the  system 
performance  of  two  basic  chilled  water  piping  schemes  in  order  to  determine  their  individual  advantages 
and  disadvantages. 

1 .  System  Descriptions 

In  the  early  1970’s,  variable  speed  pumping  was  rarely  used  in  HVAC  applications.  Constant  speed 
pumps  were  used  throughout  chilled  water  loops.  This  resulted  in  inefficient  part  load  operation,  because 
the  constant  speed  pumping  scheme  maintains  a  large  flow  rate  and  pressure  drop  at  all  loads.  In  the  late 
1970’s,  the  primary-secondary  chilled  water  system  was  developed  to  take  advantage  of  variable  speed 
pumping.  This  system  maintains  constant  flow  rate  in  the  primary  loop.  It  has  a  bypass  loop  that 
uncouples  the  primary  loop  from  the  distribution  loop.  This  allows  for  variable  speed  in  a  distribution 
(secondary)  loop.  These  changes  resulted  from  increased  energy  costs,  making  the  larger  first  cost  of  a 
more  efficient  system  worth  while.  Over  the  past  20  years,  numerous  problems  with  the  primary-secondary 
arrangement  have  been  discovered.  These  problems  have  prompted  the  study  of  different  chilled  water 
arrangements.  One  of  the  proposed  arrangements  is  a  system  that  varies  the  water  flow  rate  through 
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chillers  as  well  as  the  load.  It  uses  variable  speed  primary  pumps  on  chillers  and  two  way  valves  in  all  air 
handlers.  These  three  systems  are  described  in  greater  detail  below, 
a)  Constant-flow  rate  primary  systems 

The  constant- flow  rate  primary  chilled  water  system  maintains  constant  flow  throughout  the  entire 
system.  All  air  handlers  have  three-way  modulating  valves  that  provide  constant  flow.  The  temperature 
difference  across  the  water-side  of  an  air  handler  changes  with  load.  This  design  was  common  when 
energy  was  relatively  cheap.  At  that  time,  variable  speed  drives  were  not  used  in  HVAC  applications  and 
thus  constant  flow  rate  was  necessary.  Figure  1  shows  the  layout  of  the  constant-flow  rate  primary  system. 


Constant  Speed 
Pump 


Chilled  water 
coils 


Balancing 

valves 

Three-way 

valves 


Figure  1.  Constant-volume  primary  system 

The  primary  constant- volume  system  had  several  wasteful  conditions  that  are  no  longer  acceptable 
(Rishel  1996).  For  example, 

1)  Constant  flow  means  the  system  maintains  the  same  head  in  all  situations.  This  translates 
to  identical  pump  power  requirements  on  the  hottest  and  coolest  days. 

2)  Balancing  valves  were  used  to  increase  the  head  in  coils  near  the  chiller  in  order  to 
maintain  adequate  flow  in  distant  coils. 

3)  As  the  system  load  decreases  the  AT  across  each  chiller  is  reduced,  but  the  chillers  must 
remain  on-line  to  maintain  system  flow. 


b)  Primary-secondary  systems 
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Currently,  the  most  popular  chilled  water  design  is  the  primary-secondary  system.  It  separates  the 
chiller  zone  from  distribution  zone,  thus  taking  advantage  of  variable  speed  pumping.  It  uses  two-way 
valves  in  the  air  handlers  that  modulate  secondary  loop  flow  rate  with  load  requirements.  A  bypass  line  is 
included  between  the  chilled  water  supply  and  return  to  permit  constant  flow  through  the  chillers  (Rishel, 
1996). 


Variable  speed 


Constant  speed 


pump 

Figure  2.  Primary-secondary  system 

While  the  primary-secondary  system  takes  advantage  of  the  reduced  pumping  cost  offered  by 
variable  speed  secondary  loop  pumping,  it  still  has  drawbacks. 

1)  The  increased  first  cost  of  acquiring  primary  and  secondary  pumps. 

2)  Constant  evaporator  water-side  pressure  drop. 

3)  Mixing  from  supply  to  return  across  the  bypass  in  part  load  situations.  This  creates  a  low 
temperature  differential  in  the  primary  loop. 

4)  Mixing  from  the  return  to  supply  while  operating  above  design  conditions.  This  raises 
the  coil’s  chilled  water  inlet  temperature. 

5)  Chillers  need  to  stay  on-line  in  order  to  prevent  detrimental  mixing  when  low  system  AT 
problem  exists.  These  problems  can  be  caused  by  air  handler  fouling,  poorly  controlled 
air  handler  set  points  and  a  variety  of  other  causes. 

c)  Variable-flow  rate  primary  (Single  loop)  systems 
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Variable-flow  rate  primary  systems  take  advantage  of  modem  control  systems  which  allow 
variable  chilled  water  flow  through  the  chiller’s  evaporator.  This  is  done  by  removing  the  bypass  and  using 
variable  speed  drives  with  primary  pumps.  Two-way  valves  modulate  the  entire  system  flow.  This 
variable  flow  allows  for  reduced  pumping  power  due  to  a  smaller  evaporator  head  at  part  load.  The  system 
also  prevents  mixing  between  the  chilled  water  supply  and  return  lines.  Despite  these  advantages  this 
system  is  not  widely  accepted  primarily  because  it  presents  two  problems. 

1)  The  possibility  of  laminar  water  flow  in  the  evaporator  at  small  loads  and  low  flow  rates. 
This  situation  is  possible  in  single  chiller  facilities  where  one  chiller  is  required  to  meet 
all  loads.  Laminar  evaporator  flow  is  far  less  likely  in  multiple  chiller  plants. 

2)  Difficulty  starting  up  lag  chillers. 


Variable  speed 
pump 


Figure  3.  Single  loop  system 


Control  strategies  are  available  to  rectify  these  problems.  Chilled  water  reset  is  an  option  which 
prevents  laminar  flow  in  the  evaporator.  At  low  system  loads  the  chiller  set  point  is  raised.  This  change 
has  two  results.  First,  the  system  flow  is  increased.  Air  handlers  require  more  flow  to  increase  the  UA, 
because  the  logarithmic  mean  temperature  difference  (LMTD)  between  the  air  and  water-sides  is  smaller. 
Secondly,  the  chiller’s  power  requirement  is  reduced  due  to  a  smaller  compressor  lift.  One  can  see  that  this 
technique  is  a  trade-off.  Its  use  increases  pumping  while  reducing  chiller  power  requirements.  Under  most 
situations  this  strategy  increases  efficiency  and  as  such  it  is  common  practice  by  the  HVAC  community 
(Erth,  1987).  Therefore,  not  only  does  it  prevent  laminar  flow  but  saves  energy  as  well. 
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Some  facilities  require  a  constant  chilled  water  supply  temperatures.  In  this  case  chilled  water 
reset  is  not  an  option.  Although  variable  flow  may  be  attained  by  placing  a  control  valve  in  the  bypass  line 
that  opens  to  prevent  laminar  evaporator  flow  (Coad,  1998). 

Centrifugal  chillers  are  protected  from  harm  by  a  number  of  safeties.  One  such  device  is  an 
evaporator  limit  safety  that  takes  a  chiller  off  line  if  the  refrigerant  temperature  falls  below  30°F.  This 
situation  often  occurs  as  a  result  of  an  abrupt  chiller  load  variation.  For  example,  a  multiple-chiller,  single 
loop  system  may  be  operating  with  a  single  chiller  at  a  relatively  large  load.  The  load  increases  until  the 
lag  chiller  is  brought  on  line.  The  lead  chiller  flow  rate  and  load  are  rapidly  reduced  as  the  lag  chiller  is 
energized.  The  refrigerant  temperature  in  the  evaporator  falls  below  30®F  because  the  compressor  hasn’t 
had  time  to  respond  to  the  reduced  load.  The  lead  chiller  drops  off  line  (Eppelheimer,  1996).  In  older 
chillers  the  evaporator  limit  safety  is  a  mechanical  thermostat  that  senses  the  instantaneous  refrigerant 
temperature.  Modem  low  evaporator  temperature  control  uses  more  sophisticated,  solid-state  technology. 
Modem  chillers  have  a  microprocessor  based  control  that  averages  over  a  50  second  interval.  This  control 
allows  the  refrigerant  temperature  to  momentarily  fall  below  the  minimum  set  point  without  tripping  the 
chiller  off  line.  The  delay  gives  the  compressor  time  to  adjust  to  abrupt  load  variations  that  accompany 
chiller  scheduling  (Kirsner,  1996) 

2.  Low  chiller  AT  issues 

Kirsner  (1996)  describes  a  problem  afflicting  most  primary-secondary  chilled  water  systems  in  service 
today.  He  describes  a  situation  where  the  chilled  water  temperature  difference  between  the  supply  and 
return  habitually  runs  at  half  the  design  value.  This  creates  three  problems.  First,  at  a  fixed  load,  the 
secondary  loop  requires  twice  the  flow  rate.  This  flow  rate  requires  the  use  of  additional  chillers  to  prevent 
return  to  supply  mixing.  This  mixing  is  a  detrimental  situation  because  it  increases  the  air  handler  water 
supply  temperature,  reducing  heat  transfer.  These  chillers  are  then  running  at  or  below  50%  load.  As  a 
result,  the  system  capacity  is  diminished  to  half  of  its  designed  level.  Furthermore,  the  system  is  operating 
inefficiently,  chillers  are  not  optimally  loaded  and  increased  pump  power  is  required.  Kirsner  states  “  To 
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my  knowledge,  every  chilled  water  plant  serving  distributed  loads  is  afflicted  with  it  to  some  degree.” 
Sauer  (1989)  describes  a  number  of  the  causes. 

Poorly  controlled  or  calibrated  air  handler  controls  are  often  the  culprit.  Air  handler  water  flow 
rate  is  controlled  by  the  cold  deck  temperature.  If  this  temperature  is  adjusted  too  low  the  system  will 
attempt  to  meet  the  set  point  by  increasing  water  flow.  Typically,  the  valves  are  forced  wide  open.  The 
result  is  a  primary-secondary  system  with  constant  flow  in  both  loops. 

When  valves  are  not  properly  calibrated,  hot  and  cold  water  coils  may  be  in  use  simultaneously. 

In  some  air  handlers,  these  coils  share  the  same  deck.  The  reheat  here  directly  increases  the  cold  deck 
temperature.  This  again  leads  to  an  increased  chilled  water  flow. 

Poor  chilled  water  coil  selection,  also,  leads  to  low  temperature  differential.  Since  it  is  difficult  to 
determine  a  particular  zone’s  precise  cooling  requirement  the  load  may  be  underestimated.  The  engineer 
selects  a  coil  that  is  too  small  to  handle  the  load.  Later,  in  the  air  handler  performance  analysis,  it  will  be 
shown  that  an  overloaded  air  handler  delivers  a  reduced  temperature  differential. 

Fouling  is  another  culprit.  Over  time  dirt  builds  on  the  air-side  of  a  chilled  water  coil.  This 
buildup  impedes  heat  transfer  through  the  tube  wall.  The  air  handler  compensates  by  increasing  the  water 
flow  rate  attempting  to  raise  the  UA.  Poor  water  treatment  causes  water-side  fouling  which  has  a  similar 
effect  on  air  handler  performance. 

A  major  design  issue  with  the  primary-secondary  system,  as  Kirsner  (1996)  points  out,  is  that  it  is 
“blinded”  by  low  temperature  differential.  In  some  primary-secondary  systems  chillers  are  scheduled  as  a 
function  of  bypass  flow.  When  secondary  flow  rate  exceeds  primary  flow  rate  another  chiller  is  energized. 
When  primary  flow  rate  exceeds  secondary  flow  rate  by  one  chiller’s  flow  rate  a  chiller  is  turned  off.  This 
strategy  assumes  system  flow  rate  varies  appropriately  with  the  load.  When  a  facility  is  afflicted  with  low 
AT  chillers  are  required  to  stay  in  use  to  handle  the  flow  rather  than  the  cooling  load.  At  times  all  chillers 
may  be  energized  at  a  40%  to  50%  load  in  order  to  prevent  secondary  from  exceeding  primary  flow. 
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Engineers  have  developed  a  number  of  methods  to  deal  with  low  chiller  AT.  These  methods 
include;  over-sizing  air  handlers,  inserting  a  check  valve  in  the  bypass  loop  and  closing  the  bypass 
altogether.  These  design  options  will  be  discussed  in  greater  detail  in  the  body  of  this  text. 

Oblective 

The  purpose  of  this  study  is  to  investigate  bypass  loop  effects  on  the  total  chilled  water  system 
performance.  This  is  accomplished  by  modeling  each  system  component  and  combining  them  into 
comprehensive  system  models.  The  components  include  air  handler,  chiller,  pump  and  cooling  tower 
models.  Combination  of  these  models  allows  the  study  of  system  performance  effects  due  to  changing  a 
single  component. 

Single  loop  and  primary-secondary  systems  are  compared  under  different  situations  to  determine  their 
individual  strengths  and  weaknesses.  Comparisons  were  made  using  the  following  system  pairs: 


System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Single  chiller 

Single  chiller 

Primary  Pumps 

Variable  speed 

Constant  speed 

Secondary  pumps 

N/A 

Variable  speed 

Air  handler  design 
capacity 

100  tons 

100  tons 

Air  handler  quantity 

10 

10 

Air  handler  status 

Unfouled 

Unfouled 

Table  1.  Base  case  comparison 


System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Single  chiller 

Single  chiller 

Primary  Pumps 

Constant  speed 

Constant  speed 

Secondary  pumps 

N/A 

Constant  speed 

Air  handler  design 
capacity 

100  tons 

100  tons 

Air  handler  quantity 

10 

10 

Air  handler  status 

Unfouled 

Unfouled 

Table  2.  Constant  speed  pumping  comparison 


System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Two  chillers 

Two  chillers 

Primary  Pumps 

Variable  speed 

Constant  speed 

Secondary  pumps 

N/A 

Variable  speed 

Air  handler  design 
capacity 

1 00  tons 

100  tons 

Air  handler  quantity 

10 

10 

Air  handler  status 

Unfouled 

Unfouled 

Table  3.  Multiple  chiller  system  comparison 


System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Single  low 
efficiency  chiller 

Single  low 
efficiency  chiller 

Primary  Pumps 

Variable  speed 

Constant  speed 

Secondary  pumps 

N/A 

Variable  speed 

Air  handler  design 
capacity 

1 00  tons 

100  tons 

Air  handier  quantity 

10 

10 

Air  handier  status 

Unfouled 

Unfouled 

Table  4.  Low  efficiency  chiller  comparison 


System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Single  chiller 

Single  chiller 

Primary  Pumps 

Variable  speed 

Constant  speed 

Secondary  pumps 

N/A 

Variable  speed 

Air  handler  design 
capacity 

145  tons 

145  tons 

Air  handler  quantity 

10 

10 

Air  handler  status 

Unfouled 

Unfouled 

Table  5.  Increased  air  handler  capacity  comparison 
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System  Design 

Single  loop 

Primary-secondary 

Plant  type 

Single  chiller 

Single  chiller 

Primary  Pumps 

Variable  speed 

Constant  speed 

Secondary  pumps 

N/A 

Variable  speed 

Air  handler  design 
capacity 

100  tons 

100  tons 

Air  handler  quantity 

10 

10 

Air  handler  status 

Fouled 

Fouled 

Table  6.  Fouled  air  handler  comparison 


Additionally,  a  variety  of  load  profiles  were  tested  on  each  pair  to  ensure  consistent  results.  Each  system 
pair’s  performance  is  compared  by  two  methods.  The  first  compares  the  system  relative  efficiencies  when 
the  actual  load  is  smaller  than  the  design  value.  The  second  compares  the  system  capacities  when  the 
actual  load  exceeds  the  design  value. 

It  is  assumed  that  each  system  is  fully  instrumented  and  calibrated.  This  instrumentation  allows 
the  operator  precise  temperature,  pressure  and  flow  rate  readings  throughout  the  system.  Furthermore,  it  is 
assumed  that  the  operator  executes  the  chosen  control  strategy  perfectly.  In  reality,  this  detailed 
instrumentation  and  precise  control  is  not  possible.  However,  the  primary-secondary  system  is  more  prone 
to  mismanagement.  As  stated  before  the  bypass  loop  blinds  the  plant  to  the  actual  load  conditions.  As  a 
result  optimal  chiller  scheduling  is  more  difficult.  Ideal  operation  is  assumed  in  this  study  in  order  to 
obtain  an  objective  comparison  of  system  performance. 
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CHAPTER  II 

MODELING 

Component  models  were  developed  using  fundamental  principles  and  equations  found  in 
literature.  These  equations  were  compiled  using  the  Engineering  Equation  Solver  software  package  (EES) 
(Klein,  1999).  Engineering  Equation  Solver  was  chosen  for  this  research  due  to  its  ability  to  call  a  variety 
of  thermodynamic  functions.  For  example,  with  the  enthalpy  and  pressure  of  water  known,  all  other 
thermodynamic  quantities  are  readily  available  for  use.  Additionally,  psychrometric  functions  of  air-water 
mixtures  are  available  and  used  in  air  handler  and  cooling  tower  models.  Furthermore,  EES  is  able  to 
iterate  large  numbers  of  transcendental  equations,  an  ability  essential  to  the  system  models. 

The  outputs  from  component  models  have  been  compared  with  actual  manufacturer’s  data  to 
verify  model  accuracy.  The  system  models  were  manipulated  with  different  conditions  and  compared  to 
determine  trends  in  system  performance. 
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Air  Handler 


1.  Model 


The  air  handler  chosen  for  this  study  is  a  variable  air  volume  system  shown  in  Figure  4  (Lorsch 
1993).  This  air  handler  has  only  a  cold  deck.  The  bypass  and  hot  deck  were  not  considered  for  simplicity, 
as  they  would  have  no  effect  on  the  chilled  water  system.  The  air  handler  operates  by  modulating  both  air 
and  water  flow  in  response  to  temperature  change. 


Motor 

Figure  4.  SimpliOed  variable  volume  air  handier 

The  inlet  water  temperature  is  essentially  set  by  the  chiller  evaporator  exit  temperature.  The  pipe 
flow  from  the  chiller  to  the  air  handler  is  assumed  to  be  a  constant  enthalpy  process.  Since  the  air  handler 
water  inlet  pressure  is  smaller  that  the  evaporator  exit  pressure,  the  air  handler  inlet  temperature  is  elevated 
by  a  fraction  of  a  degree.  A  mixture  of  15%  outdoor  to  85%  return  air  is  maintained.  Therefore,  the  inlet 
air  temperature  is  fixed  by  the  room  and  ambient  temperature.  A  constant  75°F  room  temperature  and  55% 
maximum  relative  humidity  are  used  to  ensure  thermal  comfort  (Lorsch  1993).  In  order  to  maintain  this 
condition  the  cold  deck  is  set  to  a  constant  58°F.  The  air  handler  varies  air  flow  by  changing  fan  speed  and 
altering  damper  positions  to  different  zones  in  order  to  maintain  a  constant  room  temperature.  If  a  room 
starts  to  overheat,  more  cold  air  is  directed  into  that  area.  The  chilled  water  flow  rate  adjusts  to  ensure  the 
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cold  deck  remains  58°F.  If  the  cold  deck  temperature  increases,  the  air  handler’s  two-way  valve 
compensates  by  opening  to  allow  more  chilled  water  into  the  coil.  The  two-way  valves,  controlled  by  cold 
deck  temperatures,  dictate  the  entire  system’s  chilled  water  flow  rate. 

The  air  handler  is  modeled  using  NTU-effectiveness  relations.  Air  handler  load  is  assumed  to  be  a 
function  of  outdoor  dry  bulb  temperature.  This  assumption  results  from  Shabo  (1992),  where  a  similar  load 
profile  was  developed  using  data  from  a  commercial  building  in  Atlanta.  The  assumed  load  function  for  an 
Atlanta  office  building  is  shown  below.  Notice  that  at  95^  the  system  is  at  100%  load.  This  corresponds 
to  the  design  conditions  for  Atlanta,  95^  dry  bulb  and  78®F  wet  bulb.  The  load  decreases  linearly  until  it 
reaches  20%  at  60^.  The  load  then  stays  constant  until  the  ambient  temperature  reaches  45^  when  the 
facility  is  assumed  to  switch  to  a  water-side  economizer. 


Atlanta  Load  Profile 


Figure  5.  Air  Handler  Load  Profile 
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The  air  handler  is  modeled  using  three  heat  transfer  equations  (Incropera  and  Dewitt  1985). 


Qoh  =  {T..OU.  - T^.in  )  Equation  1 

Qah='f^'^CP^(Tajn-To.ou,)  Equation  2 

e«/.=^minC^n  ft,,.- Equation  3 

Here  Cpeff  is  the  effective  specific  heat  of  the  air-water  vapor  in  the  coil.  This  specific  heat  accounts  for  the 
latent  as  well  as  sensible  heat.  Its  value  is  given  by: 


CPeff 


^a,m  ^a,oui 

T  .  -T 

a,m  a, out 


Equation  4 


The  enthalpies  are  evaluated  by  use  of  the  psychrometric  charts  in  EES.  The  value  Cmin  is  determined  by 
comparing  the  relative  heat  capacities  of  the  air-side  and  water-side.  The  heat  exchanger  effectiveness  is 
determined  using  two  NTU  relations  for  cross-flow  heat  exchangers  (Incropera  and  Dewitt,  1985).  Here  it 
is  assumed  that  the  water  is  mixed  and  the  air  is  unmixed.  Equation  5  is  used  when  the  water  heat  capacity 
exceeds  that  of  the  air.  Equation  6  is  used  when  the  heat  capacity  of  air  is  largest. 
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Equation  5 


e  =  1  -  exp  j[l  “  exp  (-Q  •  iVrt/ )]  j 


Equation  6 


where: 


Equation  7 


The  number  of  heat  transfer  units  is  given  by: 

IJA 

NTU  =  — -  Equation  8 

c 

min 

Therefore,  the  remaining  unknowns  in  the  initial  three  heat  transfer  equations  are  water  exit  temperature, 
water  and  air  mass  flow  rates  and  the  overall  heat  transfer  conductance  (UA).  In  order  to  solve  this  set  of 
equations,  the  overall  heat  transfer  conductance  is  needed.  The  development  is  described  below. 

a)  Heat  Transfer  Calculations 

i)  Overall  Heat  Transfer  Conductance 

The  chilled  water  coil  UA  is  a  function  of  the  tube  material,  geometry  and  the  individual  heat 
transfer  coefficients.  It  is  given  by  (Oskarsson,  1990): 
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UA  = 


^  1  ^ 


hA 


+ 


In 


D 


'a 


2JZ-k^  -L-n, 


^_i_" 

KAVo 


Equation  9 


In  this  equation  hj  and  ho,  represent  the  inner  and  outer  heat  transfer  coefficients  respectively.  Ao  and  Aj  are 
the  total  inner  and  outer  surface  area.  Dq  and  Di  are  the  tube  inner  and  outer  diameters,  kt  is  the  thermal 
conductivity  of  the  tube,  L  is  the  length  of  one  pass,  and  r|o  is  the  overall  outside  air  handler  effectiveness. 
This  outer  effectiveness  is  related  to  the  fin  efficiency  and  the  air  handler  geometry  using  the  expression. 


where: 


4  = 


Equation  10 


Equation  1 1 


Equation  12 


where.  At  and  Af  are  the  tube  and  fin  surface  areas  and  T^f  is  the  fin  efficiency. 


2)  Air-side  Heat  transfer  coefficient 

Two  methods  were  used  to  model  the  air-side  heat  transfer  coefficient.  These  developments  result 
from  the  work  of  McQuiston  (1981)  and  Rich  (1973).  Each  formulation  is  described  below  as  well  as 
comparisons  with  manufacturer’s  data. 
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3)  McQuiston’s  Method 

McQuiston  finds  the  airside  heat  transfer  coefficient  by  determining  a  factor,  JP,  which  is  a 
function  of  the  air  hander  geometry  as  well  as  flow  characteristics  (McQuiston  1981). 


7P  =  Re; 


-0.4 


1-0.15 


Equation  13 


In  this  equation  Rcq  is  the  Reynolds  number  based  on  the  tube  outer  diameter,  Ao  is  the  outer  surface  area 
(fins  and  tubes),  and  is  the  outer  surface  area  of  tubes  alone. 

The  factor  JP  is  related  to  a  j4  factor  which  determines  a  four-row  coil’s  heat  transfer  performance. 

74  =0.00 14 +  0.26 18  (JP)  Equation  14 


The  variable,  jn  for  a  coil  with  other  than  four  rows  is  found  by, 

l-1280-«-Re/^ 

;4  1 -1280-4  Rej'"" 


Equation  15 


where  n  is  the  number  of  tube  rows  and  Reb  is  the  Reynolds  number  based  on  the  tube  spacing  in  the  flow 
direction.  The  jn  factor  is  then  related  to  the  heat  transfer  coefficient  by, 

Equation  16 

where  Gc  is  the  mass  velocity  through  the  minimum  flow  area,  Cp  is  the  specific  heat  of  the  air  and  Pr  is  the 
Prantl  number  of  the  air. 


4)  Rich’s  Method 

Rich’s  method  differs  in  that  it  assumes  the  coil’s  heat  transfer  coefficient  is  relatively  uniform 
from  the  first  to  the  last  row.  Rich  used  a  four  row  coil  to  obtain  experimental  data.  One  then  expects 
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fairly  consistent  results  when  modeling  a  four  row  coil  Rich  uses  a  j  factor  that  is  a  function  of  the 
Reynolds  number  based  on  the  longitudinal  tube  spacing.  This  j  factor  is  given  by  (Rich  1973): 

7  =  0. 195  •  Equation  17 

Here  the  j  factor  is  related  to  the  outer  heat  transfer  coefficient  in  the  same  manner  as  McQuistion  gives. 

5)  Water-side  Heat  transfer  coefficient 

The  water-side  heat  transfer  coefficient  is  a  function  of  fluid  and  flow  properties.  If  the  flow  in 
the  tubes  is  laminar  then  the  heat  transfer  coefficient  is  given  by  (Bejan,  1995): 

Nu  —  4.36  Equation  18 

Here,  Nu  is  the  Nusselt  number  and  is  related  to  the  inner  heat  transfer  coefficient  by: 


Equation  19 


where,  hj  is  the  inner  heat  transfer  coefficient,  Dj  is  the  tube  inner  diameter,  and  kw  is  the  water’s  thermal 
conductivity.  When  flow  is  turbulent,  the  heat  transfer  coefficient  is  determined  by  the  relation  (Bejan 
1997): 

Nu  “  0.023  •  Equation  20 

where  the  Reynolds  number  here  is  based  on  tube  inner  diameter. 

6)  Manufacturer’s  data 

The  two  methods  presented  above  have  been  compared  with  manufacturer’s  data  to  evaluate  their 
merit.  The  selected  chilled  water  coil  is  shown  in  Figure  6  (Colmac  Coil  Inc.,  1993). 
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Tube  Size 

Outer  Diameter  0.75” 
Inner  Diameter  0.666” 


^ _  Four 

Rows 


Figure  6.  Chilled  Water  Coil 

The  following  figure  shows  how  well  the  calculated  overall  heat  transfer  coefficients  compare  with  the 
manufacturer’s  results  for  the  coil  shown  above.  The  plot  depicts  the  UA  with  air  velocities  from  400  to 
700  ft/min  and  water  velocities  of  two  to  eight  ft/sec. 
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Figure  7.  Accuracy  of  UA  calculation 


It  is  evident  from  these  figures  that  both  formulations  model  the  heat  transfer  well.  In  fact,  the  maximum 
disparity  between  McQuiston’s  method  and  the  manufacturer’s  data  was  1%  and  the  average  difference 
was  0.6%.  Rich’s  model  delivered  a  maximum  difference  of  2%  and  a  average  of  1%.  McQuistion’s 
model  was  chosen  because  it  provides  more  design  flexibility  by  accounting  for  the  number  of  tube  rows, 
b)  Pressure  Drop  Calculation 
1 )  Air-side  Pressure  Drop 

The  airside  pressure  drop  model  used  was  developed  by  Donald  Rich.  It  assumes  that  the  pressure 
drop  across  the  airside  of  the  coil  can  be  broken  into  two  components.  The  first  portion  is  the  pressure  drop 
due  to  the  presence  of  the  tubes  alone,  the  other  part  is  due  to  the  fins  alone.  Rich  relates  a  fin  friction 
factor  to  the  pressure  drop  using  (Rich,  1973): 
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^fin  -  fj 


Equation  21 


where,  ff  is  the  fin  friction  factor  and  p  ^is  the  mean  air  density.  Additionally,  a  is  the  ratio  of  minimum 
air-flow  to  fi'ontal  area  and  Gc  is  the  mass  velocity  through  the  minimum  air-flow  area.  The  fin  friction 
factor  is  related  to  the  Reynolds  number  based  on  the  longitudinal  spacing  by: 


1.7 

Ref 


Equation  22 


The  pressure  drop  due  through  a  staggered  array  of  tubes  is  determined  by  a  relation  developed  by  Kakac 
and  Zukausskas  (Bejan,  1995). 


( 


Equation  23 


Here,  nt  is  the  number  of  rows  in  the  flow  direction,  pa  is  the  air  density,  and  Va^max  is  the  air  velocity  in  the 
minimum  free-flow  area.  The  two  terms  ff  and  %  are  determined  using  the  Figure  8  (Bejan,  1995). 
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Figure  8.  Friction  Factor  for  Arrays  of  Staggered  Tubes 

Finally  the  total  pressure  drop  is  found  by  summing  the  pressure  drop  due  to  the  fins  with  that  due  to  the 
tubes. 

=  AfJ  -f*  Ai^  Equation  24 

2)  Water-side  Pressure  Drop 

The  pressure  drop  across  the  inside  of  the  air  handler  is  assumed  to  be  a  function  of  flow,  tube 
roughness  and  coil  geometry.  First,  the  pressure  drop  attributed  to  the  straight  tube  walls  is  determined. 
Then  the  pressure  drop  due  to  the  tube  bends  and  the  header  are  found.  In  this  calculation,  the  pressure 
drop  due  to  tube  bends  is  a  relatively  large  value  and  thus  cannot  be  neglected.  In  some  situations  the 
pressure  drop  resulting  from  the  tube  bends  and  header  can  equal  or  even  surpass  the  pressure  drop  from 
tubes.  The  friction  factor  for  tubes  is  found  using  (White,  1986): 
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The  pressure  drop  from  the  tube  walls  alone  is  found  using: 


Equation  25 


D, 


The  pressure  drop  is  then  converted  into  a  head  using  a  common  relation: 


Head,  = 


P«,-8 


Equation  26 


Equation  27 


In  this  formula,  the  term  APt  is  the  pressure  drop  from  tubes,  Pw  is  the  density  of  the  water  inside  the  tubes 
and  g  is  the  gravitational  constant.  To  find  the  pressure  drop  due  to  tube  bends  and  header,  it  is  assumed 
that  each  header  contains  two  90°  bends  or  one  180°  bend.  A  single  air  handler  has  two  headers.  So  an  air 
handler  with  four  passes  has  five  180°  bends.  Three  of  these  bends  connect  the  air  handler’s  four  passes 
and  the  last  two  bends  result  from  inlet  and  outlet  headers.  The  head  due  to  these  bends  is  found  using  the 
relation  (White,  1986). 

K 

Head^^^=  Equation  28 


where  Kbcnd  is  a  constant  that  depends  on  tube  diameter  and  bend  geometry. 

3)  Manufacturer’s  data 

Both  air  and  water-side  pressure  drop  formulations  have  been  verified  using  manufacturer’s  data. 
The  coil  used  in  this  analysis  is  shown  in  Figure  6.  Figure  9  shows  the  correlation  between  the  calculated 
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and  actual  water-side  head.  The  maximum  deviation  between  the  calculated  and  actual  results  is  less  than 
8%  with  an  average  difference  of  4%. 


Figure  9.  Accuracy  of  Calculated  Water-side  Head 

The  air-side  pressure  drop  results  are  less  agreeable.  The  model  predicts  the  pressure  drop  fairly  well  at 
high  flow  rates  but  does  not  perform  as  well  at  low  flow  rates.  Figure  10  shows  the  model’s  agreement. 
The  disagreement  between  the  manufacturer’s  data  and  calculated  values  does  not  exceed  5%  when  the  air 
flow  rate  is  between  800  and  400  ft/min.  However,  the  disparity  grows  to  nearly  21%  when  the  flow  rate 


decreases  to  300  ft/min. 
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Pressure  Drop  [inches  H2O] 


Figure  10.  Accuracy  of  Calculated  Air«side  Pressure  Drop 


These  differences  do  not  affect  this  study’s  conclusions.  First,  the  airside  pressure  drop  is  only 
used  to  determine  the  air  handler  fan  power  requirement,  which  accounts  for  less  than  4%  of  the  annual 
energy  consumption.  Secondly,  the  air  handler  fans  are  operated  similarly  in  each  case.  Thus,  this  small 
error  exists  in  every  case  to  an  equal  degree. 

4)  Fan  model 

The  air  handler  fan  is  modeled  similar  to  a  pump.  The  isentropic  efficiency  is  equal  to  the  ratio  of 
the  enthalpy  change  of  an  isentropic  fan  to  that  of  the  actual  fan.  Assuming  a  constant  isentropic  efficiency 
allows  for  calculation  of  the  exit  enthalpy  using. 


K  = 


r  h  _u  ^ 

out, rev  *^a,in 


Vfan 


+  K,in 


Equation  29 


where  ha,out.rev  is  the  exit  enthalpy  of  the  air  if  the  fan  operated  reversibly.  This  reversible  exit  enthalpy  may 
be  determined  because  three  independent  properties  of  the  air-water  vapor  mixture  are  known.  The  exit 
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pressure  is  taken  as  the  ambient  pressure  plus  the  pressure  drop  that  must  be  overcome  through  the  chilled 
water  coil.  The  humidity  ratio  and  the  entropy  are  taken  to  equal  the  fan’s  inlet  condition. 

Once  the  exit  enthalpy  is  known  the  fan  power  requirement  is  determined  using: 


•  • 


Equation  30 


2.  Air  Handler  Performance  Characteristics 

In  order  to  properly  analyze  a  chilled  water  system  it  is  necessary  to  understand  the  dynamics  of 
the  air  handler.  The  air  handler  dictates  the  flow  rate  of  the  entire  system,  and  determines  the  maximum 
chilled  water  AT.  This  air  handler  is  designed  to  have  a  10®F  water-side  temperature  change  under  design 
conditions. 


Figure  11.  Chilled  Water  AT  versus  Load 
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Figure  1 1  shows  how  the  chilled  water  AT  increases  as  air  handler  load  is  reduced.  This  trend  is 
verified  by  research  conduced  by  Sauer  (1989).  Furthermore,  he  discovered  that  a  number  of  factors  may 
affect  the  trend.  Some  of  these  include  improper  air  flow,  improperly  adjusted  set  points  and  fouled  coils. 

Since  the  waterside  temperature  difference  increases  inversely  with  load  the  water-side  flow  rate 
has  more  than  a  linear  decrease. 


Figure  12.  System  flow  rate  versus  load 

This  sharp  decrease  in  flow  rate  has  a  significant  effect  on  the  part  load  pumping  requirements. 

Another  aspect  of  air  handler  performance  is  the  effect  of  fouling.  Air  handlers  get  fouled  on  the 
air-side  and,  to  a  lesser  degree,  on  the  water-side.  The  air-side  fouling  is  a  result  of  small  particles  that 
build  up  on  the  coil’s  outer  surfaces.  Periodic  cleaning  can  minimize  this  problem.  Water-side  fouling 
results  from  water  contaminants  being  deposited  on  tube  walls.  Most  systems  add  chemical  treatment  to 
chilled  water  to  minimize  this  build  up.  Fouling  is  added  to  the  air  handler  model  by  including  a  fouling 
factor  into  the  overall  heat  transfer  equation.  With  the  addition  of  fouling  the  equation  becomes  (Icropera 


and  Dewitt,  1985): 
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UA=f 


f  1  ^ 


^hjAj  j 


V  ) 


R 


f.o 


V  "  ) 


KKno 


Equation  3 1 


Figure  13  shows  the  effect  of  fouling  on  air  handler  performance.  This  air  handler  is  assumed  to 
be  fouled  to  the  point  that  it  has  a  7^  waterside  temperature  difference  at  design  load,  rather  than  the 
intended  10®F. 


Figure  13.  Fouled  air  handler  water-side  T  versus  load 

Notice  that  the  fouled  air  handler  operates  as  if  it  were  a  55  ton  capacity  air  handler  at  design 
conditions,  rather  than  100  tons.  The  fouling  effect  becomes  less  evident  as  load  decreases.  This  means  at 
design  conditions  a  system  with  fouled  air  handlers  requires  a  much  larger  water  flow  rate.  The  increase  in 
flow  rate  is  seen  in  Figure  14. 
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Figure  14.  Fouled  air  handler  flow  rate  versus  load 

This  change  in  flow  rate  has  a  significant  impact  on  pumping  power.  System  analysis  will  show  that  this 
affects  chiller  scheduling  as  well. 

Finally  it  is  important  to  understand  how  the  system  load  relates  to  other  system  variables.  This 
analysis  is  especially  important  when  operating  above  design  conditions.  Consider  the  equation: 

2  =  t/A  •  LMTD  Equation  32 

where  LMTD  is  the  log  mean  temperature  difference  between  the  air  and  water-sides.  The  UA  is 
predominantly  a  function  of  water  and  air  flow  rates.  Consider  four  variables:  the  flow  rates  and  coil  inlet 
temperatures  of  the  water  and  air  streams.  Holding  the  other  three  constant,  if  either  the  air  or  water  flow  is 
increased  the  heat  transfer  rate  rises.  If  the  flow  rates  and  air  inlet  temperature  are  held  constant  reducing 
the  chiller  set  point  will  increase  the  heat  transfer.  Finally,  raising  the  air  inlet  temperature,  in  effect,  the 
outdoor  dry  bulb  temperature  also  increases  the  heat  transfer. 
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Chiller 


1.  Manufacturer’s  data 


The  chiller  components  that  were  modeled  include  the  compressor,  evaporator  and  condenser. 
Manufacturer’s  performance  data  was  used  to  verify  each  model’s  accuracy.  For  proprietary  reasons,  any 
reference  to  manufacturer  or  equipment  has  been  excluded. 


Compressor 


%  Load 

KW/Ton 

100 

EBil 

0.520 

1000 

100 

0.522 

1000 

100 

ESI 

0.523 

66 

EtEI 

0.467 

66 

BB! 

0.470 

66 

nn 

0.474 

330 

33 

0.515 

330 

33 

IQ 

0.524 

330 

33 

IQ 

0.536 

1000 

100 

EHil 

0.445 

1000 

100 

444 

0.444 

100 

0.442 

66 

0.392 

66 

0.397 

1^1 

66 

^3 

0.403 

330 

33 

141 

0.427 

330 

33 

0.439 

330 

33 

IRI 

0.448 

1000 

100 

EQ 

0.398 

100 

E1BI 

0.400 

IWtl 

100 

0.402 

I  660 

66 

I3c!i! 

0.348 

66 

Bn 

0.353 

^^3 

66 

EB3 

0.359 

330 

33 

IBB 

0.364 

330 

33 

IBB 

0.373 

330 

33 

0.382 

Evaporator 


iiffn 

EWT 

Ref.  Temp 

PD(«) 

3300 

51.27 

41.12 

36.37 

El 

54.00 

41.49 

19.73 

El 

IIBB 

65.24 

42.38 

4.63 

El 

IcWiM 

48.80 

41.60 

36.43 

El 

BIBI 

50.60 

41.87 

19.77 

El 

UBil 

58.02 

42.54 

4.62 

El 

46.40 

42.03 

36.49 

El 

BjfiM 

47.30 

42.20 

19.81 

El 

BBl 

51.01 

42.66 

4.62 

El 

lcg!!l 

51.00 

41.08 

39.17 

El 

54.00 

41.49 

19.81 

El 

BBl 

65.24 

42.38 

4.63 

El 

IcliEM 

48.02 

41.47 

51.19 

El 

50.60 

41.87 

19.77 

El 

BK!!I 

58.67 

42.58 

3.89 

El 

46.01 

41.95 

51.06 

El 

KWB 

47.30 

42.20 

19.81 

El 

fBSBI 

50.34 

42.60 

5.67 

El 

50.09 

40.95 

51.12 

El 

KBBl 

54.00 

41.49 

19.73 

El 

BBl 

65.20 

42.27 

4.63 

El 

48.02 

41.47 

51.19 

El 

KMtl 

50.60 

41.87 

19.77 

El 

58.67 

42.58 

3.89 

El 

46.01 

41.95 

51.25 

El 

47.30 

42.20 

19,81 

El 

BBSI 

50.34 

42.60 

5.67 

Condenser 


i??n 

IcBTl 

LWT 

Ref.  Temp 

PD(ft.) 

85 

3000 

94.27 

95.60 

25.71 

85 

94.25 

95.60 

25.71 

85 

3000 

94.23 

95.60 

25.71 

85 

3000 

91.40 

91.87 

25.80 

85 

3000 

91.03 

91.87 

25.80 

85 

3000 

91.02 

91.85 

25.80 

85 

3000 

88.07 

88.46 

25.89 

85 

3000 

88.07 

88.46 

25.89 

85 

88.07 

88.46 

25.89 

84.10 

85.45 

26.32 

75 

EfBBl 

84.06 

85.44 

26.32 

75 

3000 

84.04 

85.41 

26.32 

mm 

IctBfil 

80.92 

81.77 

26,42 

Ire 

3000 

80.91 

81.77 

26.42 

mm 

IcMiTil 

80.91 

81.76 

26.42 

75 

3000 

78.00 

78.40 

26,52 

75 

3000 

78.00 

78.40 

26.52 

iSHil 

78.00 

78.40 

26.52 

65 

3000 

73.97 

75.39 

27.00 

65 

73.96 

75.38 

27.00 

65 

3000 

73.94 

75.37 

27.00 

65 

70.85 

71.74 

27.11 

65 

3000 

70.84 

71.73 

27.11 

65 

3000 

70.84 

71,73 

27.11 

65 

3000 

67.95 

68.36 

27.22 

65 

3000 

67.94 

68.36 

27.22 

65 

67.94 

68.36 

27.22 

Table  7.  Present  Day  Chiller  Data 


The  data  used  in  this  study  is  that  of  a  modem,  1000  ton  capacity,  electrical  chiller.  The  data 


depicts  27  different  operating  conditions,  with  varying  evaporator  flow  rates,  loads  and  condenser  inlet 


temperatures. 
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2.  Compressor  Model 

a)  Weber’s  Compressor  Model 

An  earlier  compressor  model  developed  by  Weber  (1988)  and  used  by  Joyce  (1990)  utilized  the 
ideal,  Carnot  efficiency  multiplied  by  a  constant  isentropic  compressor  efficiency  and  a  part  load  efficiency 
factor.  The  reader  should  consult  Weber  for  the  complete  model  development. 

The  model’s  purpose  is  to  predict  the  compressor  power  required  at  various  loads,  condenser 
pressures  and  evaporator  pressures.  Weber  initially  assumes  an  isentropic  compressor  that  compresses 
refrigerant  from  the  evaporator  pressure  to  the  higher  condenser  pressure.  The  model  then  determines  the 
power  needed  by  this  reversed  Carnot  cycle  using  the  relation: 


camot  — 


evap 


^cond 


evap 


-1 


Equation  33 


The  manufacturer’s  power  requirement  is  then  related  to  that  of  the  reversed  Carnot  cycle  using  an 
isentropic  compressor  efficiency,  ilisent^  design  load.  This  factor  takes  into  account  the  electric  motor 
efficiency  as  well  as  internal  compressor  irreversibility.  Finally,  Weber  includes  an  efficiency  load  factor, 
Titon»  that  decreases  with  increasing  load.  The  final  relationship  is  given  as: 


f  •  ^ 

W  camot 


actual 


Q, 


evap 


•'evap 


(  Visent  Vton  ) 


Equation  34 


Figure  15  shows  present  day  chiller  data  modeled  using  Weber’s  formulation.  The  greatest  error 
encountered  using  this  model  is  9%  while  the  average  error  is  3%. 
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Actual  Compressor  Power  [KW] 


Figure  15.  Weber’s  compressor  model  accuracy 
b)  Liu’s  Compressor  Model 

Liu  (1997)  revised  Weber’s  compressor  model  by  assuming  a  certain  amount  of  the  compressor 
power  is  used  to  overcome  constant  frictional  power.  The  remainder  of  the  power  corresponds  to  the 
Carnot  efficiency. 

•  •  • 

W  actual  =  A  W  carnal  +  W  fiict  Equation  35 

Figure  16  shows  the  present  day  chiller  data  modeled  using  Liu’s  formulation.  The  largest  error  using  this 
approach  is  14%  and  the  average  error  is  5%. 
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Figure  16.  Liu’s  Compressor  Model  Accuracy 
c)  Revised  Compressor  Model 

In  order  to  improve  the  chiller  model’s  ability  to  replicate  present  day  chiller  performance,  a  new 
compressor  model  was  designed.  Like  Liu’s  model,  it  is  assumed  that  a  portion  of  the  power  input  to  the 
compressor  is  needed  to  overcome  frictional  power.  This  frictional  power,  however,  are  assumed  to  be  a 
function  of  chiller  load.  The  remaining  compressor  power  corresponds  to  the  Carnot  power  requirement, 
although  this  is  also  taken  to  be  a  function  of  load. 
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Carnot  Power  [Btu/hr] 


Figure  17.  Compressor  Power  Actual  versus  Carnot 

Figure  17  shows  that  compressor  power  is  not  purely  a  function  of  the  Carnot  efficiency,  but  a  function  of 
load  as  well.  As  a  compressor  proceeds  from  a  large  to  a  smaller  load,  inlet  vanes  in  the  compressor  are 
turned  adding  swirl  to  the  refi-igerant  entering  the  compressor.  This  allows  the  compressor  to  run  at  a 
constant  speed,  while  controlling  the  refrigerant  flow  rate.  However,  adding  this  swirl  increases  frictional 
losses.  The  compressor  motor  efficiency  is  also  a  function  of  load.  A  chiller  is  normally  designed  for  its 
highest  efficiency  at  part  load,  because  it  operates  there  most  of  its  life.  As  a  result,  the  compressor  motor 
appears  to  operates  optimally  between  75%  and  85%  of  full  load.  As  the  load  increases  or  decreases  the 
compressor’s  efficiency  decreases.  The  relation  used  for  finding  compressor  power  is  given  by: 


2 

actual  “  camot  *  fill 


^  • 

+  W  camot  •  fil2 

/ 


J 


Equation  36 


These  individual  functions  are: 


.  2 


Equation  37 
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Where  the  heat  transfer  rate  is  in  units  of  tons  and  the  constants  are  given  in  the  table  below. 


1 

2 

3 

1.31E-12 

2.03E-06 

0.952 

B 

3.15E-09 

4.39E-03 

705 

C 

2.79E-06 

1.23 

3040 

Tables.  Compressor  constants 

Figure  18  shows  the  revised  model’s  accuracy.  The  largest  error  is  3%  and  the  average  is  1.5%. 
The  revised  model  appears  to  replicate  the  present  day  compressor’s  performance  somewhat  better  than 
Weber’s  formulation  and  much  better  than  Liu’s  technique. 


Figure  18*  Revised  Compressor  Model  Accuracy 
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3.  Evaporator  Heat  Exchanger  Model 

a)  Heat  transfer  calculation 

In  a  previous  study  Weber  (1988)  modeled  the  affect  that  varying  condenser  water  flow  rate  has  on 
condenser  heat  exchanger  UA.  Weber  assumed  that  at  a  fixed  load,  the  UA  is  only  a  function  of  flow  rate. 
He  determined  a  relationship  between  UA  and  condenser  flow  rate  at  a  constant  load.  He  assumed  that  at  a 
fixed  flow  rate  the  overall  UA  is  only  a  function  of  load.  Thus,  Weber  simply  shifts  the  UA  versus  GPM 
curve  up  and  down  to  accommodate  different  loads.  This  development  seems  logical,  however,  Weber  had 
only  two  data  points  at  a  loads  below  the  design  load.  As  a  result,  it  was  impossible  for  him  to  validate  his 
assumptions.  This  study  attempts  a  similar  approach  to  model  variable  evaporator  flow. 


Flow  Rate  [gpm] 

Figure  19.  Evaporator  UA  versus  Flow  Rate 

Figure  19  represents  the  relationship  between  evaporator  UA  and  water  flow  rate  for  three 
different  load  conditions.  This  figure  clearly  shows  that  at  lower  flow  rates  the  curves  approach  each  other, 
while  they  separate  at  the  flow  rate  increases. 

The  inner,  water-side,  heat  transfer  coefficient  is  a  function  of  flow  and  fluid  properties,  namely 
the  Reynolds  and  Prantl  numbers.  The  outer,  refrigerant  side,  heat  transfer  coefficient  is  a  function  of  load, 
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or  rate  of  boiling.  In  the  evaporator  refrigerant  is  boiling  on  the  outside  of  tubes.  At  very  small  loads 
natural  convection  starts  around  the  tube  bundle.  As  the  load  increases,  nucleate  boiling  begins.  Initially, 
there  is  only  isolated  bubbles.  As  load  increases  the  boiling  refrigerant  forms  slugs  and  columns  of  vapor. 
While  boiling  transitions  from  isolated  bubbles  to  rapid  boiling  the  refrigerant  becomes  more  agitated, 
increasing  the  outer  heat  transfer  coefficient.  This  explains  why  the  refrigerant  side  heat  transfer 
coefficient  is  tied  to  load. 

Refrigerant  Side 

Tube 
wall 

hAouter=fn(Load) 

Figure  20*  Cross-section  of  a  Single  Tube  in  the  Evaporator  Tube  Bundle 

For  the  given  chiller,  the  inner  tube  diameter  as  well  as  water  velocity  at  a  specified  flow  rate  were 
known.  With  this  knowledge,  the  Reynolds  number  was  found  as  a  function  of  water  flow  rate.  Evaporator 
UA  was  determined  directly  from  the  manufacturer’s  data.  Using  these  values  a  relationship  was 
developed  for  evaporator  UA  as  a  function  of  Reynolds  number  at  a  fixed  load.  This  function  was  then 
evaluated  as  the  Reynolds  number  approached  zero.  It  was  assumed  that  the  UA  below  the  intercept  of  this 
function  and  the  Y-axis  is  due  to  refrigerant  side  conditions.  While  the  portion  of  UA  above  the  intercept  is 


due  to  water-side  conditions. 
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Water  side  Reynolds  number 

F^re  21.  Evaporator  UA  versus  Reynolds  Number 


A  ratio  of  these  components,  ^vap,  was  used  to  reduce  the  number  of  unknowns  in  the  overall  heat  transfer 
coefficient  calculation. 


Equation  38 


then 


= 


1  1 


H 


Equation  39 


It  was  assumed  the  refrigerant  side  hA  remains  constant  with  fixed  load.  Since  the  UA  is  known,  as  well  as 
the  refrigerant  side  hA  at  the  specified  load,  it  was  possible  to  solve  for  the  water-side  hA  as  a  function  of 
Reynolds  number.  This  function  takes  the  form. 
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hA^-K-  Re^  Equation  40 

where  K  and  a  are  constants.  This  relationship  was  refined  knowing  the  heat  transfer  coefficient  is  also  a 
function  of  fluid  properties.  The  heat  transfer  coefficient  for  a  fluid  losing  heat  is  proportional  to  Pr® 

Thus,  the  equation  becomes, 

hAi=D-  Pr^  ^  -  Re"  Equation  41 

where  the  constant  D  depends  on  the  geometry  of  the  tube  bundle.  The  constants  D  and  a  were  determined 
from  experimental  data  and  take  the  values  6210  and  0.533  respectively.  This  value  for  the  constant,  a,  is 
relatively  close  to  the  expected  value,  0.8.  The  difference  is  assumed  to  result  from  heat  transfer 
augmentation. 

With  the  water-side  hA  known  the  refrigerant  side  hA  was  found  as  a  function  of  load,  taking  the 
form  of  a  quadratic. 

f^=‘^-Qevap^  +  B-Qe.„p+C  Equation  42 

Here  the  values  A,  B  and  C  were  found  to  be  1.47, 9,240  and  181,000  respectively. 
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The  method  described  above  models  the  evaporator  heat  transfer  very  well  as  seen  in  Figure  22. 
The  maximum  error  using  this  method  is  2%. 


Figure  22.  Accuracy  Results  for  Evaporator  UA  Calculation 
b)  Pressure  drop  calculation 

The  evaporator  heat  exchanger  water-side  pressure  drop  calculation  is  critical  to  this  study  because 
the  evaporator  makes  a  major  contribution  to  the  chilled  water  system’s  pressure  drop.  Decreasing  water 
flow  rate  through  the  evaporator  at  part  load  reduces  the  pressure  drop  and  with  it  the  pump  power.  Since  a 
primary-secondary  system  keeps  a  constant  flow  rate  through  the  evaporator,  this  savings  is  not  available. 
Bypass  loop  removal,  however,  permits  this  savings. 

Due  to  proprietary  reasons  the  manufacturer  was  unwilling  to  release  any  information  pertaining 
to  tube  bundle  inner  surfaces.  It  is  expected  that  there  is  augmentation  on  the  tubes  inner  surfaces  which 
increases  heat  transfer  between  the  chilled  water  and  tube  wall.  This  augmentation  affects  evaporator 
pressure  drop.  In  order  to  model  the  evaporator  pressure  drop,  an  effective  length  and  roughness  were 
determined.  In  actuality,  the  evaporator  is  much  shorter  and  rougher  than  this  effective  length  and 
roughness  indicate.  Also,  much  of  the  evaporator  pressure  drop  results  from  tube  bends  and  header  losses. 
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However,  since  the  actual  geometry  is  not  available,  effective  values  were  used.  The  following  relations 
were  used  in  this  calculation  (White  1985): 


ftj 


Yi 


=  -1.81og, 


10 


6.9 

_1_ 

(  e  ^ 

1.11“ 

3.7  D, 

^  '  J 

Equation  43 
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Equation  44 


In  order  to  solve  for  Leff  and  Ecff  two  data  points  are  used,  one  was  chosen  at  maximum  evaporator 
flow  and  another  at  minimum  flow.  These  points,  with  the  equations  above,  were  used  to  solve  for  two 
unknowns,  Leff  and  8eff.  It  is  apparent  from  Figure  23  that  the  model  replicates  the  manufacturer’s  data  very 
well  throughout  the  range  of  flow  rates. 


Figure  23.  Accuracy  Results  for  Evaporator  Pressure  Drop 
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4.  Condenser  Heat  Exchanger  Model 

a)  Heat  transfer  calculation 

The  following  equations  were  used  to  determine  the  condenser’s  heat  transfer  performance: 
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Figure  24.  Condenser  UA  versus  Load 

Figure  24  shows  that  the  condenser  UA  decreases  with  an  increased  load.  In  the  presence  of  very  small 
loads,  only  a  small  amount  of  condensate  covers  the  tube  surfaces  in  the  condenser  tube  bundle.  This 
leaves  the  majority  of  the  tube  surface  free  for  direct  contact  with  refrigerant  vapor.  As  load  increases. 
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more  condensate  builds  on  the  tube  walls.  As  a  result,  the  liquid  film  layer  tends  to  insulate  the  refrigerant 
from  the  tube  wall.  When  load  continues  to  increase  this  film  grows  thicker,  hampering  heat  transfer. 
Therefore,  the  condenser  UA  has  an  inverse  relation  with  load. 


Figure  25  shows  that  the  condenser  UA  is  a  function  of  condensing  temperature.  This  is  because  the 
refrigerant  film  thickness  is  affected  by  the  condensing  temperature.  When  the  refrigerant  is  condensing  at 
a  high  temperature,  it  is  less  viscous  and  does  not  cling  to  the  tube  walls  well.  So  it  is  not  able  to  build  a 
thick  film.  When  the  condensing  temperature  decreases,  the  refrigerant  becomes  more  viscous  and  is  able 
to  develop  a  thicker  layer.  This  thick  film  impairs  heat  transfer,  so  the  condenser  UA  decreases  with 
condensing  temperature.  Taking  the  chiller  load  and  condensing  temperature  into  account  the  following 
relation  was  developed. 

=  f{Load  )  ■  f  )  Equation  48 

The  two  functions  were  determined  from  curve  fits  of  manufacturers  data.  The  load  function  is  a  linear 
function,  and  the  temperature  factor  is  assumed  to  be  second  order.  These  functions  are: 
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.  .2 

/(r„^)  =  185000  +  22800.Q.w-96.2-a<«.  Equation  49 

/(load)  =  1 .06  -  9.23  •  10'®  •  Equation  50 


Figure  26  shows  the  model’s  accuracy.  The  maximum  error  was  1.3%  and  the  average  error  was  0.6%. 


Figure  26,  Accuracy  Results  For  Condenser  UA  Calculation 
b)  Pressure  drop  calculation 

The  chilled  water  system  considered  in  this  study  assumes  a  constant  condenser  water  flow  rate  of 
3  gpm  per  ton  capacity.  For  example,  a  single  chiller  plant  with  a  1000  ton  capacity  chiller  has  a  constant 
3000  gpm  condenser  flow  rate.  As  a  result,  the  condenser  pressure  drop  is  almost  constant.  Since  the 
condenser  flow  rate  is  constant ,  the  condenser  water  temperature  is  the  only  variable  that  affects  condenser 
head.  At  cooler  condenser  water  temperatures  the  water  will  become  more  viscous,  causing  a  greater  head. 
These  variations  are  small;  about  a  10%  change  from  the  maximum  to  minimum  head  over  the  range  of 
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temperatures  considered.  Since  the  head  varies  indirectly  with  condenser  temperature,  the  head  is  modeled 
using  the  relation. 

Equation  51 

Where  the  constants  A  and  B  are  60.6  and  0.191respectively.  This  model  results  in  maximum  errors  of  less 
than  0.06%  when  compared  with  manufacturer’s  data. 


Figure  27.  Accuracy  of  condenser  head  model 

5.  Chiller  Comparison 

Data  from  three  chillers  has  been  compared  in  order  to  examine  the  differences  in  chiller  performance. 
The  three  chillers  studied  are  all  1000  ton  electric  chillers.  One  is  the  current  high  efficiency  chiller  shown 
in  Table  7.  For  the  purpose  of  this  discussion  this  chiller  is  labeled  chiller  A.  The  other  two  include  a  high 
efficiency  chiller  and  a  lower  efficiency  chiller,  both  available  and  popular  in  1995.  These  chillers  are 
labeled  as  chillers  B  and  C  respectively.  Data  from  these  older  chillers  is  presented  below. 
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75 

44 

3000 

2400 

304 

83.14 

41.46 

600 

75 

44 

3000 

2400 

261 

91.96 

41.61 

500 

75 

44 

3000 

2400 

222 

80.79 

41.76 

400 

75 

44 

3000 

2400 

185 

79.65 

41.91 

Table  9.  Liu’s  High  Efficiency  Chiller  (B) 


Table  10.  Liu’s  Low  Efficiency  Chiller  (C) 
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Performance  comparisons  were  made  on  the  compressors,  evaporators  and  condensers.  The 
figures  below  show  the  relative  performance  of  these  three  chillers.  Figure  28  shows  the  relative 
performance  of  the  three  chiller  compressors.  All  three  compressors  have  similar  performance  at  full  load. 
In  fact,  the  chiller  B  appears  to  be  more  efficient  than  chiller  A  full  load.  As  load  decreases,  the  chillers  A 
and  B  continue  to  perform  similarly.  At  33%  load  the  chiller  A’s  compressor  is  somewhat  more  efficient 
than  that  of  chiller  B  but  by  only  a  small  amount.  Chiller  C’s  compressor  performance  drops  off 
dramatically  as  load  decreases.  This  performance  trend  may  be  attributed  to  increased  design  focus  on 
achieving  a  high  integrated  part  load  value  (IPLV),  rather  than  full  load  performance.  The  IPLV  places 
great  importance  on  part  load  efficiency. 


Figure  28.  Compressor  Performance  Comparison 

Figure  29  shows  that  the  evaporator  heat  exchangers  have  a  similar  UAs  at  33%  load.  As  loads 
increase,  the  chillers  B  and  C  continue  to  have  similar  UA’s.  The  chiller  A’s  UA  increases  dramatically 
with  load.  At  full  load  the  current  chiller’s  evaporator  performs  much  better  than  both  older  chillers. 
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Figure  29.  Evaporator  Performance  Comparison 

Figure  30  shows  that  chiller  A’s  condenser  UA  is  much  higher  than  the  other  two  chillers  at  all 
loads  depicted,  while  the  chillers  A  and  B  have  similar  U A’ s. 
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Figure  30.  Condenser  Performance  Comparison 
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These  figures  indicate  that  chiller  A’s  heat  exchangers  are  significantly  improved  over  the  chillers 
used  by  Liu.  The  differences  may  include  different  tube  surface  treatments  or  larger  heat  exchangers.  Both 
variations  not  only  increase  heat  transfer  but  also  could  be  expected  to  increase  the  chiller’s  first  cost. 


6.  Performance  Characteristics 

The  following  discussion  refers  to  the  performance  characteristics  of  chiller  A.  In  order  to 
appreciate  chilled  water  system  performance  it  is  necessary  to  grasp  how  different  conditions  affect  chiller 
operation.  The  figure  below  shows  the  effect  that  changing  load  has  chiller  efficiency. 


Figure  31.  Chiller  efficiency  versus  load 

With  high  condenser  water  temperatures  the  chiller  is  most  efficient  at  about  75%  of  maximum 
load.  As  the  condensing  water  temperature  is  decreased  the  chiller  becomes  more  efficient  because 
compressor  lift  is  reduced.  Additionally,  the  chiller’s  best  efficiency  point  shifts  to  lower  loads  as 
condenser  temperature  decreases.  At  a  65®F  condenser  inlet  temperature  the  best  efficiency  is  down  to 


about  65%. 
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Increasing  flow  rate  through  the  evaporator  at  fixed  load  has  two  different  affects.  First,  the 
evaporator  UA  increases  which  tends  to  increase  heat  transfer  rate.  Meanwhile,  at  a  fixed  chilled  water  exit 
set  point,  the  LMTD  between  the  water  and  refrigerant  is  decreased.  This  tends  to  decrease  the  heat 
transfer  rate.  George  Redden  (1996)  studied  this  phenomenon  to  determine  variable  flow  effects  on  chiller 
performance.  This  was  accomplished  by  modeling  and  bench  testing  a  225  ton  capacity  centrifugal  chiller. 
In  one  case  the  evaporator  flow  rate  was  fixed  while  the  load  was  varied.  In  another  case  the  flow  rate  was 
varied  proportional  to  the  desired  load  reduction.  Comparisons  were  made  data  at  71%,  45%  and  25%  load 
conditions.  Similar  chiller  efficiencies  were  attained  by  either  varying  flow  rate  or  inlet  temperature.  In 
fact  the  variable  and  constant  flow  results  agreed  to  within  1%  of  one  another.  Redden’s  work  concluded 
that  varying  evaporator  water  flow  rate  has  little  effect  on  chiller  performance. 
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Figure  32.  Effect  of  variable  flow  on  chiller  performance 

Figure  32  shows  that  varying  flow  rate  at  a  constant  load  has  very  little  affect  on  this  model’s 
performance,  as  well.  If  the  flow  is  kept  within  acceptable  bounds,  the  chiller  efficiency  is  a  function  only 
of  load  and  lift.  For  this  study’s  purpose  maximum  and  minimum  evaporator  water  velocities  are  taken  to 
be  three  and  1 1  ft/sec  respectively  (Eppelheimer,  1996).  These  values  agree  with  manufacturer’s  data.  The 
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minimum  flow  rate  is  set  to  prevent  laminar  flow  in  the  evaporator.  Transition  to  laminar  flow  would 
significantly  hinder  water-side  heat  transfer  coefficient.  A  velocity  of  three  ft/sec  corresponds  to  a 
Reynolds  number  of  about  15,000,  which  is  turbulent  assuming  the  tube  bundle  is  composed  of  rough 
pipes.  The  maximum  flow  limit  is  imposed  to  prevent  tube  bundle  erosion.  Operation  with  velocities  in 
excess  of  1 1  ft/sec  also  significantly  increases  evaporator  head.  The  increased  pressure  drop  can  erode  the 
tubes  inner  surfaces,  shortening  chiller  life. 

Another  important  characteristic  of  chiller  performance  is  the  effect  that  set  point  has  on 
efficiency.  If  one  increases  the  set  point  with  a  fixed  load,  the  compressor  lift  is  reduced.  This,  in  turn, 
decreases  the  power  requirement.  Figure  33  how  a  chillers  efficiency  increases  as  the  chiller  set  point  is 
raised.  This  figure  depicts  a  the  chiller’s  efficiency  rating  in  kW/ton  while  operating  at  1000  tons. 


Figure  33,  Efficiency  versus  set  point 
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Pump  Model 

Pumps  for  all  models  were  selected  using  curves  found  in  pump  selection  software  (Goulds, 

1995).  An  effort  was  made  to  ensure  all  pumps  had  similar  best  efficiency  points  (BEPs).  For  example, 
the  best  pump  available  for  use  as  the  secondary  pump  in  a  primary/secondary  system,  had  a  83%  BEP.  As 
a  result,  a  variable  speed  primary  pump  with  a  83%  BEP  was  selected  for  the  single  loop  system  even 
though  more  efficient  options  were  available.  Pumps  were  selected  by  first  determining  the  flow  rate  and 
head  required  under  the  conditions  of  greatest  annual  system  energy  consumption.  Pumps  with  an  83% 
BEP  at  this  head  and  flow  rate  were  chosen.  Next,  the  maximum  system  flow  rate  and  head  were 
determined.  All  selected  pumps  were  required  to  meet  these  maximum  values.  For  example,  the 
maximum  flow  required  by  the  system  may  be  2800  gpm,  however  it  operates  most  often  at  2200  gpm. 
Thus,  the  selected  pump  would  have  a  best  efficiency  point  of  83%  at  2200  gpm  and  may  have  and 
efficiency  as  low  at  76%  at  maximum  conditions. 

1 .  Constant  Speed  Pumping 

When  constant  speed  pumps  are  used,  the  air  handlers  dictate  system  water  flow  rate.  Pump  head 
is  a  function  of  the  flow  rate,  and  the  air  handler’s  two-way  valves  modulate  to  accommodate  the 
fluctuating  head.  Figure  34  shows  how  the  system  head  is  forced  to  increase  as  flow  decreases. 
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Figure  34.  Centrifugal  Pump  and  System  Curves 

Three  separate  pump  equations  are  used  to  model  constant  speed  pumps.  These  equations  determine,  pump 
head,  efficiency  and  shaft  power,  all  as  a  function  of  flow  rate. 


^head  ^head  ■GPM  +  C^GRM^ 


Equation  52 


Vpump  =  A,  +  ■  GPM  ^  Equation  53 

Equation  54 


Here  the  constants  in  each  equation  are  chosen  to  describe  the  particular  pump.  These  values  are  given  in 
the  tables  below. 


Head 

eta 

■SOI 

A 

71.4 

0.036 

71600 

B 

0.0154 

7.33E-04 

16.9 

C 

-8.46E-06 

-1.70E-07 

N/A 

Table  11.  Single  loop,  constant  speed,  primary  pump 
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Head 

eta 

bhp 

A 

103 

0.0859 

94O0O 

B 

0.0105 

5.97E-04 

31.8 

C 

-6.31  E-06 

-1.20E-07 

N/A 

Table  12.  Primary-secondary  system,  constant  speed,  secondary  pump 


Here  the  pump  head  and  pump  shaft  power  have  units  of  feet  and  Btu/hr  respectively.  In  order  to  find  the 
total  power  required,  the  pump  shaft  power  is  divided  by  the  constant  motor  efficiency,  taken  as  90%.  This 
gives  the  relation  for  the  total  pumping  power  requirement. 


W 


pump 


shaft 

^motor 


Equation  55 


Since  the  pressure  at  the  pump  exit  is  determined  by  equation  52  only  one  other  quantity  is  needed  to  fix 
the  pump  exit  state.  The  enthalpy  is  chosen  and  is  determined  by: 


h  +  .  Equation  56 

*^w,out  •  'Vm  ^ 

2.  Variable  Speed  Pumping 

When  variable  speed  pumps  are  used,  the  air  handler  dictates  flow  rate.  System  head  changes  with 
flow  rate. .  The  pump  speed  is  controlled  by  a  pressure  sensor  in  a  remote  air  handler.  The  pump  speed  is 
adjusted  to  keep  this  differential  pressure  between  the  chilled  water  coil  inlet  and  exit  constant.  Even 
though  this  distant  chilled  water  coil’s  pressure  drop  is  constant,  the  pump  head  varies.  This  is  because 
there  are  other  components  that  contribute  to  the  total  system  head.  These  include  piping,  tube  bends  and 
evaporator  head  losses.  As  a  result,  the  pump  operates  at  the  speed  which  allows  it  to  maintain  the  flow 
rate  and  head  determined  by  the  system.  Figure  35  describes  a  variable  speed  pump  under  normal 
operating  conditions. 
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A  pump  is  selected  to  operate  most  efficiently  at  some  specific  flow  rate  and  head.  These 
conditions  are  defined  as  the  pump’s  best  efficiency  point.  This  point  is  labeled  as  1  in  the  figure.  If  the  air 
handler  requires  more  flow  the  pump  speed  increases,  pumping  more  water.  The  pump  head  increases  and 
the  pump  operates  at  a  new  point,  On  the  other  hand,  if  the  air  handler  requires  less  flow,  the  pump  will 
slow,  diminishing  the  pump  head.  In  this  case  the  new  operating  point  is  denoted  by  1". 

Design 

Decreasing  System 


Figure  35.  Effects  of  Changing  Speed  on  the  Pump  Curve 

In  order  to  determine  the  pump  operating  characteristics  two  equations  are  needed.  The  first  determines  the 
head  as  a  function  of  flow  rate  when  the  pump  is  operating  at  its  design  speed.  This  design  speed  is  chosen 
for  the  condition  that  the  pump  will  operate  most  frequently.  This  equation  describes  a  curve  called  the 
design  pump  curve  (DPC).  The  second  equation  describes  pump  efficiency  as  a  function  of  flow  rate  when 
operating  at  design  speed.  This  equation  describes  the  curve  called  the  design  efficiency  curve  shown 


below. 
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Flow  Rate  [gpm] 


Figure  36.  Design  efficiency  curve 

However,  the  pump  may  not  be  operating  at  the  design  speed.  Therefore  it  is  first  necessary  to  determine 
the  actual  speed  of  the  pump.  This  is  done  by  use  of  the  affinity  laws. 


N 


actual 


Equation  57 


In  this  equation,  the  term  GPMdpc  and  Ndcsign  are  the  flow  rate  and  speed  of  the  pump  when  operating  on 
the  design  pump  curve.  The  corresponding  design  pump  curve  head  is  also  found  using  an  affinity  relation. 


Hecidpp^  — 


^^C^^actual 


Equation  58 


The  design  pump  curve  head  and  efficiency  are  known  as  a  function  of  the  design  curve  flow  rate  given  by 
two  equations. 
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Vdpc  =  a,  +  •  GPM opc  +  C,  •  GPMopc  Equation  59 

Head^PC  =  \ead  +  ^head  '  GPM  dpc  +  C^ad  '  ^PM  Equation  60 

The  constants  in  these  equations  are  chosen  to  describe  a  particular  pump.  These  values  are  displayed  in 
the  tables  below. 


Head 

eta 

A 

103 

0.0859 

B 

0.0105 

5.97E-04 

C 

-6.31  E-06 

-1.20E-07 

Table  13.  Single  loop,  variable  speed,  primary  pump 


Head 

eta 

A 

71.4 

0.036 

B 

0.0154 

7.33E-04 

C 

-8.46E-06 

-1.70E-07 

Table  14.  Primary-secondary  system,  variable  speed,  secondary  pump 


The  four  preceding  equations  can  be  used  to  define  the  pump’s  operating  conditions.  However,  actual 
pump  efficiency  decreases  as  pump  speed  is  reduced.  This  reduction  in  efficiency  is  small  when  the  pump 
is  operating  within  45%  of  its  design  speed.  At  slower  speeds  the  efficiency  starts  to  diminish  more  rapidly 
(Rishel,  1996).  The  actual  efficiency  is  determined  using  a  relation  taking  the  form 


pump 


f 

“  Vdpc  ’  / 

V 


\ 

/ 


where: 


Equation  61 


=  0.815 +  0.403  • 


actual  j 


N 

design 


-0.218 


N. 


actual  / 


'N. 


design 


Equation  62 


The  exit  conditions  are  then  determined  by  equation  60. 
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^w,out 


^w,out,rev  ^w,in 


pump 


wan 


The  total  pumping  power  is  calculated  using: 


W  pump  — 


^ motor  fvsd 


Equation  63 


Equation  64 


where  the  term  Timotor/vsd  is  assumed  to  be  a  function  of  speed.  It  is  a  combination  of  motor  and  variable 
speed  drive  efficiency.  This  hardware  is  selected  to  operate  most  efficiently  at  the  BEP  and  will  perform 
less  efficiently  at  greater  or  lesser  speeds.  The  expression  for  this  combined  efficiency  is  taken  fi*om  Joyce 
(1990). 


ri^,or,ysd  =0.0498  +  1.85 


N 

actual 

N 

design 


\ 


-1.03 


N 

^  ’  actual 

N 

design 


Equation  65 


The  motor  and  variable  speed  drive  efficiency  is  shown  in  Figure  37. 
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Cooling  Tower  Model 

It  was  necessary  to  include  a  cooling  tower  as  a  component  of  the  system,  because  cooling  tower 
operation  dictates  condenser  water  inlet  temperature,  which  affects  compressor  lift.  As  a  result,  the 
ambient  weather  conditions  affect  system  performance  through  the  cooling  tower  as  well  as  air  handlers. 

The  cooling  tower  used  in  this  study  was  adapted  from  the  model  developed  by  Weber(1988). 
However,  it  includes  the  counterflow  NTU-effectiveness  model  developed  by  Braun(1988)  then  used  by 
Liu  (1997)  and  Lilly(1998).  The  reader  should  consult  these  works  for  the  complete  development  of  the 
mass,  energy  and  diffusion  equations  used  to  model  the  cooling  tower. 

The  following  assumptions  were  made  in  modeling  the  cooling  tower: 

1)  Heat  and  mass  transfer  only  in  the  direction  normal  to  water  and  air  flow. 

2)  Heat  and  mass  transfer  through  the  tower  walls  to  the  environment  are  negligible. 

3)  Constant  air  and  water  specific  heats  for  a  given  operational  situation. 

4)  Heat  transfer  from  the  tower  fans  to  the  air  is  negligible. 

5)  Uniform  temperature  throughout  the  water  stream  in  any  cross  section  normal  to  the  flow. 

6)  Uniform  tower  cross  sectional  area. 

Using  steady  state  energy  and  mass  balances  on  a  differential  control  volume,  Braun  developed 
three  differential  equations.  The  first  of  these  is  derived  from  an  energy  balance  on  the  water  stream. 


^  ^  ^  Bqmtion  66 


The  next  equation  is  taken  from  an  energy  balance  on  the  air  stream. 
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Equation  67 


The  third  equation  is  derived  from  an  overall  mass  balance. 
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Equation  68 


where: 


Equation  69 
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Equation  70 


Given  the  inlet  conditions  equations  63, 64  and  65  can  be  solved  for  the  water  stream  exit  condition.  This 
solution  requires  integration  over  the  entire  tower  volume  from  inlet  to  exit.  The  process  is  simplified 
using  Merkel’s  assumptions.  Merkel  neglects  the  effect  the  water  loss  due  to  evaporation  on  water  flow 
rate  and  sets  the  Lewis  number  equal  to  unity.  Equations  63  and  64  reduce  to  equation  68. 
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Equation  72 


Braun  then  introduced  a  saturation  specific  heat,  Cs.  The  saturation  specific  heat  is  the  derivative  of 
saturation  air  enthalpy  with  respect  to  the  air  temperature  at  the  air/water  interface.  At  this  interface  the 
temperature  is  assumed  equal  to  the  water  temperature.  So  Cs  is  defined  as. 
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Equation  73 


Equation  71  may  then  be  expressed  in  terms  of  air  enthalpies. 
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Equation  74 


The  exit  conditions  can  are  for  analytically  using  equations  69  and  71,  assuming  that  the  saturated  air- water 
vapor  enthalpy  is  a  linear  function  of  temperature.  In  other  words  the  saturation  specific  heat,  Cs,  is 
assumed  to  be  uniform  throughout  the  tower. 


Figure  38.  Defining  a  Linear  Saturation  Enthalpy 


Figure  38  depicts  how  the  saturation  enthalpy  varies  with  the  saturation  temperature.  It  is  obvious  that  the 
saturation  enthalpy  does  not  vary  linearly  with  temperature.  Although,  by  choosing  an  appropriate  average 
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slope  between  the  inlet  and  outlet  conditions,  as  seen  in  the  figure,  an  effective  relationship  is  found  in 
terms  of  Cj. 
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Equation  75 


The  air-side  effectiveness  is  defined  as  the  ratio  of  the  actual  heat  transfer  to  the  maximum  heat 
transfer  possible  if  the  exiting  air  stream  were  saturated  at  the  incoming  water  temperature.  The  actual  heat 
transfer  in  terms  of  the  effectiveness  is. 
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Equation  76 


where  the  air-side  effectiveness  is  found  by  using  the  common  counter-flow  heat  exchanger  relation. 
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Equation  77 


Here  C,  is  the  heat  capacity  ratio,  given  by  equation  13. 


Equation  78 


The  outlet  water  enthalpy  is  then  be  determined  by  performing  an  overall  energy  balance. 


Equation  79 


1 .  Cooling  Tower  NTU  Calculation 

The  cooling  tower  NTU  is  given  by  the  equation  developed  by  Lowe  and  Christie  (1960). 
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Equation  80 


In  this  expression  c  and  n  constants,  specific  to  a  particular  tower  cooling  box  geometry.  Manufacturers 
offer  towers  that  have  many  different  box  sizes,  each  having  its  own  cooling  capacity.  To  calculate  the 
coefficients  and  verify  model  accuracy,  Liu  (1997)  correlated  manufacturer’s  data  for  a  number  of  box 
sizes.  Manufacturer’s  data  gives  tower  outlet  water  temperature  for  specified  water  and  air  flow  rates.  The 
reader  should  consult  Liu’s  work  for  a  complete  development.  The  cooling  tower  box  size  used  in  this 
study  is  given  below  (Lilly,  1998): 


Length 

Width 

Height 

13’ 11  1/8" 

24’  1/2” 

18’ 7  5/8" 

Table  15.  Tower  box  geometry 


2.  Cooling  Tower  Fan  Model 

Many  different  fan  motors  are  available  for  each  box  size.  A  manufacturer  typically  gives  the  fan 
shaft  power  for  each  air  flow  rate  specified  with  a  particular  box  size.  Liu  shows  the  relationship  between 
fan  horsepower  and  air  flow  rate.  He  also  provides  curve  fits  for  each  box  size  used  in  the  tower  model  to 
translate  the  flow  rate  into  fan  power.  The  equation  for  fan  horsepower  is  (Liu  1997): 
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Equation  81 


where  A  and  B  are  constants  which  correspond  to  a  specific  box  size  and  cfm  is  the  air  flow  rate  in  cubic 
feet  per  minute.  In  this  study  values  for  A  and  B  were  taken  to  be  2.31  and  6.38  respectively  (Lilly,  1998). 

In  order  to  maintain  a  minimum  condenser  water  inlet  set  point  fans  are  assumed  to  be  turned  on 
and  off  as  appropriate.  A  70®  F  tower  minimum  set  point  was  used  in  this  study,  which  means  that  fans  run 
continuously  when  the  tower  exit  water  temperature  exceeds  70®  F.  When  the  water  temperature  is  equal  to 
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70°  F,  the  fans  are  cycled  on  and  off  to  maintain  that  temperature.  The  ratio  of  time  that  the  fans  are  in  use 
is  called  the  fan  duty.  Additionally,  it  was  assumed  that  heat  transfer  from  the  tower  with  fans  off  usage  is 
10%  of  the  heat  transfer  attained  with  all  fans  on.  The  fan  power  was  then  found  by  multiplying  the  fan 
horsepower  by  the  fan  duty. 
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Pioina  Model 

The  component  models  detailed  above  were  combined  to  form  system  models.  In  order  to 
combine  models  all  pipe  flow  was  considered  at  constant  enthalpy.  Pipe  head  was  determined  using  the 
relation: 
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the  friction  factor  was  found  using  the  equation  (White,  1986): 
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Equation  84 


Where  Rei>i  is  the  Reynolds  number  based  on  the  tube  inner  diameter,  Et  is  the  tube  roughness  and 
Dj  is  the  inner  diameter.  The  roughness  was  assumed  to  be  0.{XX)15  ft.  A  tube  diameter  was  chosen  such 
that  the  water  velocity  at  design  conditions  approximately  10  ft/sec.  As  a  result,  all  single  chiller  systems 
have  pipes  with  10”  inner  diameters  while  multiple  chiller  system  piping  has  13”  inner  diameters.  This  10 
ft/sec  value  is  a  common  rule  for  office  building  design,  where  little  growth  is  anticipated.  Industrial 
applications  normally  use  a  lower  value,  typically  as  low  as  7.5  ft/sec.  This  allows  for  growth  without 
additional  piping. 

A  10  ft/sec  water  velocity  will  result  in  head  losses  close  to  3  feet  per  100  feet  of  pipe.  The 
primary  loop  of  each  primary-secondary  system  has  500  equivalent  pipe  feet.  This  length  accounts  for  the 
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elbows  and  fittings  as  well  as  the  actual  pipe.  Thus,  at  design  conditions,  the  primary  loop  head  is 
approximately  35  feet.  This  accounts  for  15  of  pipe  head  and  a  20  foot  evaporator  head.  The  secondary 
loop  is  1500  equivalent  feet.  A  differential  pressure  sensor  across  a  distant  air  handler  maintains  a  constant 
20  foot  head  across  its  chilled  water  coil  and  two-way  valve.  As  a  result,  the  secondary  loop  head  is  65  feet 
at  design  load.  Finally,  the  single  loop  system  is  assumed  to  have  2000  equivalent  feet  of  piping. 

Therefore,  the  chilled  water  loop  head  is  100  feet  at  design  conditions.  This  value  represents  20  feet  from 
the  evaporator,  20  feet  from  the  air  handler  and  the  remaining  60  feet  comes  from  pipe  losses.  All  air 
handlers  are  assumed  to  be  at  the  same  elevation  as  the  chiller.  The  condenser  loop  is  assumed  to  be  500 
feet  long.  The  cooling  tower  is  elevated  200  feet  above  the  chiller.  The  cooling  tower  is  assumed  to  be 
elevated  because  these  towers  are  normally  placed  on  roofs  of  large  office  buildings.  These  lengths  were 
imposed  to  add  realism  to  the  system  models.  However,  the  actual  lengths  make  little  difference  in  the 
comparative  results  found  in  the  study. 
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CHAPTER  Hi 

OPERATIONAL  STRATEGIES 

An  effort  has  been  made  to  ensure  the  system  designs  compared  in  this  study  use  control  strategies 

common  throughout  the  HVAC  industry.  Additionally,  similar  strategies  were  employed  in  all  designs 
» 

which  were  compared  with  one  another.  The  only  difference  between  design  pairs  was  bypass  inclusion  or 
omission.  This  was  done  in  order  to  highlight  only  how  the  bypass  loop  affects  system  performance.  The 
control  strategies  employed  were  not  optimized  because  ideal  operation  will  differ  in  each  case.  Instead 
each  model  was  controlled  consistent  with  current  industry  methods.  In  order  to  analyze  system  results,  an 
understanding  of  the  different  methods  is  vital. 

Chilled  Water  Reset 

A  conunon  chilled  water  plant  practice  is  to  raise  the  centrifugal  chiller’s  set  point  as  the  facility’s 
cooling  load  drops.  This  is  referred  to  as  chilled  water  reset.  It  is  believed  to  be  an  inexpensive  energy 
savings  method  (Erth,  1987).  By  increasing  the  chilled  water  supply  temperature  the  compressor  lift  is 
attenuated.  Since  compressor  power  is  a  function  of  load  and  lift,  decreasing  the  lift  while  maintaining  a 
constant  cooling  load  reduces  the  power  requirement.  When  the  chiller  set  point  is  raised,  the  air  handler’s 
inlet  water  temperature  is  increased  as  well.  This  decreases  the  air  handler’s  air  to  water  LMTD.  In  order 
to  maintain  a  similar  heat  transfer  rate  the  coil’s  UA  must  increase.  This  may  only  be  done  by  increasing 
the  water  and  air  velocities.  This  translates  to  an  increased  fan  and  pump  power.  Recent  studies  have 
shown  the  chiller  to  be  the  chilled  water  system’s  largest  energy  consumer.  The  chiller  is  responsible  for 
approximately  62%  of  the  energy,  while  fans  and  pumps  account  for  only  38%  (Eppelheimer,  1996). 

Reset  was  employed  in  this  study  by  increasing  the  chiller  set  point  and  as  the  outdoor  temperature 
decreases  in  order  to  fix  the  system  flow  rate.  This  strategy  is  initiated  only  when  the  outdoor  temperature 
falls  below  a  prescribed  level.  Initially,  air  handlers  control  water  flow  rate  by  modulating  two-way  valves 
in  order  to  maintain  a  58^  cold  deck  temperature.  Once  the  prescribed  outdoor  temperature  is  reached,  the 
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system  flow  rate  is  fixed  and  chiller  set  point  is  raised  to  maintain  the  deck  temperature.  In  reality  this 
ideal  control  is  not  possible.  Outdoor  temperature  affects  various  building  zones  differently.  As  a  result, 
air  handler  cold  deck  temperatures  fluctuate  differently  depending  on  location.  In  practice,  the  chiller  set 
point  is  increased  incrementally  as  outdoor  temperature  drops.  For  example,  the  set  point  may  be  raised 
one  degree  when  the  outdoor  temperature  reaches  55®F  and  then  another  degree  at  50®F.  Figure  39  shows 
this  strategy  as  employed  in  this  research.  Here,  Atlanta  weather  data  was  used  where  design  conditions 
are  95°F  dry  bulb  and  78°F  wet  bulb.  The  chiller  set  point  was  increased  starting  at  an  outdoor  temperature 
of  74®F  and  water  velocity  in  the  evaporator  was  fixed  at  the  minimum  three  ft/sec. 


Figure  39.  Chilled  water  reset 

One  can  see  the  chilled  water  supply  temperature  is  held  at  a  constant  44®F  until  the  outdoor  temperature 
reaches  74®F.  The  chilled  water  set  point  then  increases  inversely  with  cooling  load. 
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Figure  40.  Flovf  rate  variation  with  reset 

Figure  40  represents  the  system  flow  rate  response  to  this  procedure.  At  large  loads  the  system  flow 
decreases  with  load.  However,  after  reset  is  employed  the  system  flow  rate  is  fixed. 
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Figure  41  depicts  the  manner  in  which  chilled  water  AT  is  affected  by  reset.  At  design  load  the  temperature 
difference  between  the  air  handler  water  inlet  and  exit  is  10^.  As  the  load  is  reduced  the  temperature 
difference  increases.  This  trend  continues  until  the  reset  condition  is  reached,  at  which  point  the  water 
flow  rate  is  fixed,  but  the  cooling  load  requirement  continues  to  decline  with  outdoor  temperature.  As  a 
result  temperature  difference  across  the  water-side  of  the  air  handler  starts  to  decrease.  The  temperature 
difference  after  reset  decreases  because  the  amount  of  heat  being  carried  away  by  this  constant  stream  is 
reduced. 

Chilled  water  reset  was  used  in  all  systems  analyzed  in  this  study.  The  technique  was  not 
optimized  because  the  results  would  only  apply  to  a  particular  design  and  be  of  no  great  worth  to  the 
industry.  Instead,  reset  was  employed  to  ensure  the  evaporator  water  velocity  in  single  loop  systems  stayed 
within  manufacturer’s  limits.  The  same  reset  conditions  were  applied  to  all  single  loop/primary-secondary 
system  design  comparisons.  In  each  case  the  application  of  reset  improved  annual  system  efficiency. 


Figure  42  describes  the  total  power  requirement  for  a  primary-secondary  system.  The  lower  curve  is  the 
power  required  for  the  system  using  reset  while  the  upper  curve  is  the  power  required  by  the  system 
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without  reset.  At  small  loads  chilled  water  reset  can  result  in  a  19%  energy  savings.  Chilled  water  reset  is 
required  in  some  single  loop  systems  to  meet  the  minimum  evaporator  flow  requirements,  however  this 
should  not  be  viewed  as  a  deficiency.  Reset  should  be  used  in  any  chilled  water  system  in  the  interest  of 
efficiency.  As  a  result,  the  use  of  this  strategy  is  commonplace  in  industry  today.  . 

Tower  minimum  set  point 

Most  facilities  attempt  to  save  tower  fan  power  by  imposing  a  minimum  tower  set  point.  When 
the  tower  exit  water  exceeds  this  temperature,  all  tower  fans  run  continuously.  As  the  water  drops  to  the  set 
point,  the  tower  fans  are  cycled  to  maintain  the  set  point.  For  example,  if  a  tower  has  five  fans  it  will  run 
only  the  fans  required  to  maintain  the  set  point.  In  order  to  extend  the  tower  fans’  lives,  a  lead/lag  fan 
rotation  is  imposed.  If  water  temperature  continues  to  fall  eventually  all  fans  are  turned  off.  Tower 
minimum  set  point  has  two  advantages.  The  first  is  a  trade  off,  as  is  the  chilled  water  reset  procedure.  The 
tower  set  point  dictates  the  condenser  inlet  temperature.  Thus,  when  the  tower  set  point  is  increased  the 
compressor  lift  is  raised  as  well.  This  reduces  the  chiller  efficiency,  but  this  is  offset  by  diminished  tower 
fan  power  requirements.  Tower  minimum  set  point  also  lengthens  chiller  life.  When  condenser  inlet 
temperature  is  too  low,  the  compressor  lift  becomes  small.  Compressors  are  not  designed  to  operate  with 
small  lifts  and  this  use  may  damage  them.  Liu  (1997)  found  the  optimum  tower  set  point  to  be  between 
65®F  and  68^.  Joyce  (1990)  recommends  a  set  point  of  80°F.  The  optimal  set  point  is  largely  a  function  of 
the  system  design  and  component  performance. 

For  the  purposes  of  this  study,  70®F  was  used  as  a  set  point.  A  fan  duty  is  determined.  Fan  duty  is 
the  percentage  of  time  the  fans  are  operational.  For  example,  if  three  of  the  five  fans  are  running  the  fan 
duty  is  0.6.  The  fan  power  requirement  is  then  found  by  multiplying  the  fan  duty  by  the  power  required 
when  all  fans  are  on.  Additionally,  it  was  assumed  that  the  tower  heat  transfer  rate  without  fans  is  10%  the 
rate  with  all  fans  running. 
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Free  Cooling 

Many  chiller  plants  shut  down  all  chillers  when  the  outdoor  temperature  falls  below  some 
designated  level.  At  this  point,  the  heat  is  transported  to  the  outdoor  air  by  a  more  direct  path,  free  cooling. 
Free  cooling  is  divided  into  two  methods,  air-side  and  water-side. 

Air-side  free  cooling  meets  the  required  cooling  load  by  circulating  outside  air  directly  into  a 
building.  When  this  method  is  used  all  chillers  and  pumps  sit  idle.  The  only  power  requirement  is  fan 
power  needed  to  blow  the  air  into  ventilation  ducts.  This  technique  is  applicable  when  temperatures  fall 
below  about  55®F.  These  factors  make  this  a  very  efficient  method.  Nonetheless,  its  use  is  uncommon  in 
large  facilities.  Normally,  at  low  outdoor  temperatures  only  internal  zones  require  cooling.  Circulating  air 
to  these  areas  requires  a  ventilation  shaft.  The  shaft  is  usually  placed  vertically  down  the  center  of  the 
building  and  takes  outdoor  air  from  the  roof.  Building  owners  seldom  approve  this  method  because  it  is 
more  cost  effective  for  a  building  owner  to  lease  the  space  that  would  be  occupied  by  the  shaft.  These 
profits,  typically,  more  than  offset  the  additional  energy  cost  required  by  other  designs. 

Water-side  free  cooling  is  a  method  by  which  outdoor  air  cools  system  water  without  chiller  use. 

It  is  normally  used  when  ambient  temperatures  fall  below  about  45^.  This  type  of  free  cooling  may  be 
completed  directly  or  indirectly.  Direct,  water-side,  free  cooling  is  accomplished  by  taking  chilled  water 
and  passing  it  directly  from  air  handler  to  tower.  The  problem  with  this  technique  is  that  it  causes  the 
chilled  water  loop  to  be  open  to  the  atmosphere.  Since  direct,  water-side,  free  cooling  vents  the  chilled 
water  vapor  directly  to  the  atmosphere,  oxygen  is  dissolved  by  the  water.  This  oxygen  then  causes 
corrosion  in  chilled  water  system  components.  Additionally,  the  use  of  additives  in  this  design  is 
prohibited  because  they  may  be  vented  to  the  atmosphere.  These  stipulations  make  direct  water-side  free 
cooling  an  unattractive  option. 

Indirect,  water-side,  free  cooling  keeps  the  evaporator  and  condenser  streams  separate.  It  does  so 
by  using  an  additional  flat  plate  heat  exchanger,  an  economizer.  Its  use  requires  additional  pumping  power 
on  both  the  condenser  and  evaporator  sides.  The  additional  power  results  from  economizer  pressure  drop. 
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However,  this  option  maintains  a  closed  chilled  water  system  which  reduces  corrosion  and  waterside 
fouling. 

For  this  study,  it  was  assumed  the  free  cooling  energy  required  for  a  primary-secondary  and  single 
loop  system  is  identical.  As  a  result,  energy  expenditures  are  not  compared  for  ambient  temperatures 
below  45^. 

Multiple  Chiller  Scheduling 

Optimum  chiller  loading  in  multiple  chiller  plants  is  a  point  of  discussion  throughout  the  HVAC 
field.  Most  multiple  chiller  plants  operate  based  on  the  widely  accepted  rule  that  optimum  chiller  loading 
is  between  60%  and  80%  of  full  load.  As  a  result,  a  lead  chiller  is  rarely  allowed  to  exceed  80%  capacity. 
Austin  (1991)  conducted  a  study  that  disputes  this  strategy.  He  found  that  optimal  scheduling  may  involve 
operating  the  lead  chiller  at  loads  as  large  as  100%,  depending  on  system  design.  Much  of  the  confusion 
regarding  scheduling  results  from  Air-Conditioning  &  Refrigeration  Institute  (ARJ)  standards,  which  are 
applicable  only  to  single  chiller  facilities.  Chiller  efficiency  curves  are  normally  rated  at  ARI  conditions. 
These  curves  are  not  adequate  tools  by  which  to  formulate  a  scheduling  procedure.  ARI  standards  dictate  a 
85®F  condenser  inlet  temperature  at  100%  load.  The  condenser  inlet  temperature  is  then  reduced  to  75®F  at 
75%  load  and  further  to  65®  at  50%  and  25%  load.  Figure  43  shows  is  the  efficiency  curve  for  the  chiller 
used  in  this  study  rated  with  ARI  conditions. 
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Figure  43.  ARI  efficiency  curve 

Inspection  of  this  curve  would  lead  one  to  believe  that  optimal  scheduling  would  requires  energizing  the 
lag  chiller  when  the  lead  is  only  65%  loaded.  Use  of  these  standards  is  only  applicable  for  single  chiller 
system  analysis.  In  multiple  chiller  plants  condenser  water  inlet  temperatures  do  not  vary  with  load  as  they 
do  in  single  chiller  facilities.  When  focused  on  scheduling,  it  is  more  informative  to  know  relative 
efficiencies  at  a  constant  condenser  inlet  temperature.  Figure  44  is  an  efficiency  curve  for  the  same  chiller 
at  a  constant  70®F  condenser  inlet  temperature. 
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Evaporator  load  [tons] 

Figure  44.  Chiller  efficiency  curve  with  constant  condenser  inlet  temperature 

This  curve  leads  one  to  believe  the  optimal  loading  requires  energizing  a  lag  chiller  when  the  lead  reaches 
88%  load.  This  conclusion  is  more  reasonable  but  a  more  appropriate  method  still  exists. 

A  lag  chiller  should  be  energized  only  when  the  total  system  energy  requirement  is  decreased  by 
doing  so  (Austin,  1991).  In  the  case  highlighted  above,  it  is  more  efficient  to  operate  two  chillers  at  45% 
load  than  one  at  90%.  However,  one  must  consider  the  additional  auxiliary  power  requirements  associated 
with  a  lag  chiller.  Figure  45  shows  a  comparison  of  the  total  system  power  requirement  for  a  system  using 
one  and  two  chillers.  This  figure  shows  that  optimal  chiller  scheduling  requires  starting  up  the  lag  chiller 
when  the  outdoor  temperature  is  about  73°F. 
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Outdoor  temperature 


Figure  45.  Chiller  scheduling 

This  temperature  corresponds  to  a  load  of  970  tons.  Thus,  in  this  situation  optimal  scheduling  requires 
energizing  the  lag  chiller  when  the  lead  is  at  about  97%  of  its  rated  load.  Although  this  criteria  differs 
depending  on  system  design,  and  condenser  water  temperature. 


Chilled  Water  Set  Point  Reduction 

Central  plant  operators  often  reduce  chiller  set  point  in  an  attempt  to  increase  chilled  water  system 
capacity.  This  strategy  not  only  fails  to  increase  capacity,  but  also  hurts  efficiency.  When  the  chiller  set 
point  is  reduced  the  compressor  lift  is  increased.  This  raises  compressor  power  requirement.  If  a  chiller  is 
operating  at  maximum  capacity  lowering  the  set  point  decreases  its  capacity.  Figure  46  represents  the 
maximum  capacity  of  the  chiller  as  calculated  by  the  system  model  using  588  kW  as  the  maximum  chiller 


motor  power. 


Chiller  capacity  [tons] 


90 


Chiller  set  point  [F] 


Figure  46.  Effect  of  reducing  chiller  set  point 

Notice  that  this  1000  ton  nominaliy  rated  chiiier  is  able  to  deliver  more  than  1000  tons  depending  on  the 
ambient  conditions.  This  figure  also  verifies  the  trend  that  lowering  the  chiller  set  point  reduces  chiller 


capacity. 
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CHAPTER  IV 
DESIGN  OPTIONS 

Engineers  have  introduced  some  design  changes  to  the  typical  primary-secondary  system  in  an  to 
minimize  the  problems  associated  with  low  chiller  AT.  These  changes  are  discussed  below. 

Increasing  Air  Handler  Capacity 

Typically,  a  facility  is  designed  to  deliver  a  10®F  chilled  water  temperature  difference  across  the 
air  handler  at  some  specified  design  condition.  In  effect,  the  chilled  water  coil  is  sized  to  meet  the  design 
load  while  increasing  the  water  temperature  by  10®F,  Increasing  this  temperature  difference  requires  a 
larger  chilled  water  coil.  Engineers  may  design  facilities  to  have  excess  capacity.  This  is  done  for  a 
number  of  reasons.  First,  it  is  impossible  to  accurately  calculate  all  the  factors  that  contribute  to  the 
building  load.  Also,  additional  zones  may  be  added  to  a  particular  air  handler’s  load  after  construction  is 
complete.  Additionally,  pumping  power  is  reduced  with  the  reduced  flow  rate.  Over-sizing  air  handlers 
will  compensate  for  some  fouling,  minimizing  a  possible  low  temperature  differential  problem.  Finally, 
increasing  the  design  AT  may  raise  a  system’s  capacity.  By  decreasing  secondary  loop  flow  rate,  bypass 
mixing  in  the  chilled  water  supply  is  minimized.  Furthermore,  increasing  chilled  water  coil  size  raises  the 
UA. 

Increasing  design  AT  has  a  down  side  as  well.  Doing  so  means  purchasing  and  installing  much 
larger  air  handlers.  This  raises  the  equipment’s  first  cost  as  well  as  the  space  requirement. 

Check  Valve 

Kirsner  (1998)  recommends  another  solution  to  prevent  bypass  mixing  from  chilled  water  return 
to  supply.  He  proposes  placing  a  check  valve  in  the  bypass  preventing  this  flow.  Including  this  valve 
allows  primary  loop  flow  to  exceed  secondary  loop  flow,  but  not  the  converse.  The  primary-secondary 
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system  shown  in  Figure  47  operates  normally  when  the  secondary  loop  flow  rate  is  less  than  or  equal  to  the 
primary  loop  rate. 


Variable  speed 


Constant  speed 
primary  pump 


Figure  47.  Primary-secondary  system  with  check  valve 

As  the  load  requirement  increases  from  the  design  condition,  the  secondary  flow  rate  increases.  However, 
the  check  valve  prevents  secondary  flow  from  exceeding  primary  flow.  As  a  result,  both  primary  and 
secondary  flow  rates  increase  together.  The  problem  with  this  arrangement,  in  relation  to  a  single  loop 
design,  is  that  it  maintains  a  constant  evaporator  water  pressure  drop  below  design  conditions.  Additional 
savings  are  available  by  reducing  the  evaporator  pressure  drop  with  load. 

Single  Loop.  Variable  Flow  Rate 

The  primary  (single  loop)  system  referred  to  here  should  not  be  confused  with  the  old,  constant 
flow  primary  system  with  three-way  valves  in  the  chilled  water  coils.  The  single  loop  system  has  variable 
speed  primary  pumps  and  two-way  valves  in  the  chilled  water  coils.  There  are  no  secondary  pumps  and  the 
bypass  loop  is  omitted.  Single  loop  systems  have  a  number  of  advantages  over  the  primary-secondary 
system.  First  there  is  no  possibility  of  detrimental  mixing  from  return  to  supply.  Thus  for  overloaded 
situations,  this  design  operates  much  like  the  primary-secondary  arrangement  with  check  valve.  Additional 
advantages  are  realized  in  part  load  operation.  Variable  evaporator  flow  reduces  water-side  pressure  drop 
at  small  cooling  loads,  decreasing  pumping  power  requirements.  This  system  eliminates  the  need  for 
secondary  pumps,  diminishing  system  first  cost.  Finally,  since  single  loop  chiller  scheduling  is  determined 
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by  chiller  load  as  opposed  to  flow,  this  design  is  able  to  adapt  to  low  temperature  differential  without  losing 
system  capacity. 
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CHAPTER  V 

LOAD  PROFILES  AND  WEATHER  DATA 

The  analysis  in  this  study  has  been  divided  into  two  categories.  All  systems  compared  are  able  to 
deliver  the  cooling  load  at  design  conditions.  However,  some  designs  deliver  this  cooling  load  more 
efficiently  than  others.  Efficiency  is  critical  during  part  load  operation  because  a  chilled  water  system 
operates  at  part  load  the  vast  majority  of  its  life.  System  capacity  is  only  a  consideration  when  the  required 
cooling  load  is  greater  than  the  design  value.  Some  system  designs  have  greater  capacity  than  others.  The 
performance  categories  compared  in  this  study  are: 

1)  System  capacity  is  analyzed  for  all  situations  where  the  actual  cooling  load  exceeds  the 
system’s  design  load. 

2)  System  efficiency  is  studied  for  all  loads  equal  to  or  less  than  the  system  design  load. 
Efficiency  is  quantified  by  annual  energy  requirements  and  an  integrated  energy  rating  (lER).  In  order  to 
assign  values  to  these  quantities  a  load  profile  was  developed.  This  profile  expresses  system  load  as  a 
function  of  outdoor  dry  bulb  temperature.  Ambient  conditions  and  number  of  annual  hours  spent  at  a 
specified  load  is  defined  by  annual  weather  data. 

The  load  profile  used  here  is  a  modified  version  of  that  used  by  Weber  (1988),  Liu  (1997)  and 
Lilly  (1998).  Their  profiles  assumed  100%  design  load  at  95^.  The  load  drops  linearly  with  temperature 
to  20%  at  60^.  The  load  requirement  is  then  Exed  at  20%  for  all  temperatures  below  60^.  Additionally, 
this  profile  stays  a  constant  100%  at  all  temperatures  above  95®F.  This  represents  the  maximum  system 
capacity  when  operating  above  design  conditions.  The  building  load  is  assumed  to  continue  its  linear 
increase  above  95®. 
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Figure  48.  Previous  load  profile 

A  number  of  modifications  have  been  made  to  this  profile  in  the  course  of  this  research.  First,  the 
supplied  load  is  not  assumed  to  be  equal  to  the  design  load  when  operating  above  design  conditions.  The 
supplied  load,  above  design  conditions,  is  set  equal  to  the  calculated  system  capacity.  This  capacity  is  a 
function  of  ambient  conditions  and  system  design.  For  example,  a  particular  single  loop  system  may  be 
able  to  meet  the  actual  load  up  to  lOO^F.  While  a  similar  primary-secondary  system  may  only  meet  the 
load  up  to  95^.  Once  a  system  is  operating  at  full  capacity  its  maximum  load  may  increase  or  decrease 
with  variations  in  ambient  conditions.  Another  change  made  to  the  profile  is  the  addition  of  free  cooling  at 
45®F.  As  a  result,  it  is  assumed  that  there  is  no  load  below  45^.  In  actuality  there  is  some  cooling  load 
present,  although  the  energy  required  by  free  cooling  is  not  compared.  Finally,  the  design  dry  bulb 
temperature  is  adjusted  for  different  locations.  For  example,  the  design  dry  bulb  in  Atlanta  is  95^  while  in 
Phoenix  it  is  107^.  The  load  profiles  and  weather  data  are  presented  below.  Only  weather  data  at  or 
below  design  conditions  is  represented.  This  is  because  annual  efficiency  comparisons  only  deal  with 
temperature  bins  at  or  below  design  conditions. 
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50  60  70  80 

Outdoor  dry  bulb  [F] 

Figure  51.  Chicago  Load  Profile 


Dry  Bulb  [F] 

Wet  Bulb  [F] 

Total 

Observations 

[hours] 

93 

74 

N/A 

92 

74 

58 

87 

72 

165 

82 

70 

324 

77 

67 

487 

72 

64 

681 

67 

61 

759 

62 

57 

700 

57 

52 

604 

52 

47 

581 

47 

43 

565 

45 

41 

N/A 

Design 

93®  [F] 

77”  [F] 

Condition 

dry  bulb 

wet  bulb 

Table  18.  Chicago  weather  data 
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CHAPTER  VI 

PERFORMANCE  STUDIES 

Chilled  water  system  analysis  has  been  divided  into  two  areas.  These  areas  are  operation  below 
and  above  design  load.  They  are  shown  in  Figure  53.  The  analysis  was  separated  into  two  regions  because 
some  designs  are  able  to  match  load  requirements  at  high  temperatures  while  others  cannot.  All  systems 
are  able  to  meet  the  requirement  below  the  design  temperature. 


40®  50®  60®  70®  80®  90®  100®  110® 

Outdoor  temperature 

Figure  53.  LoadproHle 

This  is  because  a  system  is  designed  to  have  the  necessary  capacity  to  meet  the  design  load.  As  a 
result,  each  design  is  able  to  meet  the  required  load  when  conditions  are  equal  to  or  below  design 
conditions.  However,  some  designs  perform  more  efficiently  than  others  under  these  conditions.  Two 
parameters  are  calculated  to  compare  system  efficiency.  One  parameter  calculated  is  each  design’s  annual 
energy  requirement.  These  values  are  compared  to  determine  the  energy  saved  or  lost  as  a  result  of 
omitting  the  bypass  loop.  Integrated  energy  ratings  (EER)  are  also  compared.  These  values  include  the 
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system,  chiller  and  evaporator  auxiliary  lERs.  The  system  lER  is  the  annual  kW/ton  rating  which 
compares  the  system  power  required  to  handle  the  building  cooling  load.  Chiller  lER  is  a  similar  rating 
which  compares  compressor  power  to  evaporator  load.  Evaporator  auxiliary  BER  rates  the  efficiency  of  all 
air  handler  fans  and  evaporator  side  pumps.  lER  values  are  determined  using  the  formula. 


lER  = - ; - —  Equation  85 

hrj 

Y 

/torif 


where,  hfi,  kWi  and  toili  are  the  number  of  hours,  kilowatts  and  tons  in  the  i*  temperature  bin 


respectively.  The  variable  n  is  the  total  number  of  temperature  bins  considered. 


System  capacity  is  used  to  compare  each  design’s  ability  to  meet  load  requirements  above  design 
weather  conditions.  This  over-design  analysis  is  done  by  comparing  a  primary-secondary  system’s 


capacity  to  that  of  a  similar  single  loop  design. 


Annual  efficiency 


1 .  Base  Case 

This  case  compares  the  efficiency  of  a  single  loop  with  a  primary-secondary  system.  Atlanta 
weather  data  and  load  profile  are  used.  Design  conditions  are  at  95^  dry  bulb  and  78®F  wet  bulb.  All 
handlers  are  sized  for  a  100  ton  load  and  a  10®F  water-side  temperature  differential  at  design  conditions. 
Each  model  has  ten  air  handlers,  giving  the  system  a  1000  ton  capacity  at  design  conditions.  The  two 
systems  are  depicted  in  Figure  54. 
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pump  Single  loop  system 


Constant  speed 
pump 


Variable  speed 
pump 
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Secondary  loop 


◄ - 

Primary-secondary  system 


Figure  54.  Base  case  designs 

The  single  loop  system  gives  a  7.8%  annua!  energy  savings  over  the  primary-secondary  system 
under  these  conditions.  The  vast  majority  of  this  savings  is  derived  from  the  reduced  pump  head  at  part 
load.  The  decreased  head  is  a  result  of  lower  part  load  evaporator  flow  rates.  Figure  55  shows  the  single 
loop  pump  power  reduction. 


103 


Outdoor  Temperature  [F] 


Figure  55.  Base  case  pump  power 

The  compressor  power  required  by  the  primary-secondary  system  is  slightly  larger  as  well.  More 
pump  power  means  that  more  shaft  work  is  added  to  the  water  stream.  This  energy  has  to  leave  through  the 
evaporator  as  heat.  In  fact,  the  primary-secondary  compressor  power  is  larger  but  by  a  very  small  amount. 
The  energy  added  by  pumps  as  shaft  work  is  much  smaller  than  that  added  by  air  handlers  as  heat.  The 
relative  compressor  power  comparison  is  seen  in  Figure  56. 
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The  lER  analysis  is  seen  in  Table  20.  This  table  shows  that  the  savings  is  primarily  a  result  of 
reduced  evaporator  side  pump  power.  The  single  loop  chiller  EER  is  somewhat  larger  that  of  the  primary¬ 
secondary  system,  however  the  difference  is  less  than  one  percent.  The  values  here  are  only  carried  to  two 
significant  digits  because  the  model  is  not  meant  to  predict  absolute  values  within  one  percent  accurately. 
It  should  be  viewed  as  a  comparative  tool  designed  to  point  out  trends. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.042  [kWAon] 

0.085  [kW/ton] 

Chiller  lER 

0.37[kWAon] 

0.37  [kW/ton] 

System  lER 

0.55[kW/ton] 

0.61  [kW/ton] 

Table  20.  Base  case  lER  analysis 


The  chiller  lER  value  may  seem  smaller  than  expected.  This  indicates  very  high  chiller  efficiency. 
The  chiller  efficiency,  here,  is  increased  because  chiller  reset  was  used. 
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2.  Constant  Speed  Pumping  Design 

The  design  and  conditions  for  this  comparison  are  almost  the  same  as  the  base  case.  The  only 
change  is  the  replacement  of  all  variable  with  constant  speed  pumps.  This  situation  was  included  because 
of  the  higher  cost  of  variable  speed  drives.  Certainly  when  attempting  to  capitalize  on  savings  available  by 
varying  evaporator  flow  rate,  energy  costs  will  be  less  for  a  variable  speed  pumping  scheme.  However, 
constant  pumps  require  less  shaft  power  when  flow  rate  is  reduced.  One  would  then  expect  the  single  loop 
system  to  also  offer  a  small  savings  over  the  primary-secondary  system. 

The  single  loop  system  discussed  here  reduces  the  annual  energy  use  by  3.2%.  This  is  a  much 
smaller  savings  than  was  available  with  variable  speed  pumping.  Table  21  shows  that  the  single  loop 
efficiency  increase  is  again  derived  from  a  reduced  pump  power.  It  is  also  interesting  to  note  both  auxiliary 
efficiency  ratings  are  much  worse  than  their  variable  speed  counterparts. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.11[kW/ton] 

0.13  [kW/ton] 

Chiller  lER 

0.37  [kW/ton] 

0.37  [kW/ton] 

System  lER 

0.65  [kWAon] 

0.67  [kW/ton] 

Table  21.  Constant  speed  pumping  lER  analysis 


3.  Other  Load  Profiles 

Weather  data  and  load  profiles  were  selected  from  different  cities  within  the  United  States.  This 
was  done  to  verify  the  savings  available  in  different  HVAC  markets.  The  tables  below  show  the 
comparison  of  single  loop  to  primary-secondary  loop  performance  using  these  different  profiles.  All  single 
loop  systems  are  more  efficient  that  their  primary-secondary  loop  counterpart.  Bypass  loop  omission 
results  in  5.2%,  7.8%  and  8.8%  annual  energy  savings  in  Phoenix,  Atlanta  and  Chicago  respectively.  The 
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single  loop  system  saves  much  more  energy  in  Chicago  and  Atlanta  than  it  does  in  Phoenix.  This  is 
because  the  single  loop  saves  energy  from  the  smaller  evaporator  pressure  drop  at  diminished  loads.  When 
the  average  load  is  lower  the  system  saves  more  energy.  The  average  load  in  Phoenix  is  45%  of  design, 
while  the  average  in  Atlanta  and  Chicago  is  40.8%  and  38%  respectively. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.048  [kW/ton] 

0.084  [kW/ton] 

Chiller  lER 

0.38  [kW/ton] 

0.38  [kW/ton] 

System  lER 

0.54  [kWAon] 

0.59  [kW/ton] 

Table  22.  Phoenix  lER  analysis 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.045  [kW/ton] 

0.092  [kWAon] 

Chiller  lER 

0.37  [kWAon] 

0.37  [kWAon] 

System  lER 

0.55  [kWAon] 

0.62  [kWAon] 

Table  23.  Chicago  lER  analysis 


a)  ARI 

Comparisons  using  the  ARI  load  profile  were  included  because  it  is  an  industry  standard  for  rating 
chillers.  As  such,  these  values  may  easily  be  compared  to  real  world  systems  with  comparable 
components.  The  chiller  lER  here  is  somewhat  lower  than  the  actual  chiller  rating.  This  is  because  chiller 
reset  was  initiated  at  approximately  75®F. 
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System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.038  [kW/ton] 

0.067  [kW/ton] 

Chiller  lER 

0.37[kWAon] 

0.37[kW/ton] 

System  lER 

0.64  [kW/ton] 

0.68  [kW/ton] 

Table  24.  ARI  Load  analysis 


System  IPLV  is  much  higher  in  this  comparison  than  for  other  cooling  load  profiles.  This  is 
because  the  average  system  load  using  this  data  is  58%  of  design.  Since  this  system  weights  the  higher 
loads  more  heavily,  one  would  expect  a  smaller  savings.  The  single  loop  system  saves  5.5%  compared  to 
the  primary-secondary  system  when  using  this  profile. 

4.  Effects  of  Air  Handler  Degradation 

This  case,  again,  compares  the  annual  performance  of  single  loop  and  primary-secondary  systems. 
All  air  handlers  are  fouled  to  the  point  that  there  is  a  7^  design  temperature  differential  between  the  chilled 
water  supply  and  return.  Atlanta  weather  data  and  load  profile  are  used.  Since  there  is  a  7^  AT  at  design 
conditions,  the  system  requires  a  much  higher  flow  rate.  At  full  load  this  flow  is  greater  than  the  constant 
2400  gpm  primary  loop  flow  in  the  primary-system.  As  a  result,  the  secondary  flow  rate  exceeds  primary 
loop  flow  rate.  The  return  water  is  mixed  with  the  chilled  water  supply,  increasing  the  supply  temperature. 
This  situation  will  be  explained  further  in  the  over-design  load  analysis.  In  this  case,  the  secondary  flow 
rate  exceeds  the  primary  flow  when  the  building  load  reached  85%  of  its  design  value.  The  system  is  still 
able  to  meet  the  load  requirement  but  does  so  much  less  efficiently  it  requires  because  more  fan  and  pump 
power.  The  single  loop  design  is  able  to  increase  water  flow  rate  to  air  handlers  without  mixing. 
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System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.045  [kWAon] 

0.089  [kW/ton] 

Chiller  iER 

0.38  [kW/ton] 

0.38  [kW/ton] 

System  IER 

0.56  [kW/ton] 

0.61[kW/ton] 

Table  25.  Fouled  system  lER  analysis 


In  this  case  the  single  loop  system  requires  9.4%  less  energy  annually  than  the  primary-secondary 
system.  In  this  situation  a  fouled  air  handler  was  chosen  for  simplicity.  The  low  temperature  differential 
caused  here  by  fouling  may  result  from  a  number  of  other  factors.  This  analysis  shows  that  the  single  loop 
system  is  more  capable  of  adapting  to  load  problems  that  lead  to  low  temperature  differential. 

5.  Increased  Air  Handler  Capacity 

This  case  determines  the  effect  that  increasing  air  handler  size  has  on  system  efficiency.  This 
situation  is  the  same  as  the  base  case  except  it  employs  larger  air  handlers.  This  system  is  designed  to 
produce  a  12®F  AT  at  design  conditions.  As  a  result,  the  air  handlers  are  much  larger.  The  system’s  10  air 
handlers  are  capable  of  transferring  145  tons  each  with  a  10^  water-side  temperature  change. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
IER 

0.037  [kWAon] 

0.050  [kWAon] 

Chiller  IER 

0.37  [kWAon] 

0.37  [kWAon] 

System  IER 

0.54  [kWAon] 

0.60  [kWAon] 

Table  26.  Increased  air  handler  capacity  lER  analysis 
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In  this  case  use  of  the  single  loop  system  results  in  an  annual  8.2%  energy  reduction.  This  is 
somewhat  larger  than  the  7.8%  seen  in  the  base  case.  Most  of  the  extra  savings  is  a  result  of  reduced 
pumping  power  requirement.  It  is  interesting  to  note  that  increasing  air  handler  size  raised  the  single  loop 
and  primary-secondary  system  efficiencies  by  1.2%  and  1.0%  respectively. 

6.  Low  Efficiency  Chiller  Analysis 

The  chiller  used  in  the  analysis  to  this  point  has  been  a  modem,  high  efficiency  chiller.  Its  IPLV  rating 
is  0.38.  Few  existing  chilled  water  plants  have  high  efficiency  chillers.  Most  plants  operate  with 
equipment  that  was  purchasedTO  to  12  years  ago.  These  chillers  operate  inefficiently  by  today’s  standards. 
Since  the  single  loop  system  derives  it  savings  from  a  decreased  pump  power  at  low  loads,  similar  savings 
should  be  available  with  lower  efficiency  chillers.  Table  27  shows  the  relative  performance  of  the  two 
systems  with  an  older  chiller. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.042  [kW/ton] 

0.085  [kW/ton] 

Chiller  lER 

0.71  [kW/ton] 

0.71  [kW/ton] 

System  lER 

0.89  [kW/ton] 

0.95[kWAon] 

Table  27.  Low  efficiency  chiller  analysis 


This  system  was  modeled  using  chiller  C  from  the  proceeding  chiller  comparison  study.  It  has  a  0.59  IPLV 
rating.  Curiously,  the  chiller  lER  is  much  lower  than  the  IPLV  rating  indicates.  This  is  because  this  chiller 
performs  very  poorly  at  small  loads,  as  was  seen  in  the  chiller  comparison  study.  When  using  Atlanta 
weather  data  and  load  profile  the  average  load  and  the  average  condenser  water  inlet  temperature  are 
approximately  400  tons  and  75®F  respectively.  Under  these  conditions  the  chiller  efficiency  rating  is  0.72 
kW/ton,  much  worse  than  the  IPLV  value.  The  single  loop  system  uses  5.9%  less  energy  annually  than  the 
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primary-secondary  system  in  this  case.  Although  this  is  somewhat  less  than  the  7.8%  savings  seen  in  the 
base  case,  the  absolute  savings  is  larger. 

7.  Multiple  Chiller  Plant  Analysis 

Most  large  central  plants  have  more  than  one  chiller,  so  multiple  chiller  plant  analysis  is  critical. 
This  case  compares  the  two  plant  designs,  each  with  two  chillers.  Atlanta  weather  data  and  load  profile  are 
used.  As  discussed  earlier,  chillers  are  scheduled  to  minimize  total  system  electrical  energy  requirements. 
Optimal  scheduling  for  the  single  loop  system  is  almost  identical  to  that  of  the  primary-secondary  system. 
In  this  case  proper  scheduling  requires  starting  up  the  lag  chiller  when  the  lead  chiller  is  delivering  97%  of 
its  rated  load. 


System  type 

Single  loop 

Primary-secondary 

Evaporator  Auxiliary 
lER 

0.025  [kW/ton] 

0.044  [kW/ton] 

Chiller  lER 

0.38  [kW/ton] 

0.38  [kW/ton] 

System  lER 

0.51  [kW/ton] 

0.53  [kW/ton] 

Table  28.  Multiple  chiller  plant  lER  analysis 

It  can  be  seen  that  the  chiller  lER  values  for  the  two  chiller  system  are  larger  than  those  of  single 


chiller  systems.  This  is  because  chiller  scheduling  was  based  on  overall  system  power  rather  than  just 
compressor  power  analysis.  Auxiliary  and  total  system  BER  values  are  lower  here  than  in  single  chiller 
results.  This  is  because  the  chillers  in  these  multiple  chiller  systems  are  operated  at  a  higher  average  load 
than  the  single  chiller  systems.  In  this  example,  the  average  load  is  1382  tons  when  the  temperature  is 
greater  than  73^.  Two  chillers  are  operated  when  above  73®F,  so  the  average  load  is  69%.  When  the 
temperature  falls  below  73®F,  the  average  load  is  603  tons,  but  only  one  chiller  is  used.  Thus,  this  one 
chiller  operates  at  about  60%.  As  a  result,  the  average  annual  chiller  load  in  this  multiple  chiller  system  is 
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63%.  Since  the  average  load  is  higher  the  energy  savings  available  by  use  of  the  single  loop  system  is 
somewhat  smaller.  The  savings  here  is  3.9%. 

Over-desian  Load  Operation 

A  chilled  water  system’s  capacity  is  determined  by  a  number  of  factors.  These  include  chiller 
size,  maximum  compressor  motor  power,  maximum  water  flow  rate  and  maximum  air  flow  rate.  In  both 
the  single  loop  and  primary-secondary  systems  the  maximum  compressor  power  is  taken  to  be  588kW. 

This  value  was  acquired  from  the  chiller’s  manufacturer.  The  maximum  system  water  flow  rate  is  taken  as 
2800  gpm.  This  flow  rate  was  chosen  to  be  about  15%  greater  than  the  design  value.  Finally,  the 
maximum  air  flow  rate  into  an  air  handler  is  400  cfm  per  design  ton.  For  example,  an  air  handler  with  a 
100  ton  capacity  will  have  fans  capable  of  a  40,000  cfm  air  flow  rate.  This  is  the  standard  industry 
accepted  value  for  air  flow  rates. 

The  assumed  control  strategy  is  important  when  analyzing  its  capacity.  The  system  constantly 
adjusts  to  changes  attempting  to  meet  the  building  load.  When  a  room  starts  to  heat  up  the  air  handler  fans 
increase  speed  to  deliver  more  cool  air.  As  the  air  flow  rate  increases,  the  cold  deck  temperature  starts  to 
rise.  This  triggers  the  two-way  valves  to  open,  calling  for  more  chilled  water.  As  two-way  valves  open  the 
chilled  water  coil  head  drops,  sending  a  signal  for  the  pumps  to  increase  speed.  As  the  system  flow  rate 
increases  the  chiller  water  exit  temperature  increases  momentarily.  This  increase  signals  the  compressor  to 
increase  capacity,  cooling  the  water  to  the  chiller  set  point  again.  This  sequence  works  well  as  long  as  the 
system  is  not  above  design  load. 

1.  Base  Case 

As  the  actual  load  moves  above  the  design  value,  the  single  loop  and  primary-secondary  systems 
operate  much  differently.  The  primary  loop  of  a  primary-secondary  system  is  constrained  by  a  constant 
2400  gpm  flow  rate.  When  secondary  flow  rate  exceeds  2400  gpm  the  return  water  is  mixed  with  the 
supply  through  the  bypass. 
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Constant  speed 


pump 


Primary-secondary  system 


Figure  57.  Mixing  on  over-design  primary-secondary  system 

This  mixing  increases  the  air  handler’s  water  inlet  temperature,  which  reduces  the  air  to  water-side  LMTD 
and  with  it  the  heat  transfer.  The  air  flow  rate  increases  to  compensate,  but  this  lifts  the  cold  deck 
temperature.  The  air  handler  calls  for  more  water,  attempting  to  drop  the  cold  deck,  further  increasing  the 
inlet  temperature.  This  sequence  progresses  until  the  secondary  loop  reaches  its  maximum  flow  rate. 
Figure  58  shows  how  the  air  handler  water  inlet  temperature  increases  with  the  secondary  loop  flow  rate. 
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Secondary  loop  flow  rate  [gpm] 

Figure  58.  Over-loaded  primary-secondary  system 

In  this  figure  the  chiller  is  able  to  hold  the  44^  set  point  because  the  compressor  is  operating 
below  588kW.  The  building  load  is  met  because  the  air  handler  fans  are  running  below  maximum  levels. 
However,  the  air  handler  cold  deck  temperature  is  no  longer  maintained  at  the  constant  58®F.  Figure  59 
shows  how  the  primary-secondary  mixing  adversely  affects  system  load.  The  negative  slope  on  the 
primary-secondary  system  capacity  curve  relates  the  fact  that  this  is  an  unstable  situation.  As  soon  as  the 
secondary  flow  exceeds  primary  loop  flow  rate  the  system  will  proceed  to  a  “wide  open”  position.  The 
time  that  it  takes  this  process  to  proceed  depends  on  chilled  water  system  length,  water  velocity  and  air 
handler/pump  control  scheme.  Typically  these  controls  are  adjusted  to  be  slow  enough  that  they  do  not 
hunt  for  an  adequate  operating  point  but  fast  enough  to  adjust  to  changing  demand. 
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System  water  flow  rate  [gpm] 


Figure  59.  System  capacity  versus  water  flow  rate 


The  single  loop  system,  having  no  mixing,  maintains  a  positive  slope  throughout  the  flow  rate 
variation.  The  initial  steep  slope  represents  the  area  where  the  chiller  maintains  the  44^  set  point.  At 
approximately  2650  gpm  the  compressor  reaches  the  maximum  power  level.  The  set  point  then  rises  while 
the  compressor  continually  operates  at  maximum  power.  Notice  the  single  loop  system  curve  maintains  a 
positive  slope  after  the  chiller  reaches  full  capacity.  Since  the  chiller  is  delivering  a  warmer  supply 
temperature,  the  chilled  water  return  stream  is  warmer  as  well.  This  decreases  the  compressor  lift,  which 
raises  chiller  capacity.  Chiller  capacity  is  increased  in  the  primary-secondary  system  as  well,  but  it  is 
primary  loop  flow  rate  not  chiller  capacity  that  constrains  that  design.  The  single  loop  design  capacity  is 
constrained  by  chiller  capacity.  Thus,  increased  chiller  capacity  directly  increases  system  capacity. 

Throughout  this  process  the  system  is  attempting  to  meet  the  building  load.  Even  though  the 
chilled  water  has  reached  its  maximum  flow  rate  the  system  may  still  be  able  to  meet  the  load.  The  air 
handler  fans  increase  the  air  velocity  across  the  chilled  water  coil.  This  increases  the  airside  heat  transfer 
coefficient,  resulting  in  a  larger  UA  and  LMTD.  The  affect  of  air  speed  on  system  capacity  is  seen  in 
Figure  60.  This  graph  assumes  that  the  water  flow  rate  is  fixed  at  the  maximum  value  of  2800  gpm.  Also 
the  ambient  temperatures  are  set  at  design  weather  conditions.  The  primary-secondary  system  capacity 
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increases  steeply  until  the  air  flow  equals  33,000  cubic  feet  per  minute.  That  is  the  point  that  the  primary- 
secondary  system’s  compressor  reaches  maximum  load.  Prior  to  that  point  the  chiller  maintains  its  set 
point.  Since  the  water  flow  rate  is  fixed  at  2800  gpm  the  air  handler  inlet  temperature  is  fixed. 


Air  handler  face  velocity  [cfm] 


Figure  60.  System  capacity  versus  air  flow  rate 

The  single  loop  system’s  chiller  is  already  at  full  power,  thus  the  set  point  is  increasing.  Once  the 
compressor  is  are  operating  at  maximum  power  in  each  the  curves  become  parallel.  The  primary¬ 
secondary  system  continues  to  exhibit  lesser  capacity  because  its  air  handler  inlet  water  temperature  is 
increased  by  mixing.  Eventually,  both  systems  reach  their  maximum  capacity.  This  is  the  point  where  the 
chiller,  pumps  and  fans  are  all  running  at  maximum  power. 
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At  maximum  capacity,  the  primary-secondary  system’s  air  handler  water  inlet  temperature  is  greater  than 
the  corresponding  temperature  in  a  single  loop  design.  As  a  result,  the  primary-secondary  air  to  water 
LMTD  is  smaller.  The  chilled  water  coil  UAs  are  identical  because  the  designs  have  equal  air  and  water 
velocities.  Figure  61  represents  each  system  operating  “wide  open”.  The  smaller  air  handler  LMTD  limits 
the  primary-secondary  system  capacity.  Both  systems  exhibit  increasing  temperature  differences.  This 
shows  that  the  air  inlet  temperature  is  raising  with  outdoor  conditions  faster  than  the  water  inlet 
temperature.  Figure  62  shows  how  the  system  capacities  compare. 
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The  single  loop  system  is  able  to  meet  the  building  load  requirement  to  higher  outdoor 
temperatures.  Additionally,  this  system  has  more  capacity  when  it  is  no  longer  able  to  meet  the  actual  load. 

2.  Reducing  Chilled  Water  Set  Point 

One  might  reason  that  reducing  the  chiller  set  point  has  the  potential  to  increase  primary¬ 
secondary  design’s  capacity.  Initially,  the  system  is  constrained  by  primary  flow  rate  and  not  chiller 
capacity.  When  the  set  point  is  reduced  the  chiller  becomes  more  fully  loaded.  Also,  the  air  handler  water 
inlet  temperature  is  lowered.  This  increases  the  LMTD  and  thus  the  heat  transfer. 
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The  problem  with  this  reasoning  is  the  fact  that  the  set  point  is  no  longer  maintained  once  the  compressor 
reaches  maximum  power.  Figure  63  shows  the  effect  that  set  point  reduction  has  on  system  capacity.  Each 
system  is  able  to  meet  the  load  requirement  until  long  after  the  compressor  is  at  full  power.  As  a  result,  the 
system  capacity  remains  unchanged  from  the  original  performance  with  a  44®F  set  point. 

While  this  strategy  fails  to  increase  system  capacity,  it  hurts  system  efficiency.  Figure  64  shows 
how  chiller  efficiency  is  affected  by  reducing  the  evaporator  outlet  temperature.  Initially,  all  three  systems 
are  meeting  the  load.  However,  the  system  with  the  40®  set  point  is  the  most  inefficient.  Both  primary- 
secondary  systems  have  equal  efficiency  once  the  compressors  reach  maximum  power.  The  systems  are 
essentially  identical  thereafter  because  the  set  point  is  no  longer  maintained.  When  ail  three  systems  move 
to  maximum  capacity  the  single  loop  chiller  operates  most  efficiently.  This  is  because  has  a  higher  cooling 
capacity  using  the  same  power. 
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Figure  64.  Chiller  efficiency  with  reduced  set  point 

It  is  able  to  do  so  because  there  is  no  mixing.  Consequently,  the  evaporator  outlet  temperature  can  be 
higher  than  the  primary-secondary  system’s  temperature.  This  reduces  the  compressor  lift,  increasing  the 
chiller  capacity. 


3.  Effects  of  Air  Handler  Degradation  on  System  Capacity 

This  case  compares  the  capacity  of  a  single  loop  and  primary-secondary  system  with  fouled  air 
handlers.  The  air  handlers  are  fouled  to  the  point  that  the  chilled  water  system  AT  is  reduced  to  7®F.  As  a 
result,  system  pumps  are  operating  at  maximum  flow  rate  with  loads  considerably  smaller  than  the  design 
value.  However,  both  systems  are  still  able  to  meet  the  design  load  requirement  by  increasing  the  air 
handler  face  velocity.  Air  handler  fouling  does  not  only  affect  system  efficiency,  it  hurts  system  capacity 
as  well.  The  air  handler  UA  is  reduced  due  to  fouling,  thus  at  maximum  air  and  water  flow  rates  the  air 
handler  is  not  able  to  deliver  the  same  cooling  load.  In  this  case  air  hander  fouling  results  in  a  2%  capacity 
decrease  when  compared  with  similar  unfouled  designs.  Figure  65  shows  that  the  single  loop  system  is  still 
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able  to  deliver  a  larger  cooling  load  than  the  primary-secondary  system.  In  fact,  the  single  loop 
outperforms  the  primary-secondary  design  by  the  same  percentage  as  it  did  in  the  unfouled  analysis. 


4.  Increasing  Air  Handler  Capacity 

The  final  area  analyzed  is  the  affect  that  increasing  the  air  handler  sizes.  In  this  case  it  was 
assumed  that  the  air  handers  were  the  only  components  altered.  Therefore,  the  maximum  compressor 
power  and  water  flow  rate  remain  588kW  and  2800  gpm  respectively.  One  would  expect  bigger  fans  to 
accompany  the  larger  chilled  water  coil.  Due  to  the  larger  air  handlers  the  design  system  AT  is  increased  to 
12°F.  By  doing  so  a  145  ton  air  handler  is  needed.  Using  the  400  cfm  per  design  ton  value,  the  maximum 
air  flow  rate  is  58,000  cfm.  This  air  handler  was  given  greater  capacity  by  simply  increasing  the  face  area. 
Thus  the  larger  air  flow  rate  converts  to  a  similar  face  velocity. 
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Outdoor  Temperature  [F] 

Figure  66.  Increased  design  AT  system  capacity 

Both  system  capacities  are  larger  here  than  in  previous  examples.  This  is  because  the  air  handler 
UA  is  substantially  increased.  The  Reynolds  number  and  heat  transfer  coefficients  are  similar  to  their 
previous  values  because  the  maximum  water  and  air  velocities  remain  unchanged.  The  UA  is  raised 
because  the  air  and  water-side  heat  transfer  area  is  greater.  With  a  larger  UA  the  air  to  water  LMTD  may 
reduce.  As  a  result,  the  evaporator  exit  temperature  may  increase,  reducing  the  compressor  lift.  The 
reduced  compressor  lift  gives  the  chiller  more  capacity.  Figure  66  shows  that  the  single  loop  system  still 
outperforms  the  primary-secondary  system  by  the  same  percentage  it  did  previously. 
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CHAPTER  Vir 

CONCLUSIONS 

This  study  has  analyzed  chilled  water  system  performance  from  two  perspectives.  Using  a  variety 
of  load  profiles  and  system  designs  annual  efficiency  was  compared.  A  number  control  strategies  were 
then  compared  to  determine  a  system’s  ability  to  meet  load  demands  when  they  exceed  the  design  load. 

Air  handler  fouling  was  introduced  to  determine  each  design’s  ability  to  adapt  to  low  chilled  water 
temperature  differential,  a  common  chilled  water  system  ailment.  In  each  case  analyzed  the  single  loop 
system  outperforms  the  corresponding  primary-secondary  design. 

Part  Load  efficiency 

Single  loop,  variable  flow,  chilled  water  system  efficiency  exceeds  primary-secondary  loop 
efficiency  under  a  variety  of  conditions.  These  situations  include: 

1.  Single  and  multiple  chiller  systems 

2.  Diverse  load  profiles 

3.  Constant  and  variable  speed  pumping  schemes 

4.  Different  design  temperature  differentials 

Greater  Capacity 

Single  loop  systems  are  capable  of  meeting  larger  load  demands  at  over  design  conditions, 
because  there  is  no  mixing  from  the  retium  to  the  supply  chilled  water  across  the  bypass  leg.  Various 
design  and  control  strategies  were  studied  to  determine  their  ability  to  increase  primary-secondary  loop 
capacity  to  that  of  the  single  loop  design.  These  include: 

1.  Reducing  the  chiller  set  point 

2.  Increasing  air  handler  capacity 

3.  Inserting  a  check  valve  in  the  chiller  bypass 
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The  first  two  strategies  fail  to  close  the  performance  gap  between  single  loop  and  primary¬ 
secondary  system  capacity.  Including  a  check  valve  allows  the  primary-secondary  loop  to  have  equal 
capacity  but  fails  to  increase  annual  system  efficiency. 

Adapts  to  Load  Problems 

Single  loop  systems  are  more  efficient  and  have  greater  capacity  when  the  chilled  water  AT  is  low. 
Low  chilled  water  AT  is  a  common  problem  among  chilled  water  systems  caused  by  air  handler  problems. 
These  deficiencies  include  fouling,  uncalibrated  sensors  and  poor  air  handler  control.  The  single  loop 
system  is  able  to  adapt  to  these  conditions  by  increasing  chilled  water  flow  rate  through  evaporators  with 
little  adverse  performance  affects. 


Minimum  Evaporator  Water  Flow  Rate 

Evaporator  water  flow  rate  must  be  high  enough  for  the  flow  to  remain  turbulent.  This  is  often 
thought  to  be  a  shortcoming  of  the  single  loop,  variable  flow  rate  design.  Evaporator  flow  rate  is  not  an 
issue  because  the  problem  is  resolved  when  chilled  water  reset  is  employed.  Reset  not  only  ensures 
turbulent  flow  in  the  evaporator,  but  saves  energy  as  well.  Consequently,  reset  is  a  common  practice  used 
with  primary-secondary  systems  today.  Multiple  chiller  use  also  resolves  this  issue. 


Control  Strategy 

The  single  loop,  variable  flow  system  is  much  simpler  to  operate  that  the  primary-secondary 
design.  The  single  loop  system  eliminates  the  possibility  of  detrimental  mixing.  Pumps  are  controlled  to 
maintain  a  constant  pressure  drop  across  a  distant  chilled  water  coil,  while  the  chillers  are  scheduled  as  a 
function  of  cooling  load. 
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Summary 

The  single  loop  design  with  variable  speed  pumps  is  advantageous  to  the  popular  primary¬ 
secondary  system.  This  proposed  design  change  does  not  suggest  that  it  is  better  to  change  to  an  older 
primary  system  with  constant  speed  pumps  and  three-way  valves.  Rather,  this  change  is  the  next  logical 
step  in  chilled  water  system  design  advancement.  Primary-secondary  systems  were  designed  in  order  to 
take  advantage  of  the  reduced  pressure  drop  resulting  from  variable  water  flow  rates  in  the  air  handlers. 
The  single  loop  design  extends  this  reasoning  to  varying  water  flow  rate  in  the  evaporator,  further  reducing 
system  pressure  drop. 
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CHAPTER  Vill 

RECOMMENDATIONS 

This  study’s  results  should  be  verified  by  analyzing  a  sufficiently  instrumented  facility.  This 
instrumentation  must  be  able  to  measure  small  building  load  variations.  Additionally,  water  flow  rate  and 
temperature  readings  are  required  throughout  the  system.  This  exceptional  instrumentation  is  necessary  to 
ensure  a  similar  load  and  proper  control  strategy  in  each  case.  The  system’s  performance  should  be 
measured  with  the  bypass  loop  open  and  closed  under  similar  conditions.  These  conditions  should  include 
operation  above  and  below  design  levels. 

Additionally,  this  model  should  be  revised  as  more  efficient  chilled  water  system  components  are 
developed.  The  results  shown  in  this  study  only  reflect  system  performance  with  components  available 
today.  As  component  performance  changes  the  system  performance  will  be  affected.  This  may  affect  the 
comparative  analysis  presented  here,  but  not  the  advantages  of  the  single  loop  design  over  the  primary¬ 
secondary  system. 

Multiple  chiller  analysis  should  be  extended  to  include  three  and  four  chiller  systems.  Low 
system  AT  complicates  primary-secondary  chiller  scheduling  to  a  greater  degree  when  more  chillers  are 
involved.  Single  loop  chiller  scheduling  is  a  function  of  load,  while  primary-secondary  scheduling  is  tied 
to  system  flow  rate.  These  effect  of  these  different  strategies  would  be  pronounced  by  adding  chillers  to 
the  systems. 


APPENDIX  A 


BASE  CASE,  SINGLE  LOOP  MODEL 


"Base  Case,  Single  Loop,  Chilled  Water  System" 

{Section  I.  Procedures) 

(A.  Air  Handler  Procedures) 

{This  PROCEDURE  determines  the  load  based  on  out  door  dry  bulb  temperature) 

Procedure  loadprofile(T_amb:Q_ton_ah) 

IF  (T_amb<=60)  Then 

Q_ton_ah=20 
Endif 

IF  (T_amb>60)  Al«)  (T_amb<=105)  Then 

Q_ton_ah=I6/7*T_amb-l  17.143 
Endif 

IF  (T_amb>105)  Then 

Q_ton_ah=100 
Endif 
END 

{This  procedure  determines  the  flow  regime  and  assigns  a  Nusselt  number  and  friction  factor) 
Procedure  regime(Re_D_i_ah4*r_w_ah,D_i_ah,epsilon_pipe_ah:Nus_ah,f_i_ah) 

IF  (Re_D_i_ah<=2000)  Then 

Nus_ah=4.36  {eqn3.65  ref[l]) 

fJ_ah=64/Re_DJ_ah  {eqn3.29  refll]) 

Endif 

IF  (Re_D_i_ah>2000)  Then 

Nus_ah=0.023*Re_D_i_ah'H).8*Pr_w_ahN).4  {eqn  1 1.16  refl2]) 

f_i_ah=  l/(-l  .8*loglO(6.9/Re_D_i_ah+(epsilon_pipe_ah/(D_i_ah*3.7))''  1 . 1 1))''2 

{eqn  6.64a  ref[5] ) 

Endif 

END 

{ This  PROCEDURE  determines  the  relative  heat  capacities ) 

Procedure  heat(C_a_ah,C_w_ah:C_min_ah,C_max_ah) 

IF  (C_a_ah/C_w_ah<=l)  THEN 
C_min_ah=C_a_ah 
C_max_ah=C_w_ah 
Endif 

IF  (C_a_ah/C_w_ah>l)  THEN 
C_min_ah=C_w_ah 
C_max_ah=C_a_ah 
ENDif 
End 


{This  procedure  determines  the  heat  exchanger  efficiency  given  that  the  air  flow  is 
unmixed  and  the  water  is  mixed} 

Procedure  transfer(C_a_ah,C_w_ah,NTU_ah,C_r_ah:epsilon_ah) 

IF  (C^a^ah/C^w^ah>l)  THEN 

epsilon_ah=l-exp((-l/C_r_ah)*(l-exp(-C_r_ah*NTU_ah)))  (eqn  1 1.35a  ref[3] } 

Endif 

IF  (C_a^ah/C^w^ah<=l)  THEN 

epsilon_ah=(l/C_r_ah)*(l-exp(-C_r_ah*(l-exp(-NTU_ah))))  (eqn  1 1.34a  ref[3] } 

Endif 

END 

(This  procedure  determines  the  f  for  air  flow  over  bare  tubes. } 

Procedure  baretube(Re_D_o_ah:f_tube_ah) 

IF(Re^D^o_.ah<1000)  THEN 

Uube_ah=3.680886*Re„D_o^ah^-0.3090731)  (Figure  7.13  ref[l]} 

Endif 

IF  (Re_D_o_ah>=i(X)0)  THEN 

Ltube_ah=0.42275-5.53e-6*Re__D_o_ah-3.598e-10*Re_D„o_ah''2{Figure  7.13  ref[i] } 
Endif 
END 

(This  procedure  determines  the  effective  specific  heat  for  air  of  the  air  handler) 

Procedure  latentheat(T_a_ahJn,T_a_ah_out,T_a_ah_ave,D_a_ah_in,HMRT_a_ahJn,P_atm, 
P_a_ah_in,h_a_ah_in;HMRT_a_ah_out,h_a_ah_out,h_a_ah_out_eff,Cp_eff) 

IF  (T_a_ah_out>D_a_ah_in)  Then 
HMRT_a_ah_out=HMRT_a_ah_in 

h_a_ah_out=enthalpy(AIRH20,T=T_a_ah_out,P=P_atm,W=HMRT_a_ah_out) 

h_a_ah_out_eff=h_a_ah_out 

Cp_eff=specheat(AIRH20,T=T_a_ah_ave,P=(P_a_ahJn+P_atm)/2,D=D_a_ahJn) 

Endif 

IF  (T_a_ah_out=<D_a_ah_in)  Then 

HMRT_a_ah_out=humrat(AIRH20,T=T_a_ah_outJ^=P_atm,R=l) 
h_a_ah_out=enthalpy(  AIRH20,T=T_a_ah_out,P=P_atm,R=  1 ) 
h_a_ah_out_eff=h_a_ah_out+(HMRT_a_ahJn- 
HMRT_a_ah_out)*entha!py(WATERJ*=P_atm,T=T_a_ah_out) 
Cp_eff=(h_a_ahJn-h_a_ah_out_eff)/(T_a_ahJn-T_a_ah_out) 

Endif 

END 
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{B.  Chilled  water  reset  procedure} 

{This  procedure  first  sets  mass  flow  rate  as  a  function  of  temperature,  then  at  a  specified  outdoor  dry  bulb, 
it  switches  to  make  the  evaporator  exit  temperature  a  function  of  flow  rate. } 

Procedure 

reset(T_amb,Q_dot_evap,h_w_evap_in,h_w_evap_out,P_w_evap_out:m_dot_w_evap,T_w_evap_out) 

IF  (T^amb=>77)  Then 
T_w_evap_out=44 

m_dot_w_evap=Q_dot_evap/(h_w_evap_in-h_w_evap_out) 

Endif 

IF  (T^amb<77)  Then 
ni_dot_w_evap=538300 

T_w_evap_out=temperature(WATER,h=(h_w_evapJn-Q_dot_evap/m_dot_w_evap) 

,P=P_w_evap_out) 

Endif 

END 

( C.  Tower  procedures. } 

(This  procedure  determines  the  tower  fan  duty} 

Procedure 

to  werfan(Q_dot_to  w,epsilon_to  w,m_dot_a_to  w,h_s_w_to  w Jn,h_a_to  w Jn,T_w_to  w_out  1 : 
fanduty,T_w__tow_out) 

IF  (T_w_tow_outl>70)  Then 

fanduty=l 

T_w_tow_out=T_w_tow_outl 

Endif 

IF  (T_w_to  w_out  1<70)  AND  ((Q„dot__tow-0. 1  *epsilon_tow*m_dot_a_tow*(h_s_w_tow  Jn- 

h_a_towjn))/(0.9*epsilon_tow*m_dot_a__tow*(h_s_w_towjn-h_a_towjn))>0)  Then 

fanduty=(Q_dot_tow-0. 1  *epsilon_tow*m_dot_a_tow*(h_s_w__tow_in- 
h_a_tow__in))/(0.9*epsilon_tow*m__dot_a_tow*(h_s_w_tow_in-h_a_tow_in)) 

T_w_tow_out=70 

Endif 

IF  (T_w_tow_outl<70)  AND  ((Q_dot_tow-0. 1  *epsilon_tow*m_dot_a_tow*(h_s_w_tow  Jn- 

h_a_tow_in))/(0.9*epsilon_tow*m_dot_a_tow*(h_s_w_tow_in-h_a_tow_in))<=0)  Then 
fanduty=0 

T_w_to  w_out=T_w_to  w_out  1 
Endif 
END 

{Section  11.  Inputs} 

Hours=l  {hrs} 

T__amb=95  {F} 

ODWB=78  {F} 

{ Section  III.  Overall  Calculations } 

{A.  Energy  analysis} 

W_dot__total=W_dot_comp+W_dot_pump+W_dot_pump_cond+  10*W_dot_fan+W_dot__fan_to w  { btu/hr } 
KW_total=W_dot_total*convert(btu/hr,kW)  {kW} 

KWh_total=KWjotal*Hours  {kWh} 
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{These  calculations  sum  to  the  inverse  of  the  lER} 
IER_chiller=((Hoiirs/Hours_total)/(KW_comp/Q_ton_evap)) 

IER_sy  stem=(Hours/Hours_total)/(kW_total/(  1 0*Q_ton_ah)) 

IER_aux_evap=(Hours/Hours_total)/((W_dot_pump+10*W_dot_fan)*convert(btu/hr,kW)/(10*Q_ton_ah)) 

Hours_total=7l32 

(B.  Entropy  Analysis} 

S_dot_gen_fan_ah=10*m_dot_a_ah*(s_a_fan_out-s_a_fanjn)  {btu/hrR} 

S_dot_gen_ah=  1 0*(m_dot_w_ah  *(s_w_ah_out-s_  w_ah_in)-Q Joad_ah/(T_amb+459.67))  { btu/hr  R } 
S_dot_gen_pump=m_dot_w_pump*(s__pump_out-s_pump_in)  {btu/hrR} 

S_dot_gen_evap=m__dot_w_evap*(s_w_evap_out-s_w_evap_in)-(-Q_dot_evap/(T_evap+459.67)) 
S_dot_gen_cond=m_dot_w_cond*(s_w_cond_out-s_w_condJn)-(Q_dot_cond/(T_cond+459.67)) 
S_dot_gen_comp=(W_dot_comp-W_dot_carnot)/T_amb  { btu/hr  R } 

S_dot_gen_tower=m__dot_w_tow*(s_w_tow_out-s_w_tow_in)-(-Q_dot_tow/(T_w_tow_ave+459.67)) 
S_dot_gen_pump_cond=m_doCw_cond*(s_pump_cond_out-s_pump_cond_in)  { btu/hr  R } 

S_dot_gen__loop_pri=m_dot_w__pump*(s_pump_in-s_w_ah_out+s_w_ahJn“S_w_evap_out) 
S_dot_gen_loop_cond=m_dot_w_cond*(s__pump_cond_in-s_w_tow__out+s_w_towJn-s_w_cond_out) 
S_dot_gen_tot=S_dot_gen„fan_ah+S_dot_gen_ah+S_dot_gen_pump+S_dot_gen_evap+S_dot_gen_cond+ 
S_dot_gen__comp+S_dot_gen_tower+S_dot_gen_pump_cond+S_dot_gen_loop_pri+S_dot_genJoop_cond 
W_dotJost=S_dot_^en_tot*T_amb  {btu/hr} 


{Section  IV.  Air  handler  model } 

{A.  Air  handier  fan} 

{1.  Fluid  properties} 

T_room=75  {F} 

T_a„fanJn=temperature(AIRH20,h=:h_a__fanJnJ*=P_atm,W=HMRT_fan__in)  {F} 
T__a_fan_out=temperature(AIRH20,h=h_a_fan_out,P=P_atm+DELTAP_a_ah,W=HMRT_fan_in) 
HMRT_amb=humrat(AIRH20,T=T_amb,P=P_atm,B=ODWB) 
HMRT_room=humrat(AIRH20,T=T_room,P=P_atm,R=R_room) 

HMRT_fan Jn=0.85*HMRT_room+0. 15*HMRT_amb 
R_room=.5549 

h_amb=enthalpy(AIRH20,T=T_ambJ^=P_atm,B=ODWB)  {btu/lbm} 

h_room=enthalpy(AIRH20,T=T_room+54’=P_atm,R=R_room)  {btu/lbm} 

h_a_fanjn=0.85*h_room+0.15*h_amb  {btu/lbm} 

h_a_fan_out_rev=enthalpy(AIRH20J^=P_atm+DELTAP_a_ah,T=T_a_fan_out_rev,W=HMRT_fan_in) 
s_a_fan_in=entropy(  AIWI20,T=T_a_fan_in  JP=P_atm,W=HMRT_fan Jn)  { btu/lbm  R } 

s_a_fan_in=entropy(AIRH20,T=T_a_fan_out_revJP=P_atm+DELTAP_a__ah,W=HMRT_fan_in) 
s_a_fan_out=entropy(AIRH20,T=T_a_fan_out,P=P_atm+DELTAP_a_ah,W=HMRT_fanJn) 

{2.  Fan  Calculations} 
eta_fan=.65 

eta_fan=(h_a_fan_in-h_a_fan_out_rev)/(h_a_fan_in-h_a_fan_out) 
W_dot_fan=-m_dot_a_ah*(h_a_fan_in-h_a_fan_out)  { btu/hr } 

KW_fan=W_dot_fan*convert(btu/hr,kW)  {kW} 
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{B.  Inputs  for  Air  Handler) 

h_w_ahJn=h_w_cvap_out  {btu/Ibm} 

h_a_fan_out=h_a_ahJn  {btu/lbm} 

P_a_ahJn=P_atm+DELTAP_a_ah  {psia} 

P_w_ahJn=P_w_evap_out-(Head_pipe/2+ElevJoad+Head_valve)*rho_w_ah*g*convert(lbm/ft*s'^2,psia) 
Head_valve=20-Head_coil  {ft} 

HMRT_a_ahJn=HMRT_fanJn 

(C.  Air  Handler  overall  heat  transfer  coefficient  calculation) 

l/UA_tot_ah=l/(h_Lah*AJ_ah)+ln(D_o_ah/DJ_ah)/(2*pi*k_pipe*Length_ah*n_row*n_column) 
+l/(h_o_ah* A_o_ah*eta__o_ah)  (eqn  1 1 .5  ref[3] } 


{D.  Geometry) 

Length^ah=7.555  (ft) 

Height_ah=5  (ft) 

Depth_ah=8.5*convert(in,fl)  {ft} 

fin_pitch=l  2/con  vert(in, ft)  { l/ft) 

D_o_ah=.75*convert(in,ft)  { ft } 

D_Lah=.666*convert(in,ft)  {ft} 

t_fin=.006*convert(in,ft)  {ft} 

n_row=4 
n_column=40 

Elevjoad=0  {ft} 


{E.  Airside  calculations} 

{1.  Air  properties} 

Call  latentheat(T_a_ah_in,T_a_ah_out,T_a_ah_ave,D_a_ahJn,HMRT_a_ah_in,P_atm, 

P_a_ahJn,h_a_ahJn:HMRT_a_ah_out,h_a_ah_out,h_a_ah_out_eff,Cp__eff) 

HMRT_a_ah_ave=(HMRT_a_ahJn+HMRT_a_ah_out)/2 

k_a_ah=conductivity(AIRH20,T=T_a_ah_ave,P=P_a_ah_ave,W=HMRT_a_ah_ave) 

mu_a_ah=viscosity(AIRH20,T=T_a_ah__aveJP=P_a_ah_ave,W=HMRT_a_ah_ave) 

rho_a_ah=density(AIRH20,T=T_a_ah_aveJ*=P_a_ah_ave,W=HMRT_a_ah_ave) 

v_m_ah=  l/rho_a_ah 

Pr_a_ah=Cp_eff*mu_a_ah/k_a__ah 


{Btu/hr-ft-R} 

{Ibm/fMir} 

{lbm/ft^3} 

{ft^3/lbm} 


{2.  Air  flow  properties} 

G_a_ah=m_dot_a_ah/(sigma_ah* A_fr_ah)  { lbni/ft'^2-hr }  { eqn  11.37  ref[3] } 

Re_D_o_ah=G_a_ah*D_o_ah/mu_a_ah  {eqnl4.37  ref[4] } 

Re_L_ah=G__a_ah*b_ah/mu_a„ah  { eqn  ref[4] } 


{3.  Airside  geometry  calculations} 

a_ah=Height_ah/n_colunm  { ft } 

b__ah=Depth_ah/n_row  {ft} 

A_fi-_ah=Height_ah*Length_ah  {ft'^2} 

sigma_ah=(a_ah-D_o_ah)/a_ah 

A_f_ah=(2*a_ah*b_ah-pi*D_o_ah''2/2)*fin_pitch*Length_ah*n_row*n_coIumn  { ft'^2 } 

A_t_ah=pi*D_o_ah*Length_ah*n_row*n_column  {ft'^2} 

A_o_ah=A_f.ah+A_t_ah  {ft^2} 


{ 4.  Airside  Heat  Transfer  Calculations } 
h^o^ah=(Cp_eff*G^a_ah*j_n/Pr„a_ah''(2/3))  {Btu/ft^2-hr-R} 

j^4=0.26i8*JP+0.0014 

JP=Re_D_o_ah'^(-0.4)*(A_o_ah/A_t__ah)''(-0.15) 

Ln/j^4=(l-1280*n_row*Re_L_ah''(-L2))/(M280*4*Re_L_ah'^(-1.2)) 


{eqnl4.36  ref[4]} 
{eqn  ref[4]} 

{eqn  14.46  ref[4] 
{eqni4.50  ref[4]} 


{5.  Fin  Calculations} 

eta_o_ah=l-(A_f_ah/A_o_ah)*(l-eta_fin)  {eqn  1 1.3  ref[3] } 

eta_fin=.95 


(F.  Waterside  calculations} 

(1.  Water  properties} 

Cp_w_ah=specheat(WATER,T=T_w_ah_aveJ*=P_atm) 

k_w_ah=conductivity(WATER,T=T_w_ah_ave,P=P_atm) 

mu_w_ah=viscosity(WATER,T=T_w_ah__ave,P=P_atm) 

rho_w_ah=density(WATER,T=T_w_ah_aveJP=P_atm) 

Pr__w_ah=Cp_w_ah*mu_w_ah/k_w_ah 


{Btu/lbm-R} 

{Btu/hr-ft-R} 

{Ibm/ft-hr} 

(lbm/ft'^3} 


{2.  Flow  properties} 

G_w_ah=m_dot_w_ah/((n_column- 1  )*pi*D_i_ah'^2/4)  { lbin/ft^2-hr } 

Re_D_Lah=G_w_ah*D_i_ah/mu_w_ah  { eqn  11.10  ref[2] } 


(3.  Waterside  geometry  calculations} 

A_i_ah=pi*D_i_ah*Length_ah*n_ro  w*n_column  { ft'^2 } 

{ 4.  Waterside  heat  transfer  calculations } 
CallRegime(Re_DJ_ahJ^r_w_ah,DJ_ah,epsilon_pipe_ah:Nus_ah,fJ_ah) 
h_Lah=Nus_ah*k_w_ah/D_i__ah  {Btu/ft'^2-hr-R}  (eqn  3.66  ref[l] } 


(G.  Pipe  and  fm  properties} 
k_pipe=401  *convert(W/m*K,Btu/hr*ft*R) 
k^fin=237*convert(W/m*K.Btu/hr*ft*R) 
epsilon_pipe_ah=.000 15 


{Btu/hr-ft-R} 

{Btu/hr-ft-R} 

{ft} 


{H.  Load  conditions.} 

CALL  loadprofile(T_amb:Q_ton_ah) 

CALLTransfer(C_a_ah,C_w_ah,NTU_ah,C_r_ah:epsilon_ah) 
Call  heat(C_a_ah,C_w_ah:C_min_ah,C_max_ah) 
QJoad_ah=Q_ton_ah*convert(ton,Btu/hr) 
QJoad_ah=:m_dot_w_ah*(h_w_ah_out-h_w_ahJn) 

Q Joad_ah=m_dot_a_ah*Cp_eff*DELT  AT_a_ah 

QJoad_ah=epsilon_ah*C_min_ah*(T_a_ah_in-T_w_ahJn) 

Q_max_ah=C_min_ah*(T_a_ah_in-T_w_ah_in) 

C_w_ah=m_dot_w_ah*cp_w_ah 

C_a_ah=m_dot_a_ah*Cp_eff 

C_r_ah=C_min_ah/C_max_ah 

NTU_ah=UA_tot_ah/C_min_ah 


{Btu/hr} 

{Btu/hr} 

{Btu/hr} 

{Btu/hr} 

{Btu/hr} 

{Btu/hr  R} 

{Btu/hr  R} 


{eqn  11.25  ref[3]} 


{I.  Temperature} 

T_w_ah  Jn=temperature(W  ATER,H=h_  w_ah_in,P=P_w_ah_in)  { F } 

T_w_ah_out=temperature(WATER,h=h_w_ah_out  J^=P_w_ah_out)  {F} 

DELTAT_w_ah=T_w_ah_out-T_w_ahJn  {F} 

DELTAT_a_ah=T_a_ahJn-T_a_ah_out  {F} 

T_w_ah_ave=(T_w_ah_in-T_w_ah_out)/ln(T_w_ah Jn/T_w_ah_out)  {F } 

T_a_ah_ave=(T_a_ah_in-T__a_ah_out)/ln(T_a_ahJn/T_a_ah_out)  {F} 

T_a_ahJn=temperature(AIRH20,h=h_a_ahJnJP=P_a_ahJn,W=HMRT_a_ahJn)  {F} 

T_a_ah__out=58  {F} 

D_a_ah  Jn=dewpoint(  AIRH20,T=T_a_ah  Jn,P=P_a_ah Jn,W=HMRT_a_ah_in)  {F } 


{J.  Entropy} 

s_w__ah_in=entropy(WATER,P=P_w_ah_in,T=T_w_ahJn) 

s_w_ah_out=entropy(WATERJP=P_w_ah_out,T=T_w_ah_out) 

{K.  Flow  rates} 

{1.  Water  flow} 

V_w_ah=FPS_w_ah*(n_column- 1  )*pi  *D_Lah'^2/4*convert(  1/s,  1/min) 
GPM_w_ah=V_w_ah*convert(ft'^3,gal) 
m_dot_w_ah=V__w_ah*rho_w_ah*60 
{2.  Airflow} 

V__a__ah=FPM_a_ah*Length_ah*Height_ah 
m_dot_a_ah=V_a_ah*rho_a_ah*60 
FPH_a_ah=FPM_a_ah*convert(  1/min,  l/hr) 
FPS_a_ah=FPM_a_ah*convert(  1/min,  1/s) 


{btu/lbm  R} 
{btu/lbmR} 


{cfm} 

{gpm} 

{Ibm/hr} 


{cfm} 

(Ibm/hr} 

(ft/hr} 

(ft/s} 


{L.  Pressure  drop  calculations} 

{1.  Waterside  pressure  drop} 

DELT  AP J_ah=rho_w_ah*g*Head_coil*convert(lbm/ft*s'^2,psia)  { psia } 
DELTAP_tube_ah=f_i_ah*Length_ah*n_row/D_i_ah*(rho_w_ah*FPS_w_ah'^2/2) 
*convert(lbin/s'^2*ft,psia)  { psia }  ( eqn  3.27  ref[  1  ] } 

Head_tube_ah=DELTAP_tube_ah/(rho_w_ah*g)*convert(psia*ft'^2*s^2/lbm,ft)  {eqn  6.25  ref[6]} 
Head3end=:5*K*FPS^w_ah'^2/(2*g)  {ft} 

Head_coil=Head_tube_ah+Head_bend  {ft} 

P_w_ah_out=P_w_ah_in-DELTAP_Lah  {psia} 

P_w_ah_ave=(P_w_ah_in+P_w_ah_out)/2  {psia} 

K=.4 


{2.  Airside  pressure  drop} 

CALL  Baretube(Re_D_o_ah:fL.tube_ah) 

DELTAP_a_ah=deltaP_o_ah*convert(lbm/hr'^2*ft,psia)  { psia } 

P_a_ah_ave=(P_a_ah  Jn+P_atm)/2  { psia } 

f_f_ah=(deltaP_o_ah-deltaP_t_ah)/(v_m_ah*G_a_ah'^2/(2*sigma_ah)) 

LLah=1.7*Re_L_ah''(*0.5) 

a_ah_star=a_ah/D_o_ah 

b_ah_star=b_ah/D_o_ah 

ChLah=l 

deltaP_t_ah=n_row*Oube_ah*Chi_ah*(rho_a_ah*(FPH_a_ah/sigma_ah)'^2/2) 


{eqn  7  ref[6]} 
{eqn9ref[6]} 
{eqn7.114ref[l]} 
{eqn  7.114  ref[l]} 
{Fig  7.13  ref{l]} 
{lbm/hr^2  ft} 

{eqn  7.113  ref[l]} 
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{Section  V.  Primary  pump  model } 

(A.  Inputs  for  primary  pump} 

h_pump Jn=h_w_ah_out  { btu/lbm } 

ni_dot_w_pump=  1 0*m_dot_w_ah  { Ibm/hr } 

(B.  Pressure  calculations) 

DELTAP_pump=rho_w_pump*g*Head Joop*convert(lbm/ft*S''2,psia)  { psia } 


P_pump_out=P_min+DELT/^_pump 

{psia} 

P_pump_ave=P_min+DELT  AP_pump/2 

(psia} 

P_min=125 

{psia} 

{C.  Affinity  calculations} 

N_pump_act=N_pump_design*(GPM_pump/GPM_design) 

{rpm} 

{eqn  6.1  ref[7]} 

N_pump__act=N_pump_design*(Head_loop/Head_design)'^(0.5) 

{rpm} 

{eqn  6.2  ref[7]} 

N_pump_design=  1 1 80 

{rpm} 

{ref[8]} 

(D.  Pump  curves} 

Head_design=103.3620  +  0.01053070*GPM^design  -  0.0000063  PGPM_design'^2 
eta^pump=(0.08592879  +  0.00059656*GPM_design  -  0.00000012*GPM_design'^2)* 

(0.8158200  +  0.402540 l*N_pump_act/N_pump_design  -  0.21 8360 l*(N_pump_actW_pump_design)'^2) 
eta_motor=0.0498386+ 1 .846339*(N__pump_act/N_pump_design)- 
1 .029*(N_pump_act/N_pump_design)'^2 

eta_pump=(h_pump_out_rev-h_pumpJn)/(h_pump_out-h_pump_in) 
W_dot_pump=m_dot_w_pump*(h_pump_out-h_pump Jn)/eta_motor  { B  tu/hr } 

KW_pump=W_dot_pump*convert(Btu/^,kW)  { kW } 

(E.  Head  calculation} 

Head_pipe=f_loop*LengthJoop/D_Lpipe*(FPS_pump'^2/2)/g  {ft}  {eqn  6.30  ref[5] } 

Head_loop=Head_pipe+Head_coil+Head_evap+Head_valve  { ft } 

fJoop=l/(-1.8*logl0(6.9/Re_D_Lpump+(epsilon_pipe/(D_Lpipe*3.7))'^l.ll))'^2{eqn  6.64a  ref[5]} 


{F.  Pipe  geometry} 

D J_pipe=  10*convert(in,ft)  {ft} 

Length Joop=2000  {ft] 

epsilon_pipe=.(XX)  15  { ft } 


{G.  Row  properties} 

GPM_pump=m_dot_w_pump/rho_w_pump*convert(ft'^3/hr,gpm) 
CFM_pump=GPM_purnp*convert(gpm,cfm) 
FPM_pump=CFM_pump/(pi*DJ_pipe'^2/4) 
FPS_pump=FPM_pump*con  vert(  1/min,  1/s) 
FPH_pump=FPM_pump*convert(  l/minj/hr) 
Re_D_Lpump=rho_w_pump*FPH_pump*D_i_pipe/mu_w_pump 


(gpm} 

{ft'^3/min} 

{ft/min} 

{ft/s} 

{ft/hr} 


{H.  Ruid  properties} 

rho_w_pump=^ensity(WATER,T=T_pump_ave  J*=P_pump_ave)  { lbm/ft''3 } 

mu_w_pump=viscosity(W  ATER,T=T_pump_ave,P=P_pump_ave)  { Ibm/ft  hr } 


{I.  Temperatures} 

T_pump_in=temperature(WATERJ*=P_inin,h=h_pump_in) 

T_pump_ave=(T_pump_out-T_pump_in)/ln((T_pump_out)/(T_pump_in)) 

T_pump_out_rev=temperature(WATERJ*=Pj)ump_out,s=s_pump_in) 

T_pump_out=temperature(WATER,P=Pj)ump_out,H=h_pump_out) 

{J.  Enthalpy  and  entropy} 

h_pump_out_rev=enthalpy(WATERJ*=P_pump_out,T=T_pump_out_rev) 
s_pump_in=entropy(WATER,P=P_min,T=T_pump_in) 
s_pump_out=entropy(W  ATER  J*=P_pump_out,T=T_pump_out) 

{Section  VI.  Centrifugal  chiller  model} 

(A.  Chiller  inputs} 

{1.  Evaporator  Side} 


P__w_evap_in=P_pump_out 

(psia} 

h_w_evapJn=h_pump_out 

{btu/lbm} 

m_dot_w_evap=m_dot_w_pump 

{Ibm/lbm} 

{2.  Condenser  Side} 

P_w_cond_in=P_pump_cond_out 

{psia} 

h_w_cond_in=h__pump_cond„out 

{btu/lbm} 

{B.  Compressor} 

W_dot_camot=Q_dot_evap*((T_cond+459.67)/(T_evap+459.67)-l)  (btu/hr} 
KW_comp=W_dot_comp*convert(Btu/hr,kW)  {kW} 

W_dot_comp=a+b*W_dot_camot+c*W_dot_camot^2  ' 

a=3042.991  +  705.2175*Q_ton_evap  +  0.9521675*Q_ton_evap^2 
b=l  .228996  -  0.00438602*Q_ton_evap  +  0.00000203*Q_ton_evap''2 
c=0.00000279  -  3.146E-09*Q_ton_evap  +  1.310E-12*Q_ton_evap''2 

(C.  Evaporator} 

(1.  Evaporator  heat  transfer  calculations} 

Call 

Reset(T_amb,Q_dot_evap,h_w_evap_in,h_w_evap_outJ*_w_evap_out:m_dot_w_evap,T_w_evap_out) 
Q_dot_evap=Q_ton_evap’'‘convert(ton,Btu/hr)  (Btu/hr} 

"Q_dot_evap=m_dot_w_evap*(h_w_evap_in-h_w_evap_out)  {Btu/hr} (Blocked  out  for  reset}" 

Q_dot_evap=epsilon_evap*m_dot_w_evap’''C_w_evap*(T_w_evap_in-T_evap)  {Btu/hr} 

{eqnll.l9ref[3]} 

epsilon_evap=l-exp(-NTU_evap)  {eqn  1 1.30a  ref[3] } 

NTU_evap=UA_evap/(m_dot_w_evap*C_w_evap)  {eqn  11.25  ref[3}} 

l/UA_evap=l/(hA_i_evap)+l/(hA_o_evap) 
hAJ_evap=6213*(Re_D_i_evap'H).5334881)*Pr_w_evap'K).3 
hA_o_evap=l  80935.0  +  6243. 157*Q_ton_evap  - 1 .474305*Q_ton_evap''2 

(2.  Evaporator  pressure  drop  calculations } 

f_i_evap=l/(-1.8*logl0(6.9/Re_D_i_evap+(epsilon_pipe_evap/(D_i_evap*3.7))^l.ll))''2 


{eqn  6.64a  ref[5] } 

DELT  AP_evap=Head_evap*rho_w_evap*g*convert(lbm/ft*s''2,psia)  { psia } 

Head_evap=f_i_evap*Length_evap*n_pass_evap/D_i_evap*(FPS_evap''2/2)/g  { ft } 

(eqn  6.30  ref[5} } 

P_evap_ave=P_min+DELTAP_evap/2  (psia} 

P_w_evap_out=P_w_evap_in-DELTAP_evap  {psia} 


{F} 

{F} 

{F} 

{F} 


{Btu/lbm} 

{Btu/lbR} 

{Btu/lbR} 


{3.  Evaporator  geometry} 


D_Levap=0.824*convert(in,ft) 

(ft) 

Length__evap=14.2 

{ft} 

n_tubes_e  vap=2 1 5 

n_pass_evap=4 

epsilon_pipe_evap=.0003873 1 

(ft) 

{4.  Evaporator  fluid  properties} 

C_w_evap=specheat(  W  ATER,T =T_w_evap__ave  JP=P_evap_ave)  { B  tu/lbm  R } 

rho_w_evap=density(WATER,T=T_w_evap_aveJP=P_evap_ave)  { lbm/ft^3 } 
mu_w_evap=viscosity(WATER,T=T_w_evap_ave,P=P_evap_ave)  { Ibm/ft  hr } 
k_w_evap=conductivity(WATER,T=T_w_evap_aveJ^=P_evap_ave)  {Btu/ft  hr  R} 
Pr_,w_evap=C_w_evap*mu_w_evap/k_w__evap 
Re_D_i_evap=rho„w_evap*FPH_evap*DJ_evap/mu_w_evap 


{5.  Evaporator  flow  rate} 

GPM_w_evap=m_dot_w_evap/rho_w_evap*convert(ft'^3/hr,gpm)  { gal/min } 

CFM_w_evap=GPM_w_evap*convert(gpm,cfm)  { ft^3/min } 

FPM_evap=CFM_w_evap/((pi*D_Levap'^2/4)*n_tubes_evap)  { ft/min } 

FPS_evap=FPM_evap*con  vert(  1/min,  1/s)  { ft/s } 

FPH_evap=FPM_evap*convert(  1/min,  1/hr)  { ft/hr } 


{6.  Evaporator  Temperatures} 

"T_w_evap_out=44  { F }  { B  locked  out  for  reset } " 

T_w_evap Jn=temperature(WATER,h=h_w_evap_in  J^=P_w_evap Jn)  { F } 

T_w_evap_ave=(T_w_evap_in-T_w_evap_out)  An(T_w_evap_in/T_w_e  vap^out)  { F } 


{7.  Enthalpy  and  entropy} 

h_w_evap_out=enthalpy(W  ATER,T=T_w_evap_out  J^=:P_w_evap_out)  {btu/lbm} 

s_w_evap__in=entropy(WATER,h=h_w_evap Jn  J*=P_w_evap_in)  { btu/lbm  R } 

s_w_evap_out=entropy(W  ATER,h=h_w_evap_out  JP=P_w_evap_out)  { btu/lbm  R } 


{D.  Condenser} 

{L  Condenser  U A  calculation} 

Q_dot_cond=epsilon_cond*m_dot_w_cond*C_w_cond*(T_cond-T_w_cond__in)  { B  tu/hr } 


Q_dot_cond=m_dot_w_cond*C_w_cond*(T_w_cond_out“T_w_cond_in) 

Q_dot_cond=Q_dot_evap+W_dot_comp 

Q_ton_cond=Q_dot_cond*convert(btu/hr,ton) 

epsilon_cond=  1  -exp(-NTU_cond) 

NTU_cond=UA_cond/(m_dot_w_cond*C_w_cond) 

UA_.cond_l=1845702  +  22792.36*T^cond  -  96.24830*T_cond^2 
UA_cond_2=L059937  -  0.00009226*Q_ton_cond  +  0.00000001 *Qjon_cond'^2 
UA_cond=UA_cond„l  *UA_cond_2 


{eqnl.23ref[3]} 

{Btu/hr} 

{Btu/hr} 

{tons} 

{eqn  11.30ref[3]} 


{2-  Condenser  fluid  properties} 

rho_w_cond=density(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 

C_w_cond=specheat(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 

mu_w_cond=viscosity(WATER,T=T_w_cond__ave,P=P_w_cond_ave) 

k_w_cond=conductivity(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 


{lbm/ft'^3} 
{Btu/lbm  R} 
{Ibm/fthr} 
{Btu/hr  ft  R} 
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{3.  Condenser  Temperatures} 

T_w_cond_ave=(T_w_cond_out-T_w_cond Jn)/ln(T_w_cond_out/T_w_cond Jn)  { F } 

T_w_cond__in=temperature(W  ATER,H=h_w_cond_in,P=P_w_cond Jn)  { F } 

{4.  Enthalpy  and  entropy) 

h_w_cond_out=enthalpy(  WATER, T-T_w_cond_outJ^=P_w_cond_out)  { btu/lbm } 

s_w_cond_in=entropy(WATER,H=:h_w_cond_in,P=P_w_cond_in)  {btu/lbm  R} 

s_w_cond_out=entropy(WATER,T=T_w_cond_outJ^=P_w_cond_out)  (btu/lbm  R} 

(5.  Condenser  pressures} 

Head_cond=60.63713*T_w_cond_ave''(-0.1908161) 

DELT  AP_cond=rho_w_cond*g*Head_cond*convert(lbm/s^2*ft,psia)  { psia } 

P_w_cond_out=P_w_cond_in-DELTAP_cond  { psia } 

P_w_cond_ave=(P_w_cond_in+P_w_cond_out)/2  (psia) 

(6.  Condenser  flow  rate} 

GPM_w_cond=3000  {gpni} 

CFM_w_cond=GPM_w_cond*convert(gpm,cfm)  { cfm } 

m_dot_w_cond=CFM_w_cond*rho_w_cond*convert(  1/min,  1/hr)  { Ibm/hr } 

(Section  VIL  Condenser  pump  model} 

(A.  Inputs  for  condenser  pump} 

h_pump_cond_in=h_w„to  w_out  { btu/lbm } 

P_pump_cond Jn=P_atm-Head_to  wtopump*rho_w_pump_cond*g*convert(lbm/ft*s'^2,psia)  { psia } 

Head_to  wtopump=Head_pipe_cond/2-Ele  v_tow  { ft } 

m_dot_w_pump_cond=m_dot_w_cond  {Ibm/hr} 

(B.  Condenser  pump} 
eta_pump_cond=.  83 
eta_motor_cond=.90 

eta_pump_cond=(h_pump_cond_out_rev-h_pump_condJn)/(h_pump_cond_out-h_pump_condJn) 
W_dot_pump_cond=m_dot_w_pump_cond*(h_pump_cond_out-h_pump_cond_in)/eta_motor_cond 
KW_pump_cond=W_dot_pump_cond*convert(Btu/hj*,kW)  (kW) 


(C.  Pressure  calculation} 

DELTAP__pump_cond=rho_w_pump_cond*g*Head_loop_cond*convert(lbm/ft*S'^2,psia)  { psia } 

P_pump_cond_out=P_pump_cond_in+DELTAP_pump_cond  (psia) 

P_pump_cond_ave=P_pump_cond_in+DELTAP_pump_cond/2  { psia } 

Head_pipe_cond=f_loop_cond*Length_loop_cond/D_i_pipe_cond*(FPS_pump_cond'^2/2)/ g  { ft } 


(eqn  6.30ref[5]} 

Head Joop_cond=Head_pipe_cond+Head_cond+Head_strainer  { ft } 

Head_strainer=30  (ft) 

f_loop_cond=:l/(-1.8*Iogl0(6.9/Re_D_i_piimp_cond+(epsilon_pipe_cond/(DJ_pipe_cond*3.7))''l.ll))'^2 

(eqn  6.64a  ref[5] } 

(D.  Pipe  geometry} 

D_Lpipe_cond=12*convert(in,ft)  (ft) 

Length Joop_cond=500  {ft} 

epsilon_pipe_cond=.00015  (ft) 

Elev_tow=:200  (ft) 
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{E.  Flowrates} 

GPM_pump_cond=m_dot__w_pump_cond/rho_w_pump_cond*convert(ft'^3/hr,gpm)  { gpm } 


CFM_pump_cond=GPM_pump_cond*convert(gpm,cfin)  { ft''3/ min } 

FPM_pump_cond=CFM_pump_cond/(pi*D_i_pipe_cond'^2/4)  { ft/min } 

FPS_pump_cond=FPM_pump_cond*convert(  1/min,  1/s)  { ft/s } 

FPH_pump_cond=FPM_pump_cond*convert(l/min,l/hr)  {ft/hr} 


Re_D_Lpump_cond=rho_w_pump_cond*FPH_pump_cond*DJ_pipe_cond/mu_w_pump_cond 
{F.  Fluid  Froperties} 

rho_W_pump_cond=density(WATER,T=T_pump_cond_ave,P=P_pump_cond_ave)  { lbm/ft'^3 } 
mu_W_pump_cond= viscosity  (W  ATER,T=T_pump_cond_ave,P=P_4)ump_cond_ave)  { Ibm/ft  hr } 

{G.  Temperatures) 

T_pump_cond_in=temperature(WATER,h=h_pump_cond_in,P=P_pump_cond_ave)  {F} 
T_pump_cond_ave=(T_pump_cond_out-T_pump_condJn)/ln((T_pump_cond_out)/(T_pump_cond_in)) 
T^ump_cond_out_rev=temperature(W  AT]^,P=P_pump_cond_out,s=s_pump_cond_in)  { F } 
T_pump_cond_out=temperature(WATER J^=P_pump_cond_out,H=h_pump_cond_out)  { F } 

{H.  Enthalpy  and  entropy} 

h_pump_cond_out_rev=enthalpy(W  ATER  J^=P_pump_cond_out,T=T_pump_cond_out_rev)  { btu/lbm } 
s_pump_cond_in=entropy(WATERJ^=P_pump_condJn,T=T_pump_cond_in)  {btu/lb  R} 

s_pump_cond_jout=entropy(WATERJP=P_pump_cond_out,T=T_pump_cond_out)  {btu/lb  R} 

{Section  VIII.  Cooling  tower  model} 

(A.  Inputs  for  cooling  tower} 

h_w_to  w_in=h_w__cond_out  { btu/lbm } 

m_dot_w_to  w=m__dot_w_cond  { Ibm/hr } 

{B.  Cooling  tower  calculations} 

Call  to werfan(Q_dot_to w,epsilon_tow,m_dot_a_to w,h_s_w_to w_in,h_a_to w_in,T_w_tow_out  1 : 
fanduty,T_w_tow_out) 

Q_dot_tow=0.9*fanduty*epsilon_tow*m_dot_a_tow*(h_s_w_towjn- 
h_a_to  w_in)+0. 1  *epsilon_to  w*m_dot_a_to  w*(h_s_w_to  w_in-h_a_to  w_in)  { B  tu/hr } 

Q_dot_tow=m_dot_w_to  w*(h_w_tow Jn-h_w_to  w_out)  { btu/hr } 

epsilon_tow=(  1  -exp(-NTU_tow*(  l-Cr_tow)))/(  l-Cr_tow*exp(-NTU„tow*(  l-Cr_tow))) 

{eqnl  1.30a  ref[3]} 

NTU_tow=1.3936*(m_dot_w_tow/m_dot_a_tow)'^(l“.7044)  jeqn  11.17  ref[9]} 

Cr_tow=(m_dot_a_to  w*Cs)/(m_dot_w_to  w*Cp_w__to  w)  { eqn  n.  1 2  ref[9] } 

Cs=(h_s_w_tow_in-h_s_w_tow_out)/(T_w_towJn-T_w_tow_out){btu/lbm  R}  {eqn  11.13  ref[9] } 

{C.  Water  properties} 

rho_w_tow=density(WATER,T=T_w__tow_aveJ^=P_atm) 
Cp__w_tow=specheat(WATER,T=T_w_tow_aveJ’=P_atm) 
h_s_w_tow_in=enthalpy(  AIRH20,R=  1  ,T=:T_w_to  w_in,P=P_atm) 
h_s_w_tow_out=enthalpy  ( AIRH20  Jl=  1  ,T=T_w_tow_out  J^=P_atm) 

{D.  Air  Properties} 

rho_a_tow=density(  AIRH20,B=ODWB  ,T=T_a_tow_in  J*=P_atm)  { lbm/ft'^3 } 

Cp_a_tow=specheat(  AIRH20,B=ODWB  ,T=T_a_to  w Jn,P=P_atm)  { Btu/lbm  R } 

h_a_tow_in=enthalpy(AIRH20,B=ODWB,T=T_a_tow_in,P=P_atm)  {Btu/lbm} 


{lbm/ft'^3} 
{Btu/lbm  R} 
{Btu/lbm} 
{Btu/lbm} 


{E.  Temperatures} 

T_a_tow_in=T_amb  {F} 

T_w_to  w Jn=temperature(W  ATER,h=h_w_tow Jn  JP=P_atm)  { F } 

T_w_tow_out  l=temperature(WATER,h=h_w_tow_out,P=P_atm)  {F} 

T_w_to w_ave=(T_w_to wJn-T_w_to w_out)/ln((T_w_to w Jn)/(T_w_t03v_out))  { F } 

{F.  Entropy} 

s_w_tow_in=entropy(WATER,h=h_w_tow„inJ^=P_atm)  {btu/lbm  R} 

s_w_tow_out=entropy(WATER,h=h_w_tow_outJP=P_atm)  {btu/lbm  R} 

(G.  Flowrates} 

ni_dot_a_to  w=rho_a_to  w*cfm_to  w*con  vert(  1/min,  1/hr)  { Ibm/hr } 

cfm_tow=250*  1000 

{H.  Cooling  tower  fan  model} 

FHPCA=-2.3084  (Appendix  A  ref[10]} 

FHPCB=6.3769  (Appendix  A  ref[10]} 

W_dot_fan_tow=(exp(FHPCA+((ln(cfm_tow)«10.5)/2)*FHPCB)/.9)*fanduty*convert(hp,btii/hr) 

(btu/hr }  (Appendix  A  ref[10] } 

KW_fan_tow=:W_dot_fan_tow*convert(btu/hr,kW)  {kW} 

(Section  IX.  Constants} 

P_atm=14.7 
g=32.2 
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{psia) 

{ft/s''2} 
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APPENDIX  B 

BASE  CASE,  PRIMARY-SECONDARY  SYSTEM 


"Base  Case,  Primary-Secondary,  Chilled  Water  System" 

{Section  I.  Procedures} 

(A.  Air  handler  procedures} 

(This  PROCEDURE  determines  the  load  based  on  out  door  dry  bulb  temperature} 

Procedure  Ioadprofile(T_amb:Q_ton_ah) 

IF  (T_amb<=60)  Then 

Q_ton_ah=20 
Endif 

IF  (T_amb>60)  AND  (T_amb<=105)  Then 

Q_ton_ah=16/7*T_amb-117.143 
Endif 

IF  (T_amb>105)  Then 

Q_ton_ah=100 
Endif 
END 

( This  procedure  determines  the  flow  regime  and  assigns  a  Nusselt  number  and  friction  factor} 
Procedure  regime(Re_D_i_ahJ*r_w_ah,D_i_ah,epsilon_pipe_ah:Nus_ah,f_i_ah) 

IF  (Re_D_i_ah<=2000)  Then 

Nus_ah=4.36  {eqn  3.65  refll]} 

f_i_ah=64/Re_D_i_ah  {eqn3.29  ref[l]} 

Endif 

IF  (Re_DJ_ah>2000)  Then 

Nus_ah=0.023*Re_D_i_ah''0.8*Pr_w_ah'H).4  (eqn  11.16ref[2]} 

f_i_ah=l/(-1.8*loglO(6.9/Re_D_i_ah+(epsilon_pipe_ah/(D_i_ah*3.7))^l.l  l))''2{eqn  6.64a  ref[5] } 
Endif 
END 

{This  PROCEDURE  determines  the  relative  heat  capacities} 

Procedure  heat(C_a_ah,C_w_ah:C_min_ah,C_max_ah) 

IF  (C_a_ah/C_w_ah<=l)  THEN 
C_min_ah=C_a_ah 
C_max_ah=C_w_ah 
Endif 

IF  (C_a_ah/C_w_ah>I)  THEN 
C_min_ah=C_w_ah 
C_max_ah=C_a_ah 
ENDif 
End 
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{This  procedure  determines  the  heat  exchanger  efficiency  given  that  the  air  flow  is 
unmixed  and  the  water  is  mixed } 

Procedure  transfer(C_a_ah,C_w_ahJ^U_ah,C_r_ah:epsilon_ah) 

IF  (C_a_ah/C_w_ah>l)  THEN 

epsilon_ah=l-exp((-l/C_r_ah)*(l-exp(-C_r_ah*NTU_ah)))  (eqn  1 1.35a  ref[3] } 

Endif 

IF  (C_a_ah/C_w„ah<=l)  THEN 

epsilon_ah=(l/C_r_ah)*(l-exp(-C„r_ah*(l-exp(-NTU,ah))))  (eqn  1 1.34a  ref[3] } 

Endif 

END 

(This  procedure  determines  the  f  for  air  flow  over  bare  tubes. } 

Procedure  baretube(Re_D_o_ah:f_tube_ah) 

IF  (Re_D_o_ah<1000)  THEN 

f_tube^ah=3.680886*Re_.D_o_ah^(-0.3090731)  (Figure  7.13  ref[l] } 

Endif 

IF  (Re_D^o_ah>=1000)  THEN 

Uube_ah=0.42275-5.53e-6*Re_.D^o^ah-3.598e-10*Re^D^o^ah''2{Figure  7.13  ref[l]} 
Endif 
END 

(This  procedure  determines  the  effective  specific  heat  for  air  of  the  air  handler } 

Procedure  latentheat(T_a_ahJn,T_a_ah_out,T_a_ah_.ave,D_^a_ahJn,HMRT^a_ahJnJ»_atm, 
P_a_ah_in,h_a__ah_in:HMRT_a_ah_out,h_a_ah_out,h_a_ah_out_eff,Cp__eff) 

IF  (T_a_ah_out>D_a_ahJn)  Then 
HMRT_a_ah_out=HMRT_a_ahJn 

h_a_ah_out=enthalpy(AIRH20,T=T_a_ah_outJP=P__atm,W=HMRT_a_ah_out) 

h_a_ah_out_eff=h_a_ah__out 

Cp_eff=specheat(AIRH20J=T_a_.ah_aveJP=(P^a_ahJn+P_atm)/2,D=D_a_ahJn) 

Endif 

IF  (T_a_ah_out=<D_a_ah_in)  Then 

HMRT_a__ah_out=humrat(  AIRH20,T=T_a__ah_out  JP=:P_atm,R=  1 ) 
h_a_ah_out=enthalpy(  AIRH20,T=T_a_ah_out  J^=P_atm,R=  1) 
h_a_ah_out_eff=h_a_ah_out+(HMRT_a_ah_in- 
‘HMRT^a_ah_out)*enthalpy(WATERJ>=P^atm,T=T_a_ah^out) 

Cp_eff=(h_a_ah_in*'h_a_ah_out_eff)/(T_a_ahJn-T_a__ah_out) 

Endif 

END 
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{B.  Chilled  water  reset  procedure} 

{This  procedure  first  sets  mass  flow  rate  as  a  function  of  temperature,  then  at  a  specified  outdoor  dry  bulb, 
it  switches  to  make  the  evaporator  exit  temperature  a  function  of  flow  rate.) 

Procedure 

reset(T_amb,QJoad_ah,h_w_ah_in,h_w_ah_out,h_pump_sec_out,h_pump_sec_in,m__dot_w_ah,P_w_eva 
p_out:  m_dot_  w_pump_sec,T_w_evap_out) 

IF  (T^amb=>77)  Then 
T_w_evap_out=44 

m_dot_w_pump„sec=  1 0*Q Joad_ah/(h„w_ah_out'h_  w^ahjn) 

Endif 

IF  (T_amb<77)  Then 

m_dot_w_pump_sec=538300 

T_w„evap_out=temperature(WATER,h=(h_w_ah_out-h_pump_sec_out+ 

h_pump_sec_in-Q_load_ah/m_dot_w_ah),P=P_w_evap_out) 

Endif 

END 

{ C.  Tower  procedures. ) 

{ This  procedure  determines  the  tower  fan  duty } 

Procedure 

towerfan(Q_doMow,epsiIon_tow,m_dot_a_tow,h_s_w_tow_in,h_a_towJn,T_w_tow_outl: 

fanduty,T_w_tow_out) 

IF  (T_w_tow_outl>70)  Then 

fanduty=l 

T_w_to  w_out=T_w_to  w_out  1 
Endif 

IF  (T_w_tow_out  1<70)  AND  ((Q_dot_to  w-0. 1  *epsilon_tow*m_dot_a_to  w*(h_s_w_to  w_in- 

h_a_tow_in))/(0.9*epsilon_tow*m„dot__a_,tow*(h_s_w_tow_in-h__a_towjn))>0)  Then 

fanduty=(Q_doMow-0.1*epsilon_tow*m_dot__a_tow*(h_s_w_towJn-h_a_tow_in))/ 
(0.9*epsilon_tow*m_dot_a_tow*(h_s_w_tow_in-h„a_tow_in)) 

T_w_tow_out=70 

Endif 

IF  (T_w_tow_outl<70)  AND  ((Q_dot_tow-0. 1  *epsilon_tow*m_dot_a_tow*(h_s_w_tow_in- 

h_a_tow_in))/(0.9*epsilon_tow*m_dot_a_tow*(h_s_w__tow_in-h_a_towjn))<=0)  Then 
fanduty=0 

T_w_to  w_out=T_w_to  w_out  1 
Endif 
END 


{D.  Mixing  procedures.} 

Procedure  mixing 

(m_dot_w_bypass,m_dot_w_pump_pri,m_dot_w_pump_sec,h_w_ah_out,h_w__evap_out: 

h__w_bypass,h_pump_prijn,h_pump_sec_in) 

IF  (m_dot_w_bypass>=0)  Then 
h_w_bypass=h_w_evap_out 

h_pump_pri_in=(m_dot_w_bypass*h_w_evap_out+m_dot_w_pump_sec*h_w_ah_out)/ 

m_dot_w_pump_pri 

h_pump_sec_in=h_w_evap_out 

Endif 

IF  (m_dot_w_bypass<0)  Then 
h_w_bypass=:h_w_ah_out 
h_pump_pri_in=h_w_ah_out 

h_pump_secJn=(m_dot_w_pump_pri*h_w_evap_out-m_dot_w_bypass*h_w_ah_out)/ 

m_dot_w_pump_sec 

Endif 

END 

{ Section  11.  Inputs } 

Hours=l 

T_amb=95  {F} 

ODWB=  78  {F} 

{Section  III.  Overall  Calculations} 

{ A.  Energy  analysis } 

W_dot_total=W_dot_comp+W_dot_pump_tot+W__dot_pump_cond+10*W_dot_fan+W_dot_fan_tow 
KW_total=W_dot__total*convert(btu/hr,kW)  { kW } 

KWh_total=KW^total*hours  {kWh} 

{These  equations  sum  to  the  inverse  of  the  lER  value} 
IER_chilIer=((Hours/Hours_total)/(KW_comp/Q_ton__evap)) 

IER_sy  stem=(Hours/Hours_total)/ (kW_total/(  1 0*Q_ton_ah)) 

IER_aux_evap=(Hours/Hours_total)/((W_dot_pump_tot+ 1 0*  W_dot_fan)  *con  vert(btu/hr  ,kW)/ 
(iO*Q_ton_ah)) 

Hours_totaI=7 1 32 
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{B.  Entropy  Analysis} 

S_dot_gen_fan_ah=  10*m_dot_a_ah*(s_a_fan_out-s_a_fan  Jn)  { btu/hr  R } 

S_dot_gen_ah=  10*(ni_dot_w_ah*(s_w_ah_out-s_w_ah_in)-Q  Joad_ah/(T_amb+459,67))  { btu/hr  R } 
S_dot_gen_pump_pri=:m_dot_w_pump_pri*(s_pump_pri_out-s_pump_pri_in)  { btu/hr  R } 

S_dot_gen_pump_sec=m_dot_w_pump_sec  *(  s_pump_sec_out-s_pump_sec  Jn)  { btu/hr  R } 

S_dot_gen_mixing=m_dot_w_pump_pri*s_pump_priJn- 

(m_dot_w_bypass*s_w_bypass+ni_dot_w_puinp_sec*s_w_sec)  { btu/hr  R } 

S_dot_gen_evap=m_dot_w_evap*(s_w_evap_out-s_w_evapJn)-(-Q_dot_evap/(T_evap+459.67)) 
S_dot_gen_cond=m_dot_w_cond*(s_w_cond_out-s__w_cond_in)  -(Q_dot_cond/  (T_cond+459 .67)) 
S_dot__gen_comp=(W_dot_comp“  W_dot_camot)/T__anib  { btu/hr  R } 

S_dot_gen_tower=m_dot_w_tow*(s_w_tow_out-s_w__tow_in)-(-Q_dot__tow/(T_w_tow_ave+459.67)) 
S_dot_^en_pump_cond=m_dot_w_cond*(s_pump_cond_out-s_pump_cond_in)  { btu/hr  R } 

S_dot_genJoop_pri=m_dot_w_pump_pri*(s_pump_sec_in-s_w_evap_out)  { btu/hr  R } 

S_dot_genJoop_sec=m__dot_w_pump_sec*(s_w_sec-s_w__ah_out+s_w_ahJn-s_pump_sec_out) 
S_dot_gen_loop_cond=m_dot_w_cond*(s_pump__cond_in-s_w_tow_out+s_w_tow_in-s_w__cond_out) 
S_dot_gen__tot=S_dot_gen_fan_ah+S_dot_gen_ah+S_dot_gen_pump_pri+S_dot_gen_pump_pri+S_dot_ge 
n_mixing+S_dot_gen_evap+S_dot_gen_cond+S_dot_gen_comp+S_dot_gen_tower+S_dot_gen_pump_co 
nd+S_dot_gen_Ioop_pri+S_dot_genJoop_sec+S_dot_gen_loop_cond 

W_dot Jost=S_dot_gen_tot*T_amb  { btu/hr } 

{Section  IV.  Air  handler  model } 

(A.  Air  handler  fan} 

(1.  Fluid  properties} 

T_room=75  {F} 

T_a_fan Jn=temperature(  AIRH20,h=h_a_fan Jn,P=P_atm,W=HMRT„fan Jn)  { F } 

T_a_fan_out=temperature(AIRH20,h=h_a„fan_out,P=P_atm+DELTAP_a_ah,W=HMRT_fan_in) 
HMRT^amb=humrat(AIRH20,T=T_amb,P=P_atm,B=ODWB) 
HMRT_room=humrat(AIRF[20,T=T_room,P=P__atm,R=:R_room) 

HMRT3n_in=0.85*HMRT_room+0. 15  *HMRT_amb 
R_room=.5549 

h_amb=enthalpy(AIRH20,T=:T_amb,P=P__atm,B=0DWB)  (btu/lbm } 

h_room=enthalpy(  AIRH20,T =T_room+5  J*=P_atm,R=R_room)  { btu/lbm } 

h_a_fan_in=0.85*h_room+0.15*h_amb  (btu/lbm) 

h_a_fan_out_rev=enthalpy(AIRH20,P=P_atm+DELTAP__a_ah,T=T_a_fan_out_rev,W=HMRT_fan_in) 
s_a_fan_in=entropy  ( AIRH20,T =T_a_fan_in,P=P_atm,  W=HMRT_fan_in)  { btu/lbm  R } 

s_a_fan_in=entropy(AIRH20,T=T_a_fan_out_rev,P=P_atm+DELTAP_a_ah,W=HMRT_fan_in) 
s_a_fan_out=entropy(AIRH20,T=T_a_fan_out,P=P_atm+DELTAP_a_ah,W=HMRT_fanJn) 

(2.  Fan  Calculations) 
eta_fan=.65 

eta_fan=(h_a_fan_in-h_a_fan_out_rev)/(h_a_fan_in-h_a_fan_out) 


W_dot_fan=-m_dot_a_ah*(h_a_fanJn-h_a_fan_out)  {btu/hr) 

KW_fan=W_dot_fan*convert(btii/hr,kW)  jkW) 

{B.  Inputsfor  Air  Handler) 

h_w_ah Jn=h_pump_sec_out  { btu/lbm } 

h_a_fan_out=h_a_ah_in  { btu/lbm } 

P_w_ah_in=P„pump_sec_out- 

(Head_pipe_sec/2+EIev_load+Head_valve)*rho^w_ah*g*convert(lbm/ft*s'^2,psia) 
P_a_ah_in=:P_atm+DELTAP_a_ah  (psia) 

Head_valve=20-Head_coil  {ft} 

HMRT_a_ahJn=HMRT_fan_in 
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{C.  Air  Handler  overall  heat  transfer  coefficient  calculation} 

lAJA_tot_ah=l/(hJ__ah*AJ_ah)+ln(D_o_ah/D_i_ah)/(2*pi*k_pipe*Length_ah*n_row*n_column)+l/(h_ 
o^ah* A„o_ah*eta_o_ah)  { eqn  1 1 .5  ref[ 3] } 


{D.  Geometry} 

Length_ah=7 .555  { ft } 

Height_ah=5  {ft} 

Depth_ah=8.5*convert(in,ft)  {ft} 

fin_pitch=12/convert(in,ft)  { 1/ft} 

D_o_ah=.75  *convert(in,ft)  { ft } 

E>_Lah=.666*convert(in,ft)  { ft } 

t_fm=.006*convert(in,ft)  { ft } 

n_row=4 
n_column=40 

Elevjoad=0  {ft} 


{E.  Airside  calculations} 

{1.  Air  properties} 

Call  latentheat(T_a_ah_in,T_a_ah_out,T_a_ah_ave,D_a_ah_in,HMRT_a_ah_inJ*_atm, 

P_a_ah_in,h__a_ah_in:HMRT_a_ah_out,h_a_ah_out,h_a_ah_out_eff,Cp_eff) 

HMRT_a_ah_ave=(HMRT_a_ah_in+HMRT_a_ah_out)/2 

k_a_ah=conductivity(  AIRH20,T=T_a_ah_ave,P=P_a_ah_ave,W=HMRT_a_ah_ave)  { Btu/hr-ft-R } 

mu_a_ah= viscosity  (AJRH20,T =T_a_ah_ave,P=P_a_ah_ave,W=HMRT_a_ah_ave)  { Ibm/ft-hr } 

rho^a^ah=density(  AIRH20,T=T__a_ah_ave,P=P^a^ah_ave,W=HMRT_a_ah__ave)  { lbm/ft''3 } 

v_m_ah=  l/rho_a_ah  { ft^3/lbm } 

Pr_a_ah=Cp_eff*mu_a_ah/k_a_ah 


{2.  Air  flow  properties} 
G_a_ah=m_dot_a_ah/(sigma_ah*A_fr_ah) 
Re_D_o_ah=G_a_ah*D_o_ah/mu_a_ah 
Re_L_ah=G_a_ah*b_ah/mu_a_ah 


{ lbm/ft''2-hr }  { eqn  1 1 .37  ref[3] } 

{eqn  14.37  ref[4]} 
{eqn  ref[4] } 


{3.  Airside  geometry  calculations} 
a__ah=Height_ah/n_column 
b__ah=Depth_ah/n_row 
A_fr_ah=Height_ah*Length_ah 
sigma_ah=(a_ah-D_o_ah)/a_ah 

A_f_ah=(2*a_ah*b_ah-pi*D_o_ah'^2/2)*fin_pitch*Length_ah*n_row*n_column 

A_t_ah=pi*D_o_ah*Length_ah*n_row*n_column 

A_o_ah=A__f_ah+A_t_ah 


{ft} 

{ft} 

{ft'^2} 

{fr2} 

{ft^2} 

{ft^2} 


{ 4.  Airside  Heat  Transfer  Calculations } 
h_o_ah=(Cp^eff*G_a^ah*j_n/Pr_a_ah''(2/3))  {Btu/ft^2-hr-R} 

j_4=0.2618*JP+0.0014 

JP=Re^D^o_.ah^(-0.4)*(A^o^ah/A_,Lah)''(-0.15) 

j^n/j_4=(M280*n_row*Re_L^ah''(-1.2))/(M280*4*Re^L_ah'^(-1.2)) 


{eqn  14.36  ref[4]} 
{eqnref[4]} 

{eqn  14.46  ref[4]} 
{eqn  14.50  ref[4]} 


{5.  Fin  Calculations} 

eta_o_ah=l-(A_f_ah/A_o_ah)*(I-eta_fin)  (eqn  1 1.3  ref[3] } 

eta_fin=.95 
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{ F.  Waterside  calculations } 

{1.  Water  properties} 

Cp_w_ah=specheat(WATER,T=T_w_ah_ave,P=P_atm) 

k_w_ah=conductivity(WATER,T=T_w_ah_ave,P=P_atm) 

mu_w_ah=viscosity(WATER,T=T_w_ah_ave,P=P_atm) 

rho_w_ah=density(WATER,T=T_w_ah_ave,P=P_atm) 

Pr_w_ah=Cp„w_ah*mu_w_ah/k_w__ah 

{2.  Flow  properties} 

G_w_ah=m_dot_w_ah/((n_column-l)*pi*DJ_ah'^2/4) 

Re_D_i_ah=G_w_ah*D_i_ah/mu_w_ah 

{3.  Waterside  geometry  calculations} 
A_Lah=pi*D_Lah*Length_ah*n__row*n_column 


{Btu/Ibm-R} 

{Btu/hr-ft-R} 

{Ibm/ft-hr} 

{lbm/ft'^3} 


{Ibm/ft^2-hr} 

{eqn  11.10ref[2]} 


{ft'^2} 


{4.  Waterside  heat  transfer  calculations} 

h_Lah=Nus_ah*k_w_ah/D_i_ah  {Btu/ft'^2-hr-R}  {eqn  3,66  ref[l]} 

Call  Regime(Re_D_Lah,Pr_w_ah,D_Lah,epsilon_pipe_ah:Nus_ah,f_Lah) 


(G.  Pipe  and  fin  properties} 
k_pipe=40 1  *convert(W/m*K,Btu/hr  *ft*R) 
k^fin=237*convert(W/m*K,Btu/hr*ft*R) 
epsilon_pipe„ah=.000 1 5 


(Btu/hr-ft-R} 

{Btu^r-ft-R} 

{ft} 


{H.  Load  conditions.} 

Call  loadprofile(T_amb:Q_ton_ah) 

CALL  heat(C_a_ah,C_w_ah:C_min_ah,C_max_ah) 

CALL  Transfer(C_a_ah,C_w_ah,NTU„ah,C_r_ah:epsilon_ah) 

Q_load_ah=Q_ton_ah*convert(ton,Btuyhr) 

"Q_load_ah=m_dot_w_ah*(h_w_ah_out-h_w_ahJn) 

Q_load_ah=m_dot_a_ah*Cp_eff*DELTAT_a_ah 

Q_max_ah=C_min_ah*(T_a_ahJn-T_w_ahJn) 

C_w_ah=m_dot_w_ah*cp_w_ah 

C_a__ah=m_dot_a_ah*Cp_eff 

C_r_ah=C_min_ah/C_max_ah 

Q_load_ah=epsilon_ah*C_min_ah*(T_a_ah_in-T_w_ah_in) 

NTU_ah=UA_tot_ah/C_min_ah 


{Btu/hr} 

{Btu/hr}  {Blocked  out  for  reset}” 
{Btu/hr} 

{Btu/hr} 

{Btu/hr  R} 

{Btu/hr  R} 

{Btu/hr} 

•{eqn  11.25  ref[3]} 


{I.  Temperature} 

T_w_ah Jn=temperature(WATER,H=h_w„ah_in,P=P_w_ah  Jn)  { F } 

T_w_ah__out=temperature(WATER,H=h_w_ah_out,P=P_w_ah_out)  { F } 

T_w__ah_ave=(T_w_ah__in-T_w_ah_out)/ln(T_w_ah_in/T_w_ah_out)  {F} 

T_a_ah_a  ve=(T_a_ah_in-T_a_ah_out)/ln(T_a_ah  Jn/T_a_ah_out)  { F } 

T_a_ah_in=temperature(  AIRH20,h=h_a_ah Jn,P=P_a_ah_in,W=HMRT_a_ah_in)  { F } 

T_a_ah_out=58  { P } 

D_a_ah__in=dewpoint(AJRH20,T=T__a_ah_in,P=P_a_ah_in,W=HMRT_a_ah_in)  {F} 

DELTAT_w_ah=T_w_ah_out-T_w_ahJn  {F} 

DELTAT_a_ah=T_a_ah_in-T_a_ah__out  {F} 
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{J.  Entropy} 

s_w_ahJn=entropy(WATER,P=P_w_ahJn,T=T_w_ah_in)  {btu/lbm  R} 

s_w_ah_out=entropy(W  ATER,P=P_w_ah_out,T=T_w_ah_out)  { btu/lbm  R } 


{K.  Flow  rates} 

{1.  Water  flow} 

V_w_ah=FPS_w_ah*(n_column- 1  )*pi*D_Lah^2/4*convert(  1  /s,  1/min)  { cfm } 


GPM_w_ah= V_w_ah  *con  vert(ft^3  ,gal)  { gpni } 

m_dot_w_ah=V_w_ah*rho_w_ah*60  {Ibm/hr} 

{2.  Airflow} 

V__a_ah=FPM_a_ah*Length_ah*Height_ah  { cfm } 

m_dot_a_ah= V_a_ah*rho_a_ah*60  { IbmAir } 

FPH_a_ah=FPM_a_ah*convert(  1/min,  l/hr)  { ft/hr } 

FPS_a_ah=FPM_a_ah*convert(  1/min,  1/s)  { ft/s } 

{L.  Pressure  drop  calculations } 

{1.  Waterside  pressure  drop} 

DELTAP_Lah=rho_w_ah*g*Head_coil*convert(lbm/ft*s'^2,psia)  { psia } 

DELTAP_tube_ah=LLah*Length_ah*n_row/DJ_ah*(rho_w_ah*FPS_w_ah'^2/2) 
*convert(lbm/s'^2*ft,psia)  { psia } 

Head_tube_ah=DELTAP__tube_ah/(rho_w_ah*g)*convert(psia*ft'^2*s^2/lbm,ft)  {ft} 
Head_bend=5*K*FPS_w^ah''2/(2*g)  { ft } 

Head_coil=Head_tube_ah+Head_bend  { ft } 

P_w_ah_out=P_w__ah_in-DELTAP_Lah  {psia} 

P_w_ah_ave=(P_w_ah_in+P_w_ah_out)/2  { psia } 

K=.4 


{eqn3.27ref[l]} 
{eqn  6.25  ref[6]} 


{2.  Airside  pressure  drop} 

CALL  Baretube(Re_D_o_ah:f_tube_ah) 

DELTAP_a_ah=deltaP_o„ah*convert(lbni/hr^2*ft,psia) 

P_a_ah_ave=(P_a_ah_in+P_atm)/2 

f_f.ah=(deltaP_o_ah-deltaP_t_ah)/(v_m_ah*G_a_ah'^2/(2*sigma_ah)) 

f-f-ah=l  .7*Re^L^ah''(-0.5) 

a_ah_star=a_ah/D_o_ah 

b_ah_star=b__ah/D_o_ah 

ChLah=l 

deltaP_t_ah=n__row*f_tube_ah*ChLah*(rho_a_ah*(FPH__a_ah/sigma__ah)^2/2) 


{psia} 

{psia} 

{eqn  7  ref[6]} 
{eqn9ref[6]} 
{eqn7.114ref[l]} 
{eqn7.114ref[l]} 
{Fig  7.13  ref[l]} 
{lbm/hr^2  ft} 

{eqn  7.113  refll]} 


{Section  V.  Pump  models} 

{A.  Primary/Secondary  Mixing  Calxulations} 

CALL  mixing  (m_dot_w_bypass,m_dot_w_pump_pri,m_dot_w__pump_sec,h_w_ah__out,h_w_evap_out: 
h_w_bypass,h_pump_prLin,h_pump_secJn) 

m_dot_w_bypass=m_dot_w_pump__pri-m__dot_w_pump_sec  { Ibm/hr } 


{ B.  Pumping  power  calculations } 

KW_pump_tot=KW_pump_pri+KW_pump_sec 

W_dot_pump_tot=W_dot_pump_pri+W_dot_pump_sec 


{kW} 

{Btu/hr} 
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{C.  Primary  pump  model } 

{ I .  Inputs  for  primary  pump  model } 

s_w_bypass=entropy(WATER,P=P_min,h=h_w_bypass)  {btu/lbm  R} 

s_w_sec=entropy(  WATER, P=P_min,h=h_w_ah_out)  { btu/lbm  R } 

{2.  Pressure  calculations} 

DELT  AP_pump_pri=rho_w_pump_pri*g*Head  Joop_pri*convert(lbm/ft*S'^2,psia)  { psia } 

P_pump_pri_out=P_min+DELT  AP_pump_pri  { psia } 

P_pump_pri_ave=P_min+DELT  AP_pump_pri/2  { psia } 

P_min=i25  {psia} 


(3.  Pump  effeciencies} 
eta_pump_pri=.83 
eta_motor_pri=.90 

eta_pump_pri=(h_pump_pri_out_rev-h_pump_prLin)/(h_pump_pri„out-h_pump_pri„in) 
W_dot_pump_pri=m_dot„w_pump_pri*(h_pump__pri_out-h_pump_pri  Jn)/eta_motor_pri  { Btu/hr } 
KW_pump_pri=W_dot_pump„pri*convert(Btu/^,kW)  { kW } 

{4.  Head  calculation} 

Head_pipe_pri=Moop_pri*LengthJoop_pri/D_i_pipe_pri*(FPS_pump_pri^2/2)/g  {ft} 

Head_ioop_pri=Head_pipe_pri+Head_e  vap  { ft } 

f Joop_pri=  l/(- 1 .8*log  10(6.9/Re__D_Lpump„pri+(epsilon_pipe_pri/(D_Lpipe_pri*3.7))'^  1 . 1 1))'^2 

{5.  Pipe  geometry} 

D  J_pipe_pri=  1 0*con  vert(in,ft)  { ft } 

Length_loop_pri=500  { ft } 

epsilon_pipe_pri=,00015  {ft} 


{ 6.  Flow  properties } 

GPM_pump_pri=24CK)  { gpm } 

CFM_pump_pri=GPM_pump_pri*convert(gpm,cfm)  { ft^3/min } 

FPM_pump_pri=CFM_pump_pri/(pi*DJ_pipe_pri'^2/4) 

{ft/min} 

FPS_pump_pri=FPM_pump_pri*convert(l/min,l/s)  {ft/s} 

FPH_pump_pri=FPM_pump_pri*convert(  1/min,  1/hr)  { ft/hr } 

m_dot_w_pump_pri=rho_w_pump_pri*GPM_pump_pri*convert(gpm,cfm)*60  {Ibm/hr} 

Re_D_i_pump_pri=rho_w_pump_pri*FPH__pump_pri*DJ_pipe„pri/mu_w_pump_pri 

{7.  Ruid  properties} 

rho_w_pump_pri=density(WATER,T=T_pump_prLave  J*=P_pump_pri_ave)  { lbm/ft'^3 } 

mu_w_pump_pri=viscosity(WATER,T=T_pump_pri__ave,P=P_pump_pri_ave)  { Ibm/ft  hr } 

{8.  Temperatures} 

T_pump_pri_in=temperature(  W  ATER,h=h_pump_pri_in  J^=:P_min)  { F } 

T_pump_pri_ave=(T_pump_pri_out-T_pump_pri_in)/ln((T_pump_prLout)/(T_pump_pri  Jn))  { F } 
T_pump__pri_out_rev=temperature(W  ATER,P=P_pump_pri_out,s=s_pump_pri_in)  { F } 

T_pump_prLout=temperature(  W  ATER,P=P_pump_pri_out,H=h_pump_pri_out)  { F } 

(9.  Enthalpy  and  entropy} 

h_pump_pri_out_rev=enthalpy(WATER,P=P„pump_pri_out,T=T_pump_pri_out_rev)  { Btu/lbm } 

s_pump_pri_in=entropy(W  ATER,P=P_min,T=T_pump_pri  Jn)  {Btu/lbm  R } 

s_pump_pri_out=entropy(WATER,P=P_pump_prLout,T=T__pump_pri_out)  { Btu/lbm  R } 
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{D.  Secondary  pump  model } 

{ I .  Inputs  for  secondary  pump  model } 

P_pump_secJn=P_w__evap_out-(rho_w_pump_sec*g*Head_pipe_pri/2)*convert{lbm/s^2*ft,psia) 


m_dot_w_pump_sec=10*m_dot_w_ah  {Ibm/hr} 

{2.  Pressure  calculations} 

DELT  AP_pump_sec=rho_w__pump_sec*g*(Head_loop_sec)*convert(lbm/ft*S'^2,psia)  { psia } 
P_pump_sec_out=P_pump_sec_in+DELT  AP_pump__sec  { psia ) 

P„pump_sec_a  ve=P_pump_sec_in+DELT  AP_pump_sec/2  { psia } 


{3.  Affinity  Calculations} 

N_pump_act=N_pump_design*(GPM_pump_sec/GPM_design)  {rpni}  {eqn  6.1  ref[7]} 

N_pump_act=N_pump_design*(Head Joop_sec/Head_design)'^(0.5)  { rpm }  { eqn  6. 1  ref[7] } 
N_pump_design=  1180  { rpni }  { ref[8] } 

{4.  Pump  curves} 

Head_design=:7 1.39902  +  0.01 53946  l*GPM„design  -  0.00000846*GPM^design'^2 
eta„pump_sec=(0,036037 15  +  0.00073264*GPM_design  -  0.00000017*GPM^design^2)* 

(0.8158200  +  0.402540  l*N_.pump_act/N_pumpJesign  -  0.21 8360 l*(N„pump^act/N_,pump„design)'^2) 
eta_motor_sec=0.0498386+ 1 .846339*(N_pump_act/N_pump_design)- 
1 .029*(N_pump_act/N_pump_design)'^2 

eta_pump_sec=(h_pump_sec_out_rev-h_pump_sec_in)/(h__pump_sec_out-h_pump_sec_in) 
W_dot_pump_sec=m_dot_w_pump_sec*(h_pump_sec_out-h_pump_sec_in)/eta_motor_sec  { Btu/hr } 
KW_pump_sec=W_dot_pump_sec*convert(Btu/lv,kW)  { kW } 

{5.  Head  calculation} 

Head_pipe_sec=f_loop_sec*LengthJoop_sec/DJ_pipe_sec*(FPS_pump_sec^2/2)/g 

{ft}  (eqn  6.30  ref[5]} 

Head  Joop_sec=Head_pipe_sec+Head_coil+Head_val  ve  { ft } 

f Joop_sec=  l/(- 1 .8*log  10(6.9/Re_D J_pump_sec+(epsilon_pipe_sec/(D  J_pipe_sec*3 .7))'^  1 . 1 1))'^2 

(eqn  6.64a  ref[5]} 

(6.  Pipe  geometry} 

D J_pipe_sec=  i  0*con  vert(in,ft)  { ft } 

Length Joop__sec=  1500  {ft} 

epsilon_pipe_sec=.00015  {ft} 


{7.  Flow  properties} 

CFM_pump_sec=GPM_pump_sec*convert(gpm,cfm)  { ft''3/min } 

FPM_pump_sec=CFM_pump_sec/(pi*D_i_pipe_sec'^2/4)  { ft/  min } 

FPS_pump_sec=FPM_pump_sec*convert(l/min,l/s)  {ft/s} 

FPH__pump_sec=FPM_pump_sec*convert(  1/min,  1/hr)  { ft/hr } 

m_dot_w_pump__sec=rho_w_pump_sec*GPM_pump_sec*convert(gpm,cfm)*60  { Ibm/hr } 

Re_D_i_pump_sec=rho_w_pump_sec*FPH_pump_sec*D_i_pipe_sec/mu_w_pump_sec 

{8.  Fluid  properties} 

rho_w_pump_sec=density(W  ATER,T =T_pump_sec_a  ve,P=P_pump_sec_a  ve)  { lbm/ft^3 } 

mu_w_pump_sec= viscosity  (W  ATER,T =T__pump_sec_ave,P=P_pump_sec__ave)  { Ibm/ft  hr } 

{9.  Temperatures} 

T_pump_sec_in=temperature(WATER,h=h__pump_sec_in,P=P_pump_secJn)  {F} 

T_pump_sec_ave=(T_pump_sec_out-T_pump_sec__in)/ln((T_pump__sec_out)/ (T_pump_sec  Jn))  { F } 
T_pump_sec_out_rev=temperature(W  ATER,P=P_pump_sec_out,s=s_pump_sec_in)  { F } 

T_pump_sec_out=temperature(WATER,P=P_pump_sec_out,H=h_pump_sec_out)  {F} 
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( 10.  Enthalpy  and  entropy) 

h_pump_sec_out_rev=enthalpy(W ATER,P=P_punip_sec_out,T=T_pump_sec_out_fev)  { Btu/lbm } 
s_pump_secJn=entropy(WATER,P=P_pump_secJn,T=T_pump_secJn)  { Btu/lbm  R) 

s_pump_sec_out=entropy(WATER  J>=P_pump_sec_out,T=T_pump_sec_out)  { Btu/lbm  R } 


{Section  VI.  Centrifugal  chiller  model) 
(A.  Chiller  inputs) 


{L  Evaporator  Side} 

P_w_evap_in=P_pump__pri_out 

(psia) 

h_w_evap_in=h_pump_pri_out 

{btu/lbm} 

m_dot_w_evap=m__dot_w_pump_pri 

(Ibm/hr} 

{2,  Condenser  Side} 

P_w_cond„in=P_pump_cond_out 

(psia) 

h_w_cond_in=h_pump_cond_out 

(btu/lbm) 

{B.  Compressor) 

W_dot_camot=Q_dot_evap*((T_cond+459.67)/(T_evap+459.67)-l)  (btu/hr) 
KW_comp=W_dot_comp*convert(Btu/hr,kW)  { kW ) 

W_dot_comp=a+b*W_dot_camot+c*W_dot_camot''2 
a=3042.991  +  705.2175*Q_ton_evap  +  0.9521675*Q_ton_evap^2 
b=1.228996  -  0.00438602*Q_ton_evap  +  0.00000203 *Q_ton_evap^2 
c=0.00000279  -  3.I46E-09*Q_ton_evap  +  1.310E-12*Q_ton_evap''2 

(C.  Evaporator) 

{ 1 .  Evaporator  heat  transfer  calculations ) 

Call 

reset(T_amb,Q_load_ah,h_w_ah_in,h_w_ah_out,h_pump_sec_out,h_pump_sec_in,m. 
p_out:m_dot_w_pump_sec,T_w_evap_out) 

Q_dot_evap=Q_ton_evap*convert(ton, Btu/hr) 
Q_dot_evap=m_dot_w_evap*(h_w_evap_in-h_w_evap_out) 
Q_dot_evap=epsilon_evap*m_dot_w_evap*C_w_evap*(T_w_evapJn-T_evap) 

(eqn  11.19ref[3]) 

epsilon_evap=l-exp(-NTU_evap)  (eqn  11.30a  ref[3]) 

NTU_evap=UA_evap/(m_dot_w_evap*C_w_evap)  { eqn  1 1 .25  ref[3) } 

l/UA_evap=l/(hA_i_evap)+l/(hA_o_evap) 
hA_i_evap=6213*(Re_D_i_evap''0.5334881)*Pr_w_evap''0.3 
hA_o_evap=180935.0  +  6243.157*Q_ton_evap  -  1. 474305 ♦Q_ton_evap''2 

(2.  Evaporator  pressure  drop  calculations) 
f_Levap=l/(-1.8’'‘logl0(6.9/Re_D_i_evap+(epsilon_pipe_evap/(D_i_evap*3.7))''l.  1 1))''2 

(eqn  6.64a  ref[5) ) 

DELTAP_evap=Head_evap*rho_w_evap*g*convert(lbm/ft*s''2,psia)  { psia ) 

Head_evap=f_i_evap*Length_evap*n_pass_evap/DJ_evap*(FPS_evap''2/2)/g  { ft ) 

(eqn  6.30  refI5]) 
(psia) 

(psia) 


_dot_w_ah,P_w_eva 

(Btu/hr) 
(Btu/hr) 
(Btu/hr) 


P_evap_ave=P_min+DELTAP_evap/2 

P_w_evap_out=P_w_evapJn-DELTAP_evap 


{ 3 .  Evaporator  geometry } 


D_Levap=0.824*convert(in,ft) 

{ft} 

Length_evap=  14.2 

(ft) 

n_tubes_evap=2 15 

n_pass_evap=4 

epsilon_pipe_evap=.0003873 1 

{ft} 

{4.  Evaporator  fluid  properties} 

C_w„evap=specheat(WATER,T=T_w_evap_ave,P=P_evap_ave) 
rho_w_evap=density(WATER,T=T_w_evap_ave,P=P_evap_ave) 
mu_w_e  vap=viscosity  ( W  ATER,T =T_w_e  vap_ave,P=P_e  vap_a  ve) 
k_w_e  vap=conducti  vity  ( W  ATER,T =:T_w_evap_ave,P=P_evap_a  ve) 
Pr_w_evap=C_w__evap*mu_w_evap/k_w_evap 
Re_DJ_evap=rho_w_evap*FPH_evap*D_i_evap/mu_w_evap 


{Btu/lbmR} 

{lbni/ft^3} 

{Ibm/fthr} 

{Btu/fthrR} 


{5.  Evaporator  flow  rate} 

V__w_evap=GPM_w_evap*convert(gpm,cfm) 

{cfm} 

V_w_evap=(pi*D_i_evap'^2/4)*FPM__evap*n_tubes_evap 

{cfm} 

m_dot_w_evap=V_w_evap*rho_w_evap*60 

{Ibm/hr} 

FPS_evap=FPM_evap*convert(  1/min,  1/s) 

{ft/s} 

FPH_evap=FPM_evap*convert(  1/min,  1/hr) 

{ft/hr} 

{6.  Evaporator  Temperatures} 

’T_w_evap_out=44" 

{F}  {Blocked  out  for  reset} 

T_w_evap_in=temperature(W  ATER,h=h_w__evap_in,P=P_w_evap_in)  { F } 

T_w_evap_ave=(T_w_evap_in-T_w_evap_out)/In(T_w_evap_in/T_w_evap_out)  { F } 


{7.  Enthalpy  and  entropy} 

h_w_evap__out=enthalpy(WATER,T=T„w_evap_out»P=P_w__evap_out) 

s_w_evap_in=entropy(WATER,h=h_w_evap__in,P=P_w_evap_in) 

s_w__evap_out=entropy(WATER,h=h_w_evap_out,P=P_w__evap_out) 


{btu/lbm} 
{btuflbmR} 
{btu/lbm  R} 


(D.  Condenser} 

{ 1 .  Condenser  UA  calculation } 

Q_dot_cond=epsilon_cond*m_dot_w_cond*C_w_cond*(T_cond-T_w_cond_in)  (Btu/hr } 

(eqn  1.23  ref[3]} 

Q_dot_cond=m_dot_w_cond*(h_w_cond_out-h_w__cond_in)  {Btu/hr} 

Q_dot_cond=Q_dot_evap+W_dot_comp  { Btu/hr } 

Q_ton_cond=Q_dot_cond*convert(btu/hr,ton)  ( tons } 

epsilon_cond=l-exp(-NTU_cond)  ( eqn  1 .30  ref[3] } 

NTU„cond=UA_cond/(m_dot_w_cond*C_w_cond) 

UA_cond^l=1845702  +  22792.36*T^cond  -  96.24830*T_cond'^2 
UA_cond_2=L059937  -  0.00009226*Q_ton__cond  +  0.00000001 *Q_ton__cond'^2 
UA_cond=UA_cond_l  *UA_cond_2 


(2.  Condenser  fluid  properties } 

rho_w_cond=density(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 

C_w_cond=specheat(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 

mu_w_cond=viscosity(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 

k_w_cond=conductivity(WATER,T=T_w_cond_ave,P=P_w_cond_ave) 


{lbm/ft^3} 
{Btu/lbmR} 
{Ibm/fthr} 
{Btu/hr  ft  R} 
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{3.  Condenser  Temperatures} 

T_w_cond_ave=(T_w_cond_out-T_w_condJn)/ln(T_w_cond_out/T_w_condJn)  {F} 


T_w_cond_in=temperature(WATER,h=h_w_condJn,P=P_w_condJn)  {F} 

{4.  Enthalpy  and  entropy} 

h_w_cond_out=enthalpy(WATER,T=T_w_cond_out,P=P_w_cond_out)  { btu/lbm} 

s_w_cond  Jn=entropy  ( W  ATER,H=h_w_cond_in,P=P_w_cond_in)  { btu/lbm  R } 

s_w_cond_out=entropy(WATER,T=T_w_cond_out,P=P_w_cond_out)  { btu/lbm  R } 

{ 5 .  Condenser  pressures } 

Head_cond=60.63713*T_w„cond_ave^(-0.1908161) 

DELTAP_cond=rho_w_cond*g*Head_cond*convert(lbm/s'^2*ft,psia)  { psia } 

P_w_cond_out=P_w_cond_in-DELTAP_cond  { psia } 

P_w_cond_ave=(P_w_cond Jn+P_w__cond_out)/2  { psia } 

{6.  Condenser  flow  rate} 

GPM_w_cond=3000  {gpni} 

m_dot_w_cond=rho_w_cond*GPM_w_cond*con  vert(gpm,ft'^3/hr)  { Ibm/hr } 


{Section  VII.  Condenser  pump  model } 

(A.  Inputs  for  condenser  pump} 

h_pump_condJn=h_w_tow_out  {btu/lbm} 

P_pump_cond_in=P_atm-HeadJo  wtopump*rho_w_pump_cond*g*convert(lbm/ft*s'^2,psia)  { psia } 
Head_to  wtopump=(Head_pipe_cond/2“Elev_to  w)  { ft } 

m_dot_  w_pump_cond=m_dot_w_cond  { Ibm/hr } 

{B.  Condenser  pump} 

eta_pump_cond=.83 

eta_motor_cond=.90 

eta_pump_cond=(h_pump_cond_out_rev-h_pump„cond_in)/(h_pump_cond_out-h_pump_cond_in) 

W_dot_pump_cond=m_dot_w_pump_cond*(h_pump_cond_ouFh_pump_cond_in)/eta_motor_cond 


KW_pump_cond=W_dot_pump_cond*con  vert(Btu/hr,kW)  { kW } 

{C.  Pressure  calculation} 

DELT  AP_pump_cond=rho_w_pump_cond*g*Head_loop_cond*convert(lbm/ft*S^2,psia)  { psia } 

P_pump_cond_out=P_pump_cond_in+DELT  AP_pump_cond  { psia } 

P_pump_cond_a  ve=P_pump_cond_in+DELT  AP_pump_cond/2  { psia } 

Head_pipe_cond=f Joop_cond*Length_loop_cond/D  J_pipe_cond*(FPS_pump_cond^2/2)/ g  { ft } 


{eqn630ref[5]} 

Head_loop_cond=Head_pipe_cond+Head_cond+Head_str  ainer  { ft } 

Head_strainer=30  {ft} 

fJoop_cond=l/(-1.8*loglO(6.9/Re_D_Lpump_cond+(epsilon_pipe_cond/(D_Lpipe_cond*3.7))'^l.ll))'^2 

{eqn  6.64a  ref[5]} 

{D.  Pipe  geometry} 

D_Lpipe_cond=i2*convert(in,ft)  {ft} 

Length Joop_cond=500  {ft} 

epsilon_pipe_cond=.000 15  { ft } 

Elev_tow=200  {ft} 
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{E,  Flowrates} 

GPM_pump_cond=m_dot_w_pump_cond/rho_w_pump_cond*convert(ft^3/hr,gpm)  { gpm } 


CFM_puinp__cond=GPM_pump_cond*convert(gpm,cfin)  { ft^3/min } 

FPM_pump_cond=CFM_pump_cond/(pi*D_i_pipe_cond^2/4)  { ft/min } 

FPS_pump_cond=FPM_pump_cond*convert(  1/min,  1/s)  { ft/s } 

FPH_pump_cond=FPM_pump_cond*convert(  1/min,  1/hr)  { ft/hr } 


Re_D_Lpump_cond=rho_w_pump_cond*FPH_pump_cond*D_i_pipe_cond/mu_w_pump_cond 
{F.  Fluid  Froperties} 

rho_W_pump_cond=density(WATER,T=T_pump_cond_ave,P=P_pump_cond_ave)  { lbm/ft'^3 } 

niu_W_pump_cond= viscosity  (WATER, T =T_pump_cond_ave  J^=P_pump_cond_ave)  { Ibm/ft  hr } 

{G.  Temperatures} 

T_pump_cond  Jn=temperature(WATER,h=h_pump_cond Jn,P=P_pump_cond_ave)  { F } 

T_pump_cond_ave=(T_pump_cond_out-T_pump„cond„in)/ln((T_pump_cond_out)/(T_pump_cond_in)) 
T_pump_cond_out_rev=temperature(WATER,P=P_pump_cond_out,s=s_pump_cond_in)  { F } 
T_pump„cond_out=temperature(W  ATER  J^=P_pump_cond_out,H=h_pump_cond_out)  { F } 

{H.  Enthalpy  and  entropy} 

h_pump_cond_out_rev=enthalpy(W  ATER,P=P_pump_cond_out,T=T_pump_cond_out_rev)  { btu/lbm } 
s_pump_cond Jn=entropy(W  ATER,P=P_pump_cond Jn,T=T_pump_cond Jn)  { btu/lbm  R } 

s_pump_cond_out=entropy  (W  ATER  J^=P_pump_cond_out,T =T_pump_cond_out)  { btu/l  bm  R } 

{Section  VIII.  Cooling  tower  model} 

(A.  Inputs  for  cooling  tower} 
h_w__tow__in=h_w_cond_out  {btu/lbm} 

m_dot__w_tow=m_dot_w__cond  {Ibm/hr} 

{B.  Cooling  tower  calculations} 

Call  to werfan(Q_dot_tow,epsilon_to w,m_dot__a_to w,h_s_w_to w_in,h_a_to w_in,T_w_tow_out  1 : 
fanduty,T_w_tow_out) 

Q_dot_tow=0.9*fanduty*epsilon_tow*m_dot„a_tow*(h_s_w_towjn-h_a_tow__in)+ 

0.  l*epsilon_tow*m_dot_a_tow*(h_s_w_tow Jn-h_a_towJn)  { Btu/hr } 

Q_dot_to  w=m_dot_w_tow*(h_w_tow_in-h_w_tow_out)  { btu/hr } 

epsilon_tow=(  1  -exp(-NTU_tow*(  l-Cr_tow)))/(  l-Cr__tow*exp(-NTU_tow*(  l-Cr_tow))) 

{eqnl  1.30a  ref[3]} 

NTU_tow=1.3936*(m_dot_w_tow/m_dot_a_tow)^(l-.7044)  {eqn  11.17  ref[9] } 

Cr_tow=(m_dot_a_tow*Cs)/(m_dot_w_tow*Cp_w_tow)  { eqn  II.  12  ref[9] } 

Cs=(h_s_w_to w_in-h_s_w_tow_out)/(T_w_to w_in-T_w__tow_out)  { btu/lbm  R }  { eqn  11. 1 3  ref[9] } 

{C.  Water  properties} 

rho_w_tow=density(WATER,T=T_w_tow_ave,P=P„atm) 
Cp_w_tow=specheat(WATER,T=T_w_tow_ave,P=P_atni) 
h_s_w_to  w Jn=enthalpy  ( AIRH20,R=  1  ,T =T_w_to  w_in,P=P_atm) 
h_s_w_tow_out=enthalpy(AIRH20,R=l,T=T_w_tow_out,P=P_atm) 

{D.  Air  Properties} 

rho_a_tow=:density(  AIRH20,B=ODWB,T=T_a_tow_in,P=P_atm)  { lbm/ft'^3 } 

Cp__a_to  w=specheat(  AIRH20,B=0DWB  ,T=T_a_to  w_in,P=P_atm)  {Btu/lbm  R } 

h_a_tow_in=enthalpy(AIRH20,B=ODWB,T=T_a_towJn,P=P_atm)  {Btu/lbm} 


{lbni/ft^3} 
{Btu/lbm  R} 
{Btu/lbm} 
{Btu/lbm} 


{E.  Temperatures} 

T_a_towJn=T_amb  {F} 

T_w_to  w_in=temperature(W  ATER,h=h_w_tow_in,P=P_atm)  { F } 

T_w_to  w_out  l=temperature(WATER,h=h_w_tow_out,P=P_atm)  { F } 

T_w_to w_ave=(T_w_tow_in-T_w_tow_out)/ln((T_w_to w_in)/(T_w_to w__out))  { F } 

{F,  Entropy} 

s_w_tow Jn=entropy(W  ATER,h=h_w_tow Jn,P=P_atm)  { btu/lbm  R } 

s_w_tow_out=entropy(WATER,h=h_w_tow_out  J=P_atm)  { btu/lbm  R } 

{G.  Flowrates} 

ni_dot„a_tow=rho_a_tow*cfm_to  w*convert(  1/min,  l/hr)  { Ibm/hr } 

cfm_to  w=250*  1 000 

{ H.  Cooling  tower  fan  model } 

FHPCA=-2.3084  {Appendix  A  ref[  10]} 

FHPCB=6.3769  { Appendix  A  ref[  10] } 

W_dot_fan_tow=(exp(FHPCA+((ln(cfm_tow)-i0.5)/2)*FHPCB)/.9)*fanduty*convert(hp,btu/hr) 

{ btu/hr }  ( Appendix  A  ref[  1 0] } 

KW_fan_tow=W_dot_fan_tow*convert(btu/hr,kW)  { kW } 

{Section  IX.  Constants} 

P_atm=14.7 
g=32.2 
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